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INTRODUCTION. 


I  am  of  the  opinion  that  the  present  work,  entitled  "The  Harmonic  Formula  of 
Fourier  and  Bessel  and  Its  Application  to  the  Study  of  the  Diurnal  Variation  of  the  At- 
mospheric Pressure  in  Manila  during  the  Period  1890-1909,"  prepared  by  Rev.  Antonio 
Galan,  S.  J.,  a  temporary  member  of  the  staff  of  the  Manila  Observatory,  will  be  very 
well  received  by  those  meteorologists  who  think  that  such  studies  and  discussions  should 
be  made  in  all  the  principal  observatories  where  sufficient  material  has  been  collected 
for  a  number  of  years. 

I  hope  that  in  a  short  time  we  may  be  able  to  find  time  and  opportunity  to  work 
up,  in  a  manner  similar  to  that  done  by  Father  Galan  for  the  barometric  data,  the 
observations  collected  in  the  Observatory  on  the  temperature,  humidity,  and  other  pe- 
riodic elements  of  meteorology. 

And  I  hope  also  that  the  stay  in  Manila  of  Father  Galan,  who  is  appointed  to  con- 
tinue the  splendid  work  of  Fathers  Vines,  Gangoiti,  and  Gutierrez  Lanza  of  the  Belen 
Observatory,  Habana,  may  be  as  useful  to  him  as  my  stay  in  Habana  during  the  first 
few  months  of  1892  was  to  me. 

JOSE  ALGUE, 
Director,  Weather  Bureau. 

MANILA  OBSERVATORY,  June  19,  1914. 
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THE  HARMONIC  FORMULA  OF  FOURIER  AND  BESSEL. 


The  harmonic  formula  of  Fourier  and  Bessel  is  usually  presented  in  one  or  other 
of  the  two  following  forms: 

(I)  y=M+p1  cos  X  +  Q!  sin  x+p2  cos  2x  +  q2  sin  2x+ps  cos  3*+g8  sin  3«-f-  .... 

(II)  y  =  M+UI  sin  ( U^+x)  +u2  sin  ( U2+2x)  +ua  sin  ( Z73+3a;)  + 

It  is  evident  that  these  two  formulae  are  one  and  the  same  under  different  form,  for 
developing  (II)  we  have: 

y  =  M-\-u1&m  Ul  cos  x-\-ul  cos  U:  sin  x-\-uz  sin  U2  cos  2x-\-uz  cos  U2  sin  2x-\- etc. 

and  making  in  the  expression— 

ui  sin  Ui=pi  Ui  cos  U\=qi 

Uz  sin  Uz=p2  Uz  cos  Uz=qz 

etc.  etc. 

we  find  that   (II)   is  reduced  to   (I). 

The  practical  applications  of  this  formula  are  many  and  important. 

Fourier's  theorem  (say  Thomson  and  Jait)  is  not  only  one  of  the  most  beautiful  results  of  modern 
analysis,  but  may  be  said  to  furnish  an  indispensable  instrument  in  the  treatment  of  nearly  every 
recondite  question  in  modern  physics.  To  mention  only  sonorous  vibrations,  the  propagation  of  electric 
signals  along  the  telegraph  wire,  and  the  conduction  of  heat  by  the  earth's  crust,  as  subjects  in  their 
generality  intractable  without  it,  is  to  give  but  a  feeble  idea  of  its  importance. 

In  the  present  study  we  limit  ourselves  to  the  application  which  is  made  of  the  for- 
mula in  meteorology;  and  we  will  endeavor  to  show  what  each  one  of  its  elements  repre- 
sents, and  point  out  its  utility  and  practical  application.1 

The  greater  number  of  the  phenomena  studied  in  meteorology  are  periodic,  and 
the  curves  which  represent  them  are  consequently  periodic  curves — that  is  to  say, 
curves  which  represent  the  movement  of  a  point  which  passes  an  indefinite  number 
of  times  over  the  same  path  in  the  same  time. 

Fourier  proved  that  any  periodic  curve  could  be  resolved  into  a  series  of  harmonic 
curves  of  periods  1,  £,  £,  etc.,  of  the  given  curve,  and  that  only  one  combination  of  these 
elementary  curves  was  possible  to  reproduce  a  specified  curve.  This  corresponds  with 
the  fact  observed  by  Helmholtz  that  the  same  composite  sound  is  always  resolved  into 
the  same  elementary  sounds. 

Prescinding  for  the  moment  from  the  practical  utility  which  this  resolution  of 
curves  may  have  in  meteorology,  let  us  examine  how  the  formula  which  we  are  study- 
ing is  nothing  more  than  the  analytical  expression  of  periodic  curves  in  terms  of  their 
harmonic  components. 

A  harmonic  curve  is  one  which  represents  graphically  the  harmonic  motion  of  a 
point. 


1  One  of  the  best  practical  treatises  which  we  know  on  the  employment  of  Fourier's  series  in  mathe- 
matical physics  is  "An  Elementary  Treatise  on  Fourier's  Series  and  Spherical  Harmonics"  by  W.  E. 
Byerly,  published  by  Ginn  &  Co.,  New  York. 

The  original  work  of  Fourier  himself  "Theorie  Analytique  de  la  Chaleur"  is  perhaps  as  up  to 
date  as  any  other  modern  publication  on  this  subject.  Mellor's  "Higher  Mathematics  for  Students  of 
Chemistry  and  Physics,"  Longmans,  London  and  New  York,  may  also  be  consulted. 
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HARMONIC   FORMULA   OF   FOURIER  AND  BESSEL. 


Fig  1 


Harmonic  motion  is  that  to  which  the  projection  of  a 
point,  moving  uniformly  around  a  circumference  on  to  a 
diameter,  is  subjected.  The  circumference  is  called  the 
circle  of  reference. 

If  the  point  M  (fig.  1)  moves  around  the  circumference 
uniformly,  the  motion  of  its  pro- 
jection P  on  the  diameter  BD  is 
harmonic;  while  M  describes  the 
circumference  once  in  the  direc- 
tion ABCDA,the  point  P  traverses 
the  diameter  in  the  direction 
OBODO. 

As  a  help  to  those  who  may 


not  be  familiar  with  this  class  of  motion,  the  following  defini- 
tions of  terms  are  given : 

Period  is  the  time  required  for  the  point  M  to  describe 
the  circumference,  or  for  P  in  its  double  journey  of  the 
diameter. 

Elongation  is  the  variable  distance  between  the  center  O 
and  the  moving  point  P. 

Amplitude  is  the  maximum  elongation,  or  maximum 
distance  which  the  point  P  can  be  separated  from  0;  as  is 
evident,  this  distance  is  the  radius  of  the  circle  of  reference. 

O  and  A  are  generally  taken  as  the  starting  points  of  the 
motion.  The  zero  of  time  is  arbitrary;  it  may  be  taken  on 
any  point  in  the  circumference. 

Phase  is  the  angular  distance  between  the  moving  point 
and  the  starting  place  of  the  motion.  In  speaking  of  the 
point  which  moves  on  the  diameter,  it  is  necessary  to  refer 
it  to  the  circle  of  reference  in  order  to  determine  its  phases; 
for  example,  in  figure  1  if  the  point  is  found  in  P,  its  phase 
is  described  by  the  degrees  of  the  arc  MA  or  of  the  angle  MO  A. 

Epoch  of  a  definite  phase  is  the  instant  at  which  this 
phase  takes  place;  v.  g.,  if  the  arc  MA  is  60°  and  it  is  two 
hours  since  the  moving  point  was  at  the  zero  of  time,  we  say 
that  2  o'clock  is  the  epoch  of  the  phase  of  60°. 

The  initial  phase  is  the  angular  distance  between  the  zero 
of  time  and  the  starting  place  of  the  motion.  In  the  figure, 
if  we  take  the  point  N  as  the  zero  of  time,  the  initial  phase  of 
the  motion  is  given  by  the  arc  NA. 

Suppose  that  the  point  M  requires  a  day  to  describe  the 
circumference — that  is  to  say,  that  the  period  of  its  motion 
is  twenty-four  hours.  Starting  from  the  zero  of  time  N  (fig.  2) , 
let  us  divide  the  circumference  into  24  equal  parts ;  draw  the 
straight  line  NX  equal  to  the  length  of  the  circumference; 
this  line  will  also  be  divided  into  24  equal  parts,  each  one 
of  which  represents. the  length  of  an  arc  of  15°,  or  one  hour 
of  time.  Suppose  the  circle  of  reference  to  move  along  the 
straight  line  NX  parallel  to  the  vertical  diameter  with  a  motion 
equal  to  that  of  the  moving  point  in  the  circumference;  at  the 
end  of  each  hour  the  center  0  will  pass  through  each  of  the 
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^successive  points  in  which  the  straight  line  was  divided,  and  the  point  which  moves 
on  the  diameter  will  describe  the  sine-curve  of  the  figure. 

If  it  is  desired  to  describe  another  harmonic  curve  of  a  period  i  of  the  former,  the 
parts  into  which  the  circle  of  reference  is  to  be  divided  must  contain  double  the  num- 
ber of  degrees,  for  it  is  clear  that,  if  the  circumference  is  to  be  described  in  twelve 
hours,  the  point  must  pass  through  30°  each  hour.  In  like  manner  for  curves  of 

periods  i,  i,   of  the  first,  the  respective  circles  of  reference  must  be  divided 

into  parts  containing  triple,  quadruple,  etc.,  the  number  of  degrees;  so  that  the  periods 
of  these  curves  are  inversely  to  the  number  of  degrees  contained  by  the  parts  into 
which  their  circles  of  reference  are  divided. 

Any  ordinate  of  the  curve  represents  the  elongation  (distance  of  the  moving  point 
P  from  O)  in  the  corresponding  moment.  Call  these  ordinates  ?//,  y/',  etc.  Then 
in  the  rectangular  triangle  POM  (fig.  1)  we  have  as  a  general  expression  of  the 
ordinates — 

OP=y=OM  sin  a=u  sin   (U+x) 

where  u  is  the  radius  of  the  circle  of  reference,  or  the  amplitude  of  the  curve,  U  the 
initial  p~hase,  and  x  the  parts  of  the  circumference  through  which  the  point  M  has 
moved  at  a  given  moment. 

With  these  considerations  in  view,  we  get  for  the  different  ordinates  of  the  char- 
acteristic harmonic  curves  of  any  periodic  curve  the  following  expressions: 

Fiist    ordinal  e.  Second   ordinate.  Third  ordinal e. 

First  curve  ....  ul  sin  17,  M1sin(£71+  x)  ulsm(Ul+2x) 

Second  curve  ..  u.,  sin  £7,  u2  sin  (U 2-\-2x)  u2  sin(Z72-(-4a:) 

Third  curve  ....  w..  sin  Ua  u.A  sin(Z78-f- 3x)  ua  sin(l73+6o;) 


Making  an  algebraic  sum  of  homonymous  ordinates  we  have  the  following  expres- 
sions for  the  corresponding  ordinates  of  the  resultant  periodic  curve  : 

Resultant. 


First  ordinate  =  ul  sin  U^+u.,  sin  U.2-\-u..sin  U.t+ 

Second  ordinate  =  u^  sin  (  U^x)  +u2  sin  (  U.,-\-2x)  +u.,  sin  (  U3-\-3x)  + 

Third  ordinate  =  u^  sin(f7,-f2a:)  -\-u.,  sm(t72  +  4#)  -(-w3  sin([7:i+6a:)  -(- 


which,  as  is  readily  seen,  are  nothing  more  than  particular  expressions  of  the  general 
formula — 

«,    sin(?71+a;)+M,   sm(U.z+2x)  +u,   sin  ( £7, + So:)  -f '. 

in  which  the  variable  element  x  successively  takes  the  corresponding  values. 

Deducing  Bessel's  formula  by  this  process,  what  was  said  in  the  beginning  becomes 
clear,  viz,  that  it  is  the  analytical  expression  of  a  periodic  curve  in  terms  of  its  harmonic 
components ;  each  one  of  the  terms  representing  one  of  these  components. 

The  constants  u  and  U  are  usually  called  harmonic  coefficients  (u  lineal,  and  U 
angular) ,  for  they  are  the  elements  which  characterize  the  harmonic  motion  in  question ; 
K  gives  the  amplitude,  and  consequently  the  importance  of  the  motion,  for  as  may  be 
seen  in  the  figures,  according  as  u  decreases,  the  motion  of  the  point  P  on  the  diameter 
becomes  slower;  U  (fig.  2)  fixes  the  epoch  of  the  different  phases;  an  increase  in  its 
value  accelerates  the  epoch  of  all  the  phases,  and,  on  the  contrary,  a  decrease  retards  it. 

It  is  of  great  importance  to  remember  this  property  of  the  initial  phase;  for 
example,  if  in  a  station  whose  longitude  was  15°  W.  of  Greenwich  it  was  required  to 
calculate  in  Greenwich  mean  time  instead  of  in  local  mean  time,  the  zero  of  time  is 
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advanced  one  hour,  and  it  is  hence  necessary  to  diminish  the  value  of  the  initial  phases 
15°  in  Ulf  30°  in  U2,  45°  in  £78,  etc.  Hence,  it  is  important  that  in  all  cases  the  time 
employed  should  be  explicitly  stated. 

The  quantity  M  which  appears  in  the  formula  represents  the  mean  of  the  observa- 
tions which  enter  into  the  discussion.  What  is  developed  in  the  trigonometrical  series 
of  the  formula  is  the  difference  between  each  observation  and  the  mean  of  all  of  them ; 
consequently,  in  order  that  the  observations  (or  the  ordinates  of  the  observation  curves) 
be  represented  according  to  this  formula,  it  is  necessary  to  add  the  quantity  M  to  the 
development  of  these  differences. 

What  has  been  said  is  no  doubt  sufficient  to  give  a  general  idea  of  the  mechanism 
of  the  formula,  and  in  the  application  which  will  be  made  of  it  in  the  study  of  the 
diurnal  variation  of  the  atmospheric  pressure  in  Manila  we  will  endeavor  to  make  this 
idea  more  concrete. 

APPLICATION  OF  THE  FORMULA  TO  THE  ANALYTICAL  STUDY  OF  THE  DIURNAL 

VARIATION  OF  THE  ATMOSPHERIC  PRESSURE  IN  MANILA 

DURING  THE  PERIOD   1890-1909. 

In  the  preparation  of  the  data  which  were  to  enter  directly  in  the  calculation  the 
following  method  was  followed : 

The  hourly  means  of  each  month  were  taken  from  the  official  publications  of  the 
Observatory,  and  the  mean  of  the  zero  hour,  or  first  midnight,  of  the  representative 
day  of  each  month  was  calculated  separately;  this  mean  is  the  mean  value  of  the  read- 
ings of  the  first  midnight  of  each  day  of  the  month.  As  the  zero  hour,  or  first  mid- 
night, of  a  day  is  the  last  midnight  of  the  preceding,  in  any  month  (January  for 
example)  the  readings  of  the  Oh  will  be:  For  the  1st,  that  of  the  12  p.  m.  of  the  last 
day  of  December  of  the  preceding  year ;  for  the  2d,  that  of  12  p.  m.  of  January  1,  etc., 
so  that  the  values  that  enter  into  the  mean  of  this  hour  are  the  same  as  those  of  12 
p.  m.,  except  the  reading  for  12  p.  m.  of  the  last  day  of  each  month,  which  is  substituted 
for  the  corresponding  one  of  the  previous  month. 

It  is  evident,  therefore,  that,  having  calculated  the  mean  of  the  12  p.  m.  readings, 
it  is  sufficient,  in  order  to  calculate  that  of  the  Oh  readings,  to  subtract  the  12  p.  m. 
reading  of  the  last  day  of  the  month  in  question  from  the  corresponding  reading  of 
the  preceding  month,  divide  this  difference  by  the  number  of  days  of  the  month,  and 
add  the  result  algebraically  to  the  mean  of  the  12  p.  m.  readings. 

Having  obtained  these  25  hourly  means,  the  general  hourly  means  of  the  twelve 
months  for  the  whole  period  have  been  calculated  by  grouping  together  like  months, 
and  taking  the  mean  of  the  means  obtained  for  each  hour  of  a  given  month  in  different 
years.  These  general  hourly  means  have  been  approximated  to  thousandths  and  are 
given  in  Table  I. 

Table  II  is  the  same  as  I  with  the  aperiodic  element  removed.  The  following  will 
give  an  idea  of  this  correction : 

Noncyclic  or  aperiodic  variation. — In  the  explanation  of  the  formula  we  saw  that 
the  elementary  harmonic  curves  consist  of  25  ordinates  corresponding  to  the  different 
hours  of  the  day,  from  Oh  to  24h,  both  inclusive.  The  extreme  ordinates  are  equal  for 
the  same  curve,  consequently  the  ordinates  corresponding  to  Oh  and  24h  in  the  resultant 
curve  must  be  equal,  in  that  they  are  the  algebraic  sum  of  the  corresponding  ordinates 
in  different  elementary  curves. 

Hence,  it  will  be  understood  why  it  is  necessary  to  calculate  the  mean  of  the  zero 
hours;  for  one  must  be  quite  sure  that  the  curve  to  be  analyzed  as  a  periodic  one  has 
the  essential  characteristic  of  periodicity,  viz,  that  its  extreme  ordinates  be  equal.  In 
case  this  is  not  so,  then  some  element  has  been  introduced  into  the  march  of  the  phenom- 
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enon  which  perturbs  its  periodicity.  This  perturbation  is  what  is  called  noncyclic  or 
aperiodic  variation. 

What  is  ordinarily  the  cause  and  origin  of  this  variation  will  be  understood  by  a 
familiar  illustration. 

The  temperature  varies  in  the  different  seasons  of  the  year,  rising  in  spring  and 
summer  and  falling  in  autumn  and  winter;  so  that  in  May,  for  example,  although  it 
may  not  be  said  in  general  that  every  day  is  warmer  than  the  preceding,  nevertheless 
the  general  tendency  is  toward  a  rise  of  temperature;  and  if  we  wish  to  calculate  the 
mean  diurnal  variation  of  temperature  for  this  month  in  an  average  year  by  assigning 
to  each  hour  of  the  day  the  mean  for  that  hour  derived  from  the  whole  thirty-one  days, 
we  shall  find  that  the  second  midnight  of  the  representative  day  is  somewhat  warmer 
than  the  first  midnight. 

In  the  case  of  barometric  observations  we  come  across  an  analogous  phenomenon. 
It  is  well  known  that  the  atmospheric  pressure  is  subject  to  an  annual  oscillation  over 
and  above  the  diurnal  oscillation.  In  studying,  therefore,  the  characteristics  proper  to 
the  monthly  diurnal  variation  of  atmospheric  pressure,  we  must  begin  by  eliminating 
the  aperiodic  element. 

This  elimination  is  usually  performed  by  applying  to  each  hour  the  correction — 

D(12-h) 
24 

in  which  D  is  obtained  by  subtracting  the  value  of  the  Oh  from  that  of  the  24h,  and  h 
represents  the  hour  in  which  the  correction  is  to  be  applied.  This  correction  is  logical 
only  on  the  hypothesis  that  D  comes  in  at  a  uniform  rate  throughout  the  twenty-four 
hours. 

To  prove  the  legitimacy  of  this  hypothesis  presents  difficulties ;  nevertheless,  when 
the  investigation  is  limited  to  the  twenty-four  hours  it  is  usually  accepted,  for  it  seems  to 
be  the  only  workable  hypothesis. 

Let  us  see  how  the  formula  implicitly  supposes  this  hypothesis.  If  the  variation  is 
supposed  to  be  uniform  during  the  twenty-four  hours,  it  must  be  represented  by  a 
straight  line.  The  correction  formula  will,  therefore,  be  the  equation  of  the  straight 
line — 

y=b+ax(l) 

Take  as  origin  of  rectangular  coordinates  the  zero  hour  and  transfer  as  usual  the 
hours  to  the  axis  X.  In  the  formula  (1)  x  represents  the  successive  hours ;  y  the  cor- 
rection to  be  applied  to  each.  When  the  correction  has  been  applied  to  the  extreme 
hours  their  values  must  coincide;  therefore,  the  semldifference  of  their  values  is  to  be 
added  to  one  and  subtracted  from  the  other.  The  correction  for  these  two  hours,  called 
D  their  difference,  will  be — 

for  the  24h__     ±^= 

£* 

fortheOh_.  -+§=* 
Eliminating  b,  we  have  n 

±D=24a.-.a=±~. 

and  substituting  the  values  obtained  for  a  and  b  in  the  equation   (1)   we  get — 

-D  .D      D  (12- x) 

*)J -1 ~ (— rf i . •— 

2  ~24X~          24 

Hence,  we  see  that  the  hypothesis  which  requires  the  aperiodic  variation  to  be  uniform 
during  the  twenty-four  hours  is  necessary  and  sufficient  for  the  deduction  of  the  correc- 
tion formula. 
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An  example  will  demonstrate  the  influence  which  this  correction  has  on  the  curves. 

Let  us  take  the  month  of  November,  which  in  our  case  has  the  largest  correction. 
The  smooth  curve  (fig.  3)  represents  the  variation  without  correction;  the  dotted  one 
represents  the  same  variation,  but  with  the  correction  applied. 

From  the  form  and  position  of  these  curves  it  is  clearly  seen  that  the  elimination 
of  the  aperiodic  element  is  reduced  to  revolving  the  curve  through  a  small  angle  with 
12  a.  m.  as  the  center,  the  characteristic  undulations  remaining  the  same. 

Once  this  correction  applied,  the  diurnal  means  were  calculated  for  Table  II  by 
dividing  the  24  hourly  means  of  each  month,  from  Oh  to  23,h  by  24,  and  Table  III  was 
formed  by  subtracting  the  diurnal  mean  of  each  month  from  each  hourly  mean. 

In  this  last  calculation  we  have  approximated  only  up  to  two  decimal  places;  but, 
as  in  the  previous  calculations,  the  approximation  was  carried  to  three  decimal  places, 
the  inaccuracy  of  the  data  in  the  last  table  cannot  be  greater  than  0.005  and  the  sum  of 
the  positive  values  does  not  differ  absolutely  from  that  of  the  negative  by  more  than  two 
or  three  hundredths  at  most,  which  would  not  have  been  obtained  if  in  the  earlier  cal- 
culations the  approximation  had  only  been  carried  to  two  decimal  places. 

In  the  analytical  study  which  follows  we  have  used  local  mean  time. 

All  that  is  necessary  for  the  determination  of  each  of  the  elementary  waves,  into 
which  the  diurnal  barometric  variation  is  to  be  resolved,  are  the  lineal  harmonic  coeffi- 
cients u,  and  the  angular  harmonic  coefficients  U — that  is  to  say,  the  constants  of  For- 
mula II,  for  these  are  the  ones  which  characterize  those  waves. 


Fig.  3 

These  constants  cannot  be  found  directly;  they  must  be  determined  by  means  of 
the  constants  p  and  q  of  Formula  I  which  are  calculated  directly. 
The  equations  of  these  constants  are,  as  we  saw  before: 

p=u  sin  U\  /-,\ 
q=u  cos  U)  ^ 
Hence,  we  obtain  the  relations: 

sin  U  p      .  p  q 

-r-.^tan   U=    ;  sm  U=    ;  cos  U= 
cos  U  q  u  u 

From  the  first  of  these  equations  we  can  deduce  the  angular  constant  U  in  terms  of 
p  and  q;  and  from  the  last  two  in  terms  of  p  or  q  and  u;  this  last  is  obtained  by  the 
equation  : 


which  is  easily  obtained  from  the  equation   (1)   by  squaring  and  adding.     All  the  diffi- 
culty consists  in  determining  the  constants  p  and  q. 

The  method  which  Bessel  elaborated  for  this,  and  which  we  will  explain  at  the  end 
of  this  study,  is  long,  and  tedious,  and,  moreover,  requires  a  constant  concentration  in 
the  work  if  errors  are  to  be  avoided;  this  is  owing  to  the  great  number  of  angles  to 
be  reduced,  and  the  frequency  with  which  the  signs  change  in  passing  from  one  qua- 
drant to  another.  Father  Dechevrens,  S.  J.,  says  of  this  method:  "La  methode  na- 
turelle  n'est  que  pour  les  habiles."  This  Father  has  published  two  pamphlets,  one  en- 
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titled  "Methode  simplifiee  dite  des  facteurs  pour  le  calcul  des  series  de  Fourier  et  de 
Bessel  appliquees  a  la  Meteorologie,"  in  which  by  means  of  tables  he  simplifies  consider- 
ably the  calculation.  The  other  pamphlet  is  entitled  "Note  complementaire  a  la  methode 
simplifiee  du  calcul,  etc."  l  in  which  he  explains  the  manner  in  which  the  tables  are 
formed. 

The  first  three  of  these  tables  serve  for  the  determination  of  the  coefficients  p  and 
q  in  series  of  4,  6,  8,  10,  12,  24,  and  36  observations. 

We  think  that  it.  will  not  be  without  interest  for  those  who  may  not  be  accustomed 
to  this  class  of  analysis,  to  work  out  in  detail  an  example  by  Bessel's  method,  and  see 
in  what  the  utility  of  Father  Dechevrens'  tables  consists  ;  and,  even  at  the  risk  of  being 
too  detailed,  we  will  follow  step  by  step  the  method  of  Bessel  in  the  determination  of  the 
constants  p  and  q  for  January. 

The  equations  which  give  for  the  present  case  the  values  of  these  constants  for  the 
first  eight  terms  of  Formula  (I)  are:2 


i=am  cos  ra#  qi=am  sn  race 

='%am  cosm(2x)  I2q2=2>am  sin  m(2x) 

—  2am  cos  ra  (3x)  I2qs=^<am  sin  ra  (3x) 

=3am  cos  ra  (4#)  12<74—  Sam  sin  ra  (4#) 

The  sign  2  indicates  that  the  quantity  which  it  affects  is  to  be  developed  into  a 
series  of  terms  of  the  same  form,  which  follow  according  to  the  series  of  natural  num- 
bers 0,  1,  2,  etc.,  so  that  we  write: 

Sa,u  cosmx=a,o  cos  0+a\  cos  x-i-az  cos  2x+a3  cos  3#  +  ____  :  ____ 

x  represents  the  value  of  one  of  the  parts  into  which  the  circle  of  reference  of  the  first 
elemental  curve  is  divided.  As  the  period  of  this  curve  is  twenty-four  hours,  its  circle 
of  reference  is  divided  into  24  parts,  and  consequently  #=15°.  The  periods  of  the  other 
curves  are  ±,  ^,  %  of  that  of  the  first,  and  therefore  the  divisions  of  their  respective 
circles  of  reference  equal 

2z;3x;4x;  or  30°;  45°;  60°; 

a,,,  represents  the  values  of  Table  III  corresponding  to  the  different  hours;  thus  for 
January  — 

a-o—46,  a-!—  17,  Oz=—  20,  etc. 

According  to  this,  in  order  to  determine  p,  in  January  we  have  — 
I2pi~=  46  cos  0°-h  17  cos  15°  -20  cos  30°  -45  cos  45°.  _ 

All  that  is  then  required  is  to  find  the  values  of  these  cosines  in  the  trigonometrical 
tables  and  work  out  the  operation.     The  work,  which  would  be  enormous  if  the  multi- 
plication had  to  be  made  in  each  term,  can  be  simplified  in  the  following  manner. 
The  value  of  the  cosines  in  the  different  hours  are: 


0"  ... 

cos  0° 

—   1 

12"  

—  _i 

1  

cos  15° 

=   .966 

13  

--  .966 

2  

cos  30° 

=   .866 

14  

--  .866 

3  

cos  45° 

=   .707  j 

15  

--  .707 

4  

cos  60° 

=   .5 

16  

--  .5 

5  

cos  75° 

=   .259  ; 

17  

--  .259 

6  

cos  90° 

A 

18  

-  0 

7  

....  cos  105°  —  —cos  75° 

=  -  .259 

19  

-   .259 

8  

cos  120°  —  —cos  60° 

=  -  .5 

20  

—    5 

9  

=  -  707 

21  

-   .707 

10  

=  -  .866 

22  

-   .866 

11 

=  -  .966 

23  

-   .966 

1  Tipografia  della  Pace  di  Filippo  Cuggiani,  Roma,  Via  della  Pace,  num.  35. 

2  The  mathematical  proof  of  these  equations  will  be  given  later  on. 
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Each  of  those  numbers  appears  as  a  common  factor  of  several  hours, 
therefore : 


We  have 


12pi==(0-12)l  +  (1-11-13+23).  966+  (2-10— 14+22).  866  H- (3-9-15+21).  707 

+  (4-8-16+20). 5+ (5-7— 17+19). 259 

Substituting  for  the  hours  which  appear  in  parentheses  the  corresponding  values  taken 
from  the  Table  III  we  get: 

12pi=(46-15)l+  (17-79+63  +  59).  966+_ 
In  the  determination  of  q}  we  obtain  by  the  same  method: 

0"  sin  0°zzO  12"  =     0 

1  , sin  15° zz   .259     13  =  -   .259 

2  sin  30°  zz   .5  14  zz  -   .5 

3  sin  45°  zz   .707  15  —  —   .707 

4 sin  60°  zz   .866  j  16  zz  —   .866 

5 sin  75°  zz   .966  i  17  zz  —   .966 

6  sin  90°zzl  18  zz  — 1 

7  sin  105°zzsin  75°  =   .966  |  19  =  -   .966 

8  =   .866  i  20  zz—   .866 

9  zz   .707  |  21  zz-   .707 

10  zz   .5         22  zz-   .5 

11  zz   .259    23  zz-   .259 

Taking  into  account  the  common  factor  of  several  hours  we  write: 

12gi=(l+ll-13-23).259  (-(2+10-14— 22). 5+  (3+9— 15— 21). 707 

f  (4+8-16-20).  866  +  (5-f-  7— 17-19).  966-;-  (6— 18)1 

and  substituting  their  values  for  the  hours  we  find  gv 

For  the  determination  of  p2  and  q2  we  proceed  in  the  same  way,  remembering  that 
the  divisions  of  the  circle  of  reference  equal  30°,  we  say  therefore: 


P2 


0"  cos  Oc 

1  cos  30C 

2  cos  60C 

3  cos  90C 

4  cos  120°  zz— cos  60C 

5  cos  150°  zz -cos  30C 

6  

7  

8  

9  

10  

11  . 


1 

12" 

—  1 

.866 

13 

—   .866 

.5 

14 

—   .5 

0 

15 

—  0 

—  .5 

16 

—  —  .5 

—  .866 

17 

—  —  .866 

—  1 

18 

—  .866 

19 

—  —  .866 

—  .5 

20 

—  -  .5 

0 

21 

—  0 

.5 

22 

—   .5 

.866 

23  .. 

..  -   .866 

12p2=(0-6+12-18)l+(l-5— 7+11+13-17— 19 


23). 866+ (2-4-8 

+10+14- 


16-20+22).5 


sin  0' 

sin  30C 

sin  60' 

sin  90C 

4  sin  120°zzsin  60C 

5  

6  

7  

8 

9  

10  

11  .. 


zz       0 


=     1 


12". 
.5         13 
.866     14 

15  . 
.866     16 
.5         17  . 

18  . 

19  . 


-     0 


.5 

.866 
L 

.866 
.5 


0 

-  .5 

-  .866  20  zz-   .866 

-1  21 zz-1 

-  .866  22  =-   .866 

-  .5  23  ..  -   zz-   .5 
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12<72^(l+5-7-ll+13+17-19—  23).  5-f  (2—  4-8-10  r  14  -1-16-20-22).  866 

+  (3-9+15-21)1 

The  circle  of  reference  of  the  third  curve  is  divided  at  every  45°  hence  : 

P* 

0"  ..................................  cos  0°=  1  12"  ..............................................  =  —  1 

1  ........  .  ......................  cos  45°=   .707  13  ......  ........................................  =—  .707 

2  ................................  cos  90°=  0  14  .....................  .'  ........................  i=  0 

3  .........  cos  135°  =  -cos  45°  =  -  .707  15  ..............................................  =   .707 

4  .............................................  =-1  16  ..............................................  =   1 

5  .............................................  =—  .707  17  ..............................................  =   .707 

6  .............................................  =  0  18  ....................................  .  .........  =  0 

7  ...........................................  '..  =   .707  '  19  ......................  :  .......................  =  -  .707 

8  .............................................  =   1  20  ..............................................  =-1 

9  .............................................  =   .707  I  21  .....  ,  ........................................  =—  .707 

10  .............................................  =  0     22  ..............................................  =  0 

11  .............................................  =-  .707  23  ..............................................  =   .707 

12p8=(0—  4+8—  12  +  16—  20)1+  (1—  3—54  7  +  9—  11—  13+15+17—  19-21  +  23).  707 


fore: 


=     0    • 
=  -   .707 
=  -1 
=  -   .707 
=     0 
=       .707 
=     1 
=       .707 
=     0 
=  -   .707 
=  -1 
=  —  .707 

12g3=(l+3-5-7  19+11-13—  15+17+19-21—  23).  707+  (2—  6+10—  14+18—  22)1 
The  divisions  of  the  circle  of  reference  of  the  fourth  curve  equal  60°  ;  we  have  there- 


0" 

sin  0°  — 

0 

12h  

1 

.      .               sin  45°  — 

.707 

13       

2 

....  sin  90°  — 

1 

14  

3 

sin  135°  —  sin  45°  — 

.707 

15 

4 

0 

16       

5 



—   .707 

17                      

6 



—  1 

18       

7 



—   .707 

19  

8 



0 

20                            ....     

9 

^_ 

.707 

21                                 

10  .... 



1 

22               ..... 

11  . 

— 

.707 

23  .. 

0" 

cos  0°—  1 

12" 

— 

1 

1 

cos  60°—   .5 

13 



5 

2  ..  . 

....  cos  120°  —  —cos  60°  —  —  .5 

14  .       



5 

3 

—  —1 

15 



1 

4  

—  —  .5 

16 



5 

5 

—   .5 

17 

4 

5 

6  

—  1 

18 



1 

7 

—   .5 

19 



5 

8  

•  —  -  .5 

20 



6 

9  

—  -1 

21 



1 

10  .... 

—  —  .5 

22 



6 

11  . 

..  -   .5 

23  .. 

— 

.5 

3+6-9+12— 15H8— 21)1  +  (1— 2— 4+5+7— 8-10+11  f  13-14— 16  + 

•17+19— 20— 22+23) -5 


0" sin  o°=  o   i  12". 

1  , sin  60°=   .866!  13 

2  sin  120°  =  sin  60°—   .866  j  14 

3  =  0  15 

4  =±-  .866  16 

5  =-  .866  17  . 

6  =  0  18  . 

7  =   .866  19 

8  =   .866  20  . 

9  =  0    i  21  . 

10  =—  .866  22  . 

11  .  ..-'==-  .866  23  . 


=  0 
=   .866 
=   .866 
=  0 
==-  .866 
=  -  .866 
=  0 
=   .866 
=   .866 
=  0 
=  -  .866 
--  .866 
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5-4— 5+7+8-10— 11+13+14— 16— 17+19+20— 22— 23).  866 
We  have  descended  to  so  many  particulars  because  we  wished  to  show  that  Father 
Dechevrens'  tables  for  the  determination  of  these  coefficients  could  be  thoroughly  relied 
on,  in  that  they  are  founded  on  Bessel's  method,  and;  moreover,  the  operations  most  ex- 
posed to  errors  are  already  verified.     The  formulae  given  by  Father  Dechevrens  are: 

Pi=(0-12)-083+ (1+23-11-13). 08+ (2+22-10-14). 072 

+  (3+21 -9-15). 059+ (4+20-8-16). 042  + (5+19-7— 171.022 
1  tfi= (6— 18).  083  +  (1-1-11-13-23).  022+  (2+10-14-22).  042 

+  (3+9-15-21). 059+ (4  +  8-16-20). 072+ (5  +  7-17-19). 08 

!=(0+ 12-6-18).  083  +  (1  +  11+13+23-5-7-17-19).  072 

+  (2+10+14+22-4-8-16-20). 042 
~—  'OI1C  n  °1).  083+  (1+5+13+17-7-11-  19-23).  042 

+  (2+4+14+16-8-10-20-22). 072 

rp8=(0+8+16-4-12-20). 083+  (1+7+9  1 15  +  17+23-3-5- -11-13-19-21). 059 
1  U3=(2+10+18-6— 14— 22).  083+  (1+3+9+11+17+19-5-7-13-15-21-23).  059 

rp4=(0-f  6+12+18-3-9-15-21).  083  +  (1  +  5+7  +  11+13 

IV  \  +17  +  19+23-2-4-8-10-14-16-20-22). 042 

I  <74=(l+2+7+8+13+14+19+20-4-5-10-ll-16- 17-22-23). 072 

Comparing  them  with  those  found  by  Bessel's  method  we  see  that  they  are  iden- 
tical ;  the  coefficients  of  the  parentheses  are  the  same  as  Bessel's  divided  by  12. 

By  means  of  these  tables  we  have  determined  the  constants  which  form  Table  IV. 

With  these  constants  we  proceed  to  the  determination  of  the  harmonic  coefficients  by 
means  of  the  formulae  which,  as  we  have  seen,  join  them.  We  give  below  the  ones  we 
made  use  of,  and  the  method  followed: 

2— tg  U. '.  log  p— log  g— log  tg  U 


u= 


log  sin  U 


-     TT  -  - 
sin  U 

In  this  example  for  January  we  should,  according  to  this,  have  : 

log  p!=  1.127429 
_  log  gi=l.  738463 
log  pi-toggi=1.388966=tog  tg  Ui  .  '.  C7i=13°  45'  37" 

log  pi  =1.127429 
log  sin  Ui=  1.376320 


log  pi-log  sin  t/i=1.751109=log  m  .'.  HI -0.5638 

One  must  remember  in  which  quadrant  U  lies.     This  rule  will  help  to  verify  it  easily. 

If  p  and  q  are  positive,  U  belongs  to  the  first  quadrant ;  if  p  is  positive  and  q  negative, 
to  the  second  quadrant;  if  both  are  negative,  to  the  third;  and  if  p  is  negative  and  q 
positive,  then  U  is  in  the  fourth  quadrant. 

Thus  to  determine  U2  in  this  month  we  find  directly  in  the  tables  the  angle  15° 
23';  but  as  p2  is  positive  and  q2  negative,  it  must  be  reduced  to  the  second  quadrant, 
which  gives 

t/2=164°  37' 

In  this  way  the  values  of  u  and  U  given  in  Table  V  have  been  deduced  for  each 
month. 

Although  this  calculation  actually  suffices  for  the  discussion  of  the  elemental  curves, 
nevertheless  we  have  calculated  separately  these  four  curves  .for  each  month. 
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The  values  of  the  ordinates  of  each  hour  together  with  those  of  the  resultant, 
appear  in  Table  VI,  to  which  we  have  added  another  column  containing  the  differences 
between  the  calculated  resultant  variation  and  the  observed  value.  The  smallness  of 
these  differences  seems  to  be  a  good  proof  of  the  correctness  of  the  calculation. 

Once  those  waves  calculated,  we  thought  that  it  would  not  be  altogether  without 
interest  to  add  the  12  plates  in  each  of  which  appear  the  elemental  harmonic  curves  of 
a  month,  together  with  the  resultant.  In  this  manner  it  will  be  easier  to  appreciate 
the  importance  of  each  of  the  curves.  From  them  it  will  be  seen  quite  clearly  why  it 
was  enough  to  calculate  the  first  four  terms  of  the  formula;  the  series  is  rapidly  con- 
vergent, so  that  the  fourth  sinusoid  is  almost  confounded  with  a  straight  line. 

Although  Father  Dechevrens  also  gives  tables  for  the  verification  of  this  last 
calculation,  we  thought  it  would  be  much  shorter,  if  not  easier,  to  work  it  out  by  the 
natural  method. 

As  an  example,  we  will  calculate  the  four  components  for  January.  The  general 
expression  for  the  ordinates  of  the  component  of  twenty-four  hours  is 


consequently  we  have  the  general  expression 

0.564  sin(13°46'+z) 

In  order  to  determine  the  angle  in  the  parentheses  for  the  different  hours,  it  is  sufficient 
to  note  that  in  the  first  wave,  the  point  which  passes  round  the  circle  of  reference,  has 
described  an  arc  of  15°  at  the  end  of  first  hour;  30°  at  the  end  of  the  second;  45°  at  the 
end  of  the  third,  etc.,  so  that  for  the  different  hours  we  have  : 

Oh__  __  sin  13  46'=  .238  7h  ____  sin  118"  46'=  sin  61°  14'=  .877 

1-  __  sin  28  46'=  .481  8  ____  sin  133°  46'=  sin  46  14'=  .722 

2  __________  __  sin  43  46'=  .692  9  _____  sin  148  46'=  sin  31  14'=  .519 

3  ..  ..sin  58  46'=  .855  10  _.   __sin  163  46'=  sin  16  14'=  .280 

4  ________  .  ________  sin  73  46'=  .960  11  _____  sinl78°46'=  sin    1   14'=  .022 

5  _„  _.sin  88  46'=!  12  _-  _____  sin  193  46'=  -sin  13°  46'=—  .238 

6  ____  sin  103  46'=sin  76"  14'=  .971 

After  the  twelfth  hour,  the  values  begin  to  repeat  in  the  same  order,  but  with 
opposite  sign.  It  then  suffices  to  multiply  each  one  of  the  natural  functions,  which, 
appear  in  the  last  member,  by  the  coefficient  0.564  in  order  to  obtain  the  values  which 
represent  the  different  ordinates  of  this  wave. 

Second  wave.  —  The  general  expression  is: 


which  is  reduced  to 

0.978  sin  (164°  37'-J-2z) 
The  variable  element  2x  increases  30°  every  hour;  hence 

Oh-__  _________________  sin  164   37'=    sin  15°  23'=    .265 

1  -  ....  sin  194°  37'=—  sin  14°  37'=-.  252 

..  sin  224°  37'=—  sin  44"  37'=  -.702 

3  ...  ._._  sin  254°  37'=—  sin  74°  37'=—  .964 

4  -  ______  sin  284°  37'=  -sin  75°  23'=—  .968 

__________  sin  314   37'=—  sin  45   23'=—  .712 

6  ________  ...  sin  344   37'=—  sin  15°  23'=  -.265 

After  the  sixth  hour  the  values  begin  to  repeat  but  with  opposite  sign,  so  that  in 
order  to  obtain  the  ordinates  of  this  curve,  all  that  is  necessary  is  to  multiply  these 
values  by  0.978,  and  go  on  repeating  this  series  of  products  till  the  twenty-third  hour, 
being  careful  to  change  signs  at  each  repetition. 
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The  other  two  waves  are  obtained  by  an  analogous  process,  taking  into  account  the 
value  of  the  variable  element  in  each  one. 


Third 

wave. 

Fourth 

wave. 

0" 

sin  11° 

1'  — 

191 

O'1  ..sin 

130° 

17'  = 

sin  49° 

43'  = 

.763 

1 

sin  56° 

1'  = 

829 

1    ..sin 

190° 

17'=  — 

sin  10° 

17'  =  — 

179 

2  ..sin 

101° 

r— 

sin  78° 

59'  = 

982 

2    ..sin 

250° 

17'=- 

sin  70° 

17'=- 

.941 

3  ..sin 

146° 

1'= 

sin  33° 

59'  = 

559 

3    ..sin 

310° 

17'=- 

sin  49° 

43'  =  - 

763 

4  ..sin 

191° 

1'  = 

sin  11° 

1'=- 

191 

This  is  the  natural  manner  of  proceeding  in  the  calculation;  in  reality  it  is  much 
simpler,  if  we  note  that  after  the  sixth  hour  in  the  first  wave,  after  the  third  hour  in 
the  second,  and  second  hour  in  the  third  wave  the  angles  become  the  complementary  of 
the  former  angles.  Hence  it  is  sufficient  to  determine  the  angles  up  to  the  fifth  hour, 
inclusive,  in  the  first  wave ;  up  to  the  second  hour  in  the  second,  and  the  first  hour  in  the 
third  wave;  and  without  paying  attention  to  the  rest,  take  at  once  from  the  trigonome- 
trical tables  the  value  of  the  sine  of  each  of  those  angles  and  their  complementaries, 
which  are  found  at  the  end  of  the  same  line.  For  example,  for  the  zero  hour  in  the  first 
wave  of  January  we  find  the  value  of  the  sine  of  13°  46',  namely  0.238,  and  at  the  end 
of  the  same  line  we  find  the  value  0.971  which  is  that  of  the  sine  of  the  complementary 
angle  76°  14',  which  we  write  in  front  of  the  sixth  hour. 

In  order  to  find  the  ordinate  for  the  zero  hour  of  January  by  means  of  Father  De- 
chevrens'  tables,  we  have  to  apply  the  following  formula : 


(0- 12) .083+  (1-13) .08+  (2-14)  .072+ (3-15) .059+  (4-16)  .042+  (5-17) . 022+  (7-19)  X 
+  (8-20)  X— .042+  (9-21)  X— .059+  (10-22)  X— . 072+  (11-23)  X— .08 


.022 


substituting  the  values  for  the  hours  in  parentheses,  and  performing  the  proper  operations. 
In  the  present  case  after  having  made  the  subtractions  in  the  parentheses,  we  have 

.31X.083+.8X.08+1.05X.072+1.1X.059+1.04X.042+.89X.022-.84X.022-.87X.042 

-79X.059— .58X. 72— .2X. 08=0.1341 

By  the  natural  method  we  can  obtain  the  same  result,  after  the  values  of  the  trigo- 
nometrical lines  have  been  found,  by  one  single  multiplication : 

0.238X0.564-0.1342 

Of  course,  Father  Dechevrens'  method  is  less  exposed  to  errors,  as  it  reduces  the 
work  to  simple  arithmetic ;  but  we  think  that  a  little  practice  in  trigonometrical  calcula- 
tion would  be  sufficient  to  save  time  by  the  natural  method. 

Epoch  of  the  critical  phases  (maxima  and  minima) . — The  determination  of  the 
epochs  of  the  critical  phases  offers  no  difficulty.  As  every  ordinate  is  expressed  by 

u  sin  (£/+#) 

the  necessary  and  sufficient  condition  in  order  that  this  value  be  a  maximum  is  that  the 
sine,  which  appears  in  the  general  expression,  be  equal  to  unity,  which  leaves  as  condi- 
tion for  the  maximum 

U+x=9Q°  or  450° 

according  as  U  is  less  or  greater  than  90°.     This  being  supposed,  let  us  remember  that 

the  variable  element  x  has  for  the  successive  hours  0,  1,  2,  3,  h,  in  the  first 

curve,   the  values   0°,   15°,   2x15°,   3x15°, fex!5°;  and  in  the  same  way  for 

the  second,  third,  and  fourth  curves,  this  variable  has  for  any  hour  the  value  7^x30°, 
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/&X45°,  fex60°,  respectively.     We  will  now  determine  the  epochs  of  the  first  maximum 
for  the  four  January  curves. 

For  the  first  we  have  as  condition: 


13° 


for  the  second: 

for  the  third : 
and  fourth: 


11° 


130° 


15  15 

450°-164°  37' 
30 

}°.'.h=lh  46m 
'.h=5h  20m 


When  we  have  determined  the  epoch  of  the  first  maximum,  we  have  determined 
the  epochs  of  all  the  other  critical  phases,  because  they  are  twelve  hours  distant  in  the 
first  curve,  six  in  the  second,  four  in  the  third,  and  three  in  the  fourth.  The  epochs 
of  the  first  maximum  for  the  four  curves  in  the  different  months  are  given  below. 

Epoch   of  the  first  maximum. 


Month. 

First 
wave. 

Second 
wave. 

Third 
wave. 

Fourth 
wave. 

January     

H.   m. 
5      5 

H.   m. 
9    31 

H.  m. 
1    46 

H.  m. 
5    20 

February 

5    47 

9    41 

1    44 

5    24 

March   ._ 

5    23 

9    45 

1     10 

5    37 

April  ... 

5      3 

9    44 

0    49 

5    38 

May 

4     10 

9    44 

7    15 

5    59 

June 

3     11 

9    49 

6    20 

0      0 

July  

2    57 

9    52 

6    27 

5    32 

August  -. 

3      9 

9    51 

6    52 

4    49 

September 

3    25 

9    41 

1    23 

4    38 

October 

4     16 

9    30 

1     11 

4    42 

November    .  .       

4    20 

9    25 

1    19 

4    42 

December     

4    26 

9    26 

1    40 

4    46 

DISCUSSION. 

In  what  has  preceeded  we  have  resolved  the  diurnal  oscillation  of  pressure  in 
Manila  into  a  series  of  simple  harmonic  oscillations,  the  first  two  of  which,  viz,  of 
twenty-four  and  twelve  hour  periods,  are  of  incomparably  more  importance  than  those 
of  eight  and  six  hour  periods,  so  that  the  principal  characteristics  of  the  diurnal  oscilla- 
tion of  pressure  are  represented  almost  exclusively  by  the  first  two. 

But  this  analysis  of  the  curves  is  in  itself  nothing  more  than  an  artifice  in  calcula- 
tion, for  every  periodic  curve  can  be  treated  by  it.  Hence,  the  component  harmonic 
curves  may  have  a  purely  mathematical  meaning  without  corresponding  to  anything  in 
reality. 

Helmholtz  in  treating  of  this  class  of  analysis  says:  "It  is  only  a  mathematical 
fiction;  admirable,  because  it  renders  calculation  easy,  but  not  necessarily  correspond- 
ing with  anything  in  reality."  Nevertheless,  speaking  in  general,  the  practical  utility 
of  this  analysis  cannot  be  denied. 

A  periodic  phenomenon  can  only  be  the  result  of  periodic  causes  of  periods  the 
submultiple  of  that  of  the  phenomenon;  otherwise  the  periodicity  would  be  disturbed. 
By  means  of  this  analysis,  the  curve  is  resolved  into  the  only  possible  elementary  ones 
with  periods  which  are  submultiples  of  the  original  one,  and  hence  it  is  most  probable 
that  the  partial  causes  which  influence  the  phenomenon  will  appear  either  totally, 
or  in  their  principal  characteristics  in  the  manner  in  which  the  elementary  curves  vary. 

127182 2 
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If  it  is  required  to  ascertain  whether  or  not  two  phenomena  are  related,  the  study  of 
the  characteristic  harmonic  curves  of  both  is  the  only  process  by  which  we  can  come  to 
this  knowledge;  the  fact  that  there  may  appear  no  trace  of  mutual  relation  in  the 
observation  curves  does  not  mean  very  much,  for  the  influence  which  the  one  may 
exercise  on  the  other  may  easily  be  disguised  by  other  causes  which  are  acting  at 
the  same  time. 

Even  though  we  may  deduce  nothing  more  from  the  comparison  of  two  phenomena 
by  means  of  this  analysis  than  the  persuasion  that  there  is  no  relationship  between 
them,  our  labor  would  not  have  been  thrown  away. 

Coming  to  the  particular  case  of  the  barometer  Dr.  Julius  Hann  1  says : 

I  was  myself  convinced  that  all  the  attempts  to  explain  the  diurnal  barometric  oscillation  by 
means  of  the  daily  variations  of  the  meteorological  elements  at  any  one  place,  as  attempted  by 
Kreil,  Blanford,  Renou,  and  others,  could  lead  to  no  conclusion;  and  I  have  published  a  series  of 
papers  giving  a  precise  description  of  the  phenomenon  as  manifested  over  the  whole  earth,  at  sea 
level  as  well  as  at  all  elevations  for  which  observations  exist,  and  I  have  endeavored  to  give  the 
results  in  such  a  form  as  would  be  suited  for  the  basis  of  a  physico-mathematical  theory.  With 
this  object  I  have  represented  all  the  results  of  observations  in  periodical  functions,  and  have 
calculated  the  amplitudes  and  phase  epochs  of  the  individual  waves,  whatever  was  their  period, 
and  which,  when  combined  together,  produce  the  complex  result  of  the  daily  barometric  curve,  such 
as  presents  itself  to  direct  observation.  *  *  *  This  is  the  only  method  which  has  enabled  us  to 
make  material  progress  in  the  comprehension  and  explanation  of  the  whole  problem. 

What  we  propose  to  do  principally  is  to  compare  the  results  obtained  for  Manila 
with  those  of  other  regions,  and  hence  we  will  put  the  following  discussion  in  such 
a  form  as  to  make  this  comparison  easy: 

First  wave. — One  of  the  most  fundamental  facts  that  have  been  established  in  the 
science  of  meteorology  is  that  the  lower  strata  of  the  atmosphere  receive  their  heat 
from  the  surface  of  the  earth  by  means  of  conduction,  and  then  transmit  it  to  the  layers 
immediately  above  by  convection.  This  movement  of  convection  is  chiefly  ascending; 
the  air  expands  laterally  in  the  higher  atmosphere  and  causes  a  fall  of  the  barometer. 
The  opposite  takes  place  with  the  cooling  of  the  surface  of  the  earth ;  the  air  in  contact 
with  it  is  cooled,  contracts,  thus  causing  descending  currents  and  a  rise  of  the  barom- 
eter; therefore  by  virtue  of  this  variation  of  temperature  there  must  be  produced 
a  maximum  and  a  minimum  in  the  atmospheric  pressure  about  the  hours  of  minimum 
and  maximum  temperature. 

The  exchange  of  heat  between  the  surface  of  the  earth  and  the  lower  layers  of 
the  atmosphere  is  of  a  periodic  character,  the  principal  period  being  one  of  twenty- 
four  hours.  We  say  principal  period,  because  this  variation  is  not  simply  harmonic; 
so  that  if  it  were  to  be  analyzed  separately  it  would  not  be  represented  by  a  simple 
harmonic  motion,  but,  in  addition  to  a  fundamental  component  of  a  twenty-four  hour 

period,  there  would  result  others  less  important  of  twelve,  eight,  six,  hour 

periods. 

The  first  wave  in  the  analysis  of  the  barometric  curve  is  one  of  a  twenty-four  hour 
period,  and  invariably  presents  the  maximum  and  minimum  about  the  hours  of  minimum 
and  maximum  temperature.  This  general  character  of  the  curve  makes  it  likely  that  it 
represents  the  principal  character  of  the  variation  of  pressure,  caused  by  the  exchange 
of  heat  between  the  surface  of  the  earth  and  the  lower  layers  of  the  atmosphere. 

Both  the  amplitude  and  the  epochs  of  the  critical  phases  of  this  wave  present  special 
characteristics  in  Manila,  which  it  is  well  to  note.  The  amplitude  is  less  and  the  epoch 
of  the  maximum  occurs  earlier  than  in  the  greater  number  of  the  stations  for  which 

'Further  contributions  to  the  foundation  of  a  theory  of  the  daily  barometric  oscillation.  Q. 
Journal  of  the  Roy.  Met.  Soc.,  London,  1899,  p.  40. 
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this  analysis  has  been  made.  Hann  made  note  of  these  peculiarities  for  Mauritius, 
Singapore,  and  Batavia  in  his  "Untersuchungen  iiber  den  taglichen  Gang  des  Barometers." 
Mr.  Eliot  includes  in  the  following  table  these  three  stations,  together  with  several 
others,  studied  by  him  in  India,  in  order  to  compare  them  with  Trivadrum,  the  only 
station  in  India  which  presents  these  characteristics.  To  these  we  add  Manila  for 
comparison. 


Station. 

North 
latitude. 

Ampli- 
tude Ul 
for  mean 
of  year. 

Epoch  of 
max- 
imum. 

Phase  Ui 
for  mean 
of  year. 

Station. 

North 
latitude. 

Ampli- 
tude Ml 
for  mean 
of  year. 

Epoch  of 
max- 
imum. 

Phase  Ui 
for  mean 
of  year. 

o     / 
8    31 

0.413 

H.    m. 
4    36 

O       r 

21      2 

Allahabad 

o     > 
25    26 

0  774 

H.    m. 
1    38 

o     / 
335    33 

Madras 

13      4 

.595 

6      1 

359    38 

Magpur 

21      9 

767 

7      4 

343    52 

Bombay  --  

18    54 

.472 

7    56 

330    53 

Lahore 

31    34 

.667 

8    32 

322      2 

Rangoon     

16    46 

..685 

6    20 

355      3 

Mauritius  

20      6S 

.310 

2    54 

46    35 

Chittagong       .  . 

22    21 

.582 

6    50 

347    45 

Singapore  

1    17 

.531 

4    18 

25    35 

Calcutta  (Alipore) 

22    32 

.667 

6    56 

346      5 

Batavia 

6    US 

626 

4    19 

25    17 

Kurrachee 

24    47 

.497 

7      5 

343    49 

Manila 

14    35 

.544 

4     15 

25    58 

12    45 

.779 

7    26 

338    24 

• 

Mr.  Eliot  terminates  his  discussion  of  this  first  wave  with  the  following  words: 

The  first  component  of  the  diurnal  oscillation  at  Trivandrum,  hence,  exhibits  characteristic  features 
differing  largely  from  those  of  other  stations  in  India,  but  resembling  strongly  those  at  Batavia  and 
Singapore  and  (to  a  less  extent)  Mauritius.  The  chief  features  are  its  small  amplitude  and  the  early 
epoch  of  the  maximum  phase.  These  are  probably  due  mainly,  if  not  entirely,  to  .its  maritime  position 
and  low  latitude. 

The  similarity  which  Mr.  Eliot  makes  mention  of  as  existing  between  Trivandrum, 
Batavia,  and  Singapore  appears  much  more  marked  in  the  case  of  Manila,  both  in  the 
amplitude  and  initial  phase  as  well  as  the  epoch  of  the  maximum,  and  consequently 
Manila  is  an  example  which  confirms  the  conjecture  concerning  these  properties  in  coast 
stations  of  low  latitudes. 

Annual  variation  of  u^ — The  annual  variation  of  u^  and  the  annual  variation  of  the 
absorption  of  heat  in  the  lower  layers  of  the  atmosphere  have  between  them  the  relation 
of  cause  and  effect,  as  was  pointed  out  by  Hann  ("Untersuchungen,  etc.")  and  Angot 
("Annales  du  Bureau  central  Meteorologique  de  France,"  1887,  p.  B.  297),  and  confirmed 
by  the  work  of  Mr.  Eliot  in  Indian  Meteorological  Memoirs,  Vol.  XII,  p.  374. 

A  little  consideration  will  show  that  the  same  takes  place  in  Manila.  In  the  "Climato- 
logia"  published  by  Rev.  Jose  Coronas  in  the  "El  Archipielago  Filipino"  it  is  shown  that  the 
months  of  February,  March,  and  April  are  the  driest  and  least  cloudy  months  of  the  year, 
while  July,  August,  and  September  are  the  rainiest  and  cloudiest.  These  circumstances 
are  extremely  favorable  during  the  first  three  months  for  the  absorption  of  heat  in  the 
lower  layers  of  the  atmosphere,  and  very  unfavorable  during  the  last  three  months,  in 
that  the  clouds  intercept  a  large  amount  of  the  calorific  energy,  and  of  that  which  passes 
to  the  surface  of  the  earth,  some  is  converted  into  latent  heat  of  evaporation. 

A  glance  at  the  values  of  u^  makes  it  clear  that  the  three  months  during  which  it 
reaches  its  maximum  are  February,  March,  and  April,  in  which  the  absorption  of  the 
calorific  energy  in  the  lower  layers  of  the  atmosphere  must  be  a  maximum ;  and  it  reaches 
its  minimum  values  in  July,  August,  and  September  when  the  absorption  is  least. 

Mr.  Eliot  mentions  in  the  publication  already  cited  that  during  the  transition  period 
between  the  rainy  and  dry  seasons  a  relative  maximum  of  the  value  of  u^  usually  takes 
place.  This  relative  maximum  appears  in  Manila  in  October  which  may  be  considered  as 
the  border  month  between  the  two  seasons. 

The  second  wave. — There  are  two  characteristics  which  are  noteworthy  in  this  wave ; 
first,  that  it  is  much  more  important  than  the  others  in  all  the  months  and  in  all  the 
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regions  where  it  has  been  studied,  and,  secondly,  the  constancy  of  the  epoch  of  the  critical 
phases. 

In  the  following  table  are  given  the  epochs  of  the  first  maximum  found  by  Mr.  Eliot 
for  several  stations  of  India: 


North 
latitude. 

Station. 

Janu- 
ary. 

Feb- 
ruary. 

March. 

April. 

May. 

June. 

July. 

Au- 
gust. 

Sep- 
tember. 

Octo- 
ber. 

Novem- 
ber. 

Decem- 
ber. 

O              ' 

8      31 
10      50 
13      04 
18      54 
24      47 
34      10 
Mean  vs 
14      35 

Trivandrum 

H.  m. 
9  30 
9  50 
9  45 
9  43 
9  41 
9  48 
9  42 
9  31 

H.    m. 

9    38 
9    59 
9    52 
9    51 
9    58 
10      3 
9    53 
9    41 

H.    m. 

9    34 
9    59 
9    50 
9    53 
9    58 
10      6 
9    50 
9    45 

H.    m. 

9    36 
9    58 
9    49 
9    54 
9    57 
9    58 
9    47 
9    44 

H.  m. 
9  37 
9  52 
9  48 
9  54 
9  58 
9  52 
9  47 
9  44 

H.    m. 
9    42 
9    52 
0    53 
10      2 
10    10 
9    50 
9    54 
9    49 

H.    m. 
9    47 
9    53 
9    57 
10      6 
10      9 
9    58 
9    58 
9    62 

H.    m. 

9    41 
9    53 
9    53 
9    57 
9    51 
10  .22 
9    55 
9    51 

H.    m. 
9    29 
9    48 
9    39 
9    45 
9    49 
9    47 
9    41 
9    41 

H.    m. 

9  16 
9  38 
9  28 
9  25 
9  28 
9  49 
9  29 
9  30 

H.    m. 
9    14 
9    39 
9    27 
9    23 
9    29 
9    49 
9    28 
9    25 

H.    m. 
9    22 
9    47 
9    33 
9    36 
9    35 
9    49 
9    36 
9    26 

Trichinopoly 

Madras  ._  _    

Bombay  

Leh                        .  .    . 

lue  for  tropical  India 

Manila 

This  wave  with  its  double  daily  oscillation  is  the  one  which  produces  the  char- 
acteristic tides  of  the  barometer.  It  is  therefore  of  the  greatest  possible  interest  to 
discover  its  cause. 

The  discussions  of  Hann,  Angot,  and  Lord  Kelvin  on  this  point  are  too  voluminous 
to  reproduce  here  in  all  their  details,  so  we  will  try  and  summarize  them  as  shortly 
and  clearly  as  possible. 

The  phenomenon,  represented  by  this  wave  must  have  the  two  properties,  viz,  that 
it  have  greater  influence  on  the  barometer  than  any  other,  and,  secondly,  that  it  be 
independent  of  local  influences. 

From  the  outset,  it  is  evident  that  the  cause  of  this  phenomenon  can  have  no  other 
origin  than  the  influence  of  the  sun  on  the  terrestial  atmosphere.  "But,"  says  Lord 
Kelvin,  "it  cannot  be  owing  to  the  gravitational  tide-generating  influence  of  the  sun, 
since  we  would  have  a  much  greater  influence  of  the  same  kind  due  to  the  moon, 
whereas  in  reality  the  lunar  barometric  tide  is  insensible  or  nearly  so."  It  appears 
therefore  certain  that  it  must  be  attributed  to  the  variation  of  the  temperature.  But 
this  variation  in  the  lower  layers  of  the  atmosphere  is  subject  to  all  possible  variations 
caused  by  local  circumstances  of  the  different  stations,  whereas  the  periodic  action  of 
solar  rays  on  the  higher  layers,  recurring  day  by  day,  must  necessarily  produce  periodic 
movements  of  great  regularity. 

This  easily  explains  the  regularity  of  this  wave.  What  does  not  appear  so  clearly 
is  the  reason  of  its  diurnal  double  oscillation  and  its  great  amplitude.  The  oscillation 
of  the  temperature  in  the  higher  layers  of  the  atmosphere  must  be  analogous  in  its 
general  form  to  the  oscillation  of  the  temperature  in  the  lower  layers,  since  the  causes 
which  produce  them  are  the  same  —  the  insolation  during  the  day  and  the  radiation 
during  the  night.  The  first  causes  the  maximum  to  take  place  shortly  after  the  sun 
has  reached  its  maximum  altitude,  and  the  second  produces  the  minimum  shortly  before 
sunrise.  How,  therefore,  can  this  thermic  oscillation  of  one  maximum  and  one  minimum 
give  rise  to  the  double  barometric  oscillation  which  we  are  studying? 

To  solve  this  problem,  Hann  has  recourse  to  the  only  means  which  exists  for 
ascertaining  the  relation  between  two  phenomena.  He  resolves  the  diurnal  variation 
of  the  temperature  into  its  harmonic  components,  and  finds  that  of  the  two  principal 
waves,  the  first  (of  twenty-four  hour  period)  appears  constantly  with  an  amplitude 
greater  than  the  second  (of  a  twelve  hour  period)  opposite  to  what  takes  place  with 
the  pressure.  For  exampl 


Latitude  6°  0',  equatorial  Pacific  .......................................................  0.°87  O.°ll 

Latitude  16°  0',  tropical  Pacific  .........................................................  l.°39  O.°31 

Latitude  33°  0',  subtropical  Pacific  ...................................................  0.°97  0.°25 

Latitude  48°  12',  Temperate  Zone  (Vienna)  ...................................  2.°86  O.°51 
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In  spite  of  this  disagreement  in  the  harmonic  components  of  the  two  phenomena, 
Hann  did  not  hesitate  to  assign  as  cause  of  the  double  oscillation  of  the  barometer 
this  double  oscillation  of  the  temperature. 

When  this  assertion  of  the  eminent  meteorologist  was  published,  two  difficulties 
were  raised  by  scientists:  First,  that  the  thermic  variation  analyzed  was  that  of  the 
lower  layers  of  the  atmosphere,  whereas  the  point  at  issue  was  the  variation  in  the 
higher  layers  and  of  this  it  was  questionable  whether  the  oscillation  was  of  exactly  the 
same  character  as  in  the  lower.  The  second  difficulty  was  that  this  double  daily  tem- 
perature wave  was  only  a  result  of  calculation,  a  mathematical  fiction,  which  could  not 
serve  in  the  explanation  of  a  real  phenomenon.  There  is  in  reality,  no  daily  variation 
of  temperature  with  two  maxima  and  two  minima;  and  if  we,  in  spite  of  this,  obtain 
a  double  daily  temperature  wave,  because  we  insist  on  representing  the  daily  march  of 
temperature  by  a  series  of  sines,  this  forced  mathematical  form  can  never  serve  to 
explain  an  observed  phenomenon,  whereas  for  this  some  real  natural  process  must  be 
sought  for  as  a  cause. 

The  first  difficulty  is  solved  by  what  was  said  before  concerning  the  variation  of  the 
temperature  in  the  higher  layers  of  the  atmosphere;  the  causes  which  produce  it  are 
the  same ;  therefore  it  must  be  of  the  same  general  form  as  that  observed  in  the  lower. 

The  solution  of  the  second  difficulty  requires  a  little  more  explanation.  It  is  known 
that  if  a  limited  mass  of  liquid  or  gaseous  fluid  be  set  in  simple  pendulumlike  oscilla- 
tions, the  amplitude  of  these  oscillations  are  governed  by  the  conditions  proper  to  the 
fluid  (its  dimensions,  temperature,  etc.),  so  that  if  the  mass  of  fluid  receives  a  single 
impulse,  it  will  always  move  with  that  class  of  undulations  which  is  possible  to  it  by 
reason  of  its  peculiar  conditions.  These  oscillations  are  commonly  called  free  oscillations. 

If  the  impulse  recurs  periodically,  then  the  mass  is  forced  to  take  up  oscillations  of 
that  period,  even  though  these  new  oscillations  do  not  coincide  with  any  of  the  forms 
of  oscillation  which  belong  to  the  free  waves.  If  this  periodic  impulse  is  not  of  a  simple 
harmonic  character — that  is  to  say,  if  it  cannot  be  represented  by  a  simple  sine  curve — 
then  it  can  be  resolved,  according  to  Fourier,  into  its  harmonic  components  of  periods 
1,  i,  i,  etc.,  so  that  the  effect  due  to  the  impulse  would  be  the  same  as  that  produced  by 
several  simple  harmonic  impulses  acting  together.  It  is  quite  immaterial  to  the  question 
whether  or  not  these  component  impulses  exist  in  practice;  it  is  sufficient  if  the  real 
impulse  can  be  considered  as  the  resultant  of  them. 

If  the  period  of  one  of  these  component  impulses  coincides  with  that  of  the  free 
waves  of  the  mass  of  the  fluid,  then  the  resulting  forced  oscillation  will  attain  a  dis- 
proportionately great  amplitude;  so  that  although  the  amplitude  of  this  impulse  be  rel- 
atively small  compared  with  that  of  the  other  component  impulses,  the  effect  produced 
by  it  is  much  larger;  whence  it  appears  that  if  the  resulting  complex  oscillation  be 
analyzed  by  means  of  a  sine  series,  the  harmonic  curve  corresponding  to  the  free  oscil- 
lations will  be  of  large  amplitude,  while  the  corresponding  wave  in  the  analysis  of  the 
impulse  will  be  of  small  amplitude,  and  yet,  that  is  no  reason  why  the  second  cannot  be 
the  cause  of  the  first. 

To  come  now  to  the  case  before  us :  The  curve  which  represents  the  diurnal  variation 
of  the  temperature  in  the  higher  strata  of  the  atmosphere  must  necessarily  be  asymmet- 
rical, because  its  two  causes,  insolation  and  radiation,  which  depend  on  the  duration  of 
day  and  night,  are  asymmetrical.  If,  therefore,  on  analyzing  it,  we  find  that  it  is  repro- 
duced almost  exclusively  by  the  combination  of  the  two  harmonic  components  of  periods 
cf  twenty-four  and  twelve  hours,  then  we  may  consider  the  real  variation  which  we  ac- 
tually observe  as  the  resultant  of  two  harmonic  variations  of  the  temperature;  one  of  a 
twenty-four  hour  period  and  the  other  of  twelve;  and  it  does  not  really  matter  whether 
or  not  there  actually  exist  two  causes,  which  give  rise  to  these  two  variations;  it  is 
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enough  for  our  proposition  that  the  phenomenon  be  reproduced  in  the  same  manner  as 
if  they  really  existed. 

Doctor  Margules,  following  the  indications  made  by  Lord  Kelvin,  undertook  the 
difficult  and  laborious  task  of  calculating  the  oscillations  caused  in  the  atmosphere  by 
its  periodic  heating  and  cooling.1  In  the  first  of  his  papers  he  showed  that  the  period 
of  free  oscillation  of  our  atmosphere  is  nearly  twelve  hours. 

These  studies  and  conclusions  of  Doctor  Margules  have  been  received  with  a  good  deal 
jo  tioiseooo  aqq.  uaeq  9A«q  pus  uisBisniftue  joseveral  interesting  publications  such  as  Paul 
Jaerisch's  "Zur  Theorie  der  hydrodynamischen  Gleichungen  in  spharischen  Koordinaten" 
(Meteorologische  Zeitschrift,  1907,  pp.  481-498)  ;  Trabert's  "Die  Theorie  der  taglichen 
Luftdruckschwankung  von  Margules,  und  die  tagliche  Oszillation  der  Luftmassen"  (Met. 
Zeit.,  1903,  pp.  481-501)  ;  W.  J.  Humphreys'  "On  the  diurnal  variations  of  atmospheric 
pressure"  (Bulletin  of  the  Mount  Weather  Observatory,  1912,  vol.  5,  part  2,  pp.  132-156) 
and  the  most  complete  of  all,  Horace  Lamb's  "Atmospheric  Oscillations"  (Proc.  Royal 
Society,  1911,  vol.  84,  pp.  511-572)  .2 

We  can,  therefore,  explain  the  double  oscillation  and  the  great  amplitude  of  the  second 
wave  thus :  When  the  atmosphere  receives  a  thermic  impulse  (if  we  may  make  use  of  this 
expression)  there  is  produced  in  it  the  same  effect  as  would  be  produced  by  two  simul- 
taneous impulses  which  could  be  reproduced  by  sine  curves  of  twenty-four  and  twelve 
hour  periods,  respectively.  The  amplitude  or  intensity  of  the  second  is  less  than  that  of 
the  first,  but  the  oscillations  which  it  causes  in  the  atmosphere  must  be  of  great  amplitude 
because  its  period  coincides  with  that  of  the  free  oscillations ;  while  the  oscillations  pro- 
duced by  the  first  are  considerably  reduced  because  their  period  does  not  coincide  with 
that  of  the  free  oscillations. 

The  hypothesis  that  the  second  wave  is  due  to  the  absorption  and  radiation  of  heat  in 
the  higher  strata  of  the  atmosphere  is  confirmed  by  the  comparison  of  the  annual  variation 

1  "Uber   die   Schwingungen   periodisch   erwarmter  Luft."     (Sitzungsberichte,   Wien.     Akad.   Marz, 
1890.) 

"Luftbewegungen  in  iener  rotirenden  Spharoidschale."     (Ibid.  April,  1890;  Jan.  and  Dec.,  1893.) 

2  This  work  was  already  in  the  press  when  the   Observatory  was  in  receipt  of  Prof.   Frank  H. 
Bigelow's  "The  Thermodynamics  of  the  Circulation  and  Radiation  of  the  Earth's  Atmosphere"  published 
as  Bulletin   No.  3  of  the  Argentine   Meteorological  Office.     On  page  15,  he  proposes  two  difficulties 
against  Lord  Kelvin's  theory;  first,  that  the  semidiurnal  pressure  waves  disappear  within  2,000  meters 
of  the  surface;  and,  secondly,  that  the  heat  term  (Qi— Qo)   is  not  included  in  the  general  equation  of 
motion  as  employed  by  Margules,  Jaerisch,  and  Gold.     It  is  not  the  first  time  that  Professor  Bigelow 
has  impugned  this  theory;  see,  for  example,  his  "Studies  on  the  diurnal  periods  in  the  lower  strata  of 
the  Atmosphere"  (Monthly  Weather  Review,  1905,  p.  93)   and  chapter  9,  p.  458,  of  his  "Report  on  the 
international  cloud  observations."     At  the  same  time  it  is  not  to  be  supposed  that  scientists  like  Lord 
Kelvin  and  Hann  were  ignorant  that  the  semidiurnal  oscillation   diminishes   with   altitude  till   it  dis- 
appears or  becomes  almost  imperceptible,  and  yet,  notwithstanding  this,  these  scientists  did  not  hesitate 
to  establish  their  theory.     Hence,  it  must  be  supposed  the  semidiurnal  oscillation  reappears  in  the  higher 
atmosphere.     Mr.  Ernest  Gold  seems  to  indicate  this  in  his  interesting  memoir  "The  International  Kite 
and  Balloon  ascents"  which  won  the  first  prize  offered  by  the  German  Meteorological  Society  in  1908. 
For,  on  page  120,  he  says:  "Hence  changes  of  pressure  at  the  earth's  surface  are  connected  with  changes 
in  the  opposite  direction  in  the  stratosphere,  which  give  rise  to  currents  flowing  between  the  stratosphere 
and  troposphere." 

The  second  difficulty  would  seem  to  indicate  that  the  three  eminent  mathematicians  had  over- 
looked the  influence  which  the  free  heat  term  could  have  on  the  general  equation  of  motion,  whereas 
Doctor  Margules,  after  introducing  into  the  differential  equation  of  motion  the  free  heat  term,  proceeds 
thus:  "For  our  purpose  it  will  be  more  convenient  to  consider  the  pressure  variations  as  a  consequence 
of  the  temperature  variations,  not  as  a  consequence  of  the  variable  flow  of  heat.  We,  therefore,  return 
to  equation  (4a),"  where  the  term  of  which  Professor  Bigelow  speaks  certainly  does  not  enter. 
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of  these  two  factors  and  that  of  the  amplitude  «2.  A  brief  explanation  of  this  will  suffice 
for  our  purpose.1 

One  of  the  most  important  factors  in  the  influence  which  the  absorption  and  radiation 
may  have  on  changes  of  pressure,  is  the  presence  or  absence  of  clouds  in  the  atmosphere. 
The  amount  of  heat  absorbed  by  a  dry  and  clear  atmosphere  is  employed  in  expanding  it, 
thus  causing  lateral  and  ascending  currents,  which  are  shown  on  the  barometer  by  changes 
of  pressure.  On  the  other  hand,  the  calorofic  energy  which  reaches  an  atmosphere 
charged  with  clouds  is  converted,  in  great  part,  into  latent  heat  of  evaporation,  which 
can  not  have  any  effect  on  the  pressure.  Hence,  in  the  dry  season  the  changes  of  pressure 
due  to  these  causes  will  be  greater  than  during  the  wet  season,  on  account  of  the  great 
difference  of  cloudiness  between  these  two  seasons. 

This  difference  in  the  values  of  u2,  which  are  constantly  greater  during  the  months 
constituting  the  dry  season  than  during  the  months  of  wet  season,  is  seen  quite  clearly 
in  Table  V,  as  is  also  the  difference  in  values  in  the  perfectly  defined  separation  of  the 
two  seasons  from  May  to  June. 

We  will  terminate  the  discussion  of  this  wave  by  comparing  its  principal  characters 
in  Manila  with  those  of  the  stations  mentioned  in  connection  with  the  first  wave. 

The  following  table  is  copied  from  Volume  X  of  Indian  Meteorological  Memoirs,  page 
74,  with  the  exception  that  we  have  converted  the  original  observations  in  inches  to  milli- 
meters in  order  to  facilitate  the  comparison : 
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From  this  table  one  can  appreciate  in  a  general  way  the  constancy  of  the  principal 
characteristics  of  this  wave  in  the  different  stations;  the  almost  perfect  agreement  in 

1  See  a  more  detailed  discussion  of  this  point  in  Vol.  XII  of  "Indian  Meteorological  Memoirs," 
pages  328-347. 

The  following  may  also  be  consulted:  "Researches  on  solar  heat  and  its  absorption  by  the  earth's 
atmosphere"  by  Langley  in  No.  15  of  "Professional  Papers  of  the  Signal  Service,"  and  in  Vol.  I  of 
"Annals  of  Astrophysical  Observatory  of  the  Smithsonian  Institution,"  and  "Preparatory  Studies  for 
Deductive  Methods  in  Storm  and  Weather  Predictions"  by  Cleveland  Abbe. 
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the  ratio  between  the  maximum  and  minimum  amplitude,  and  the  epoch  of  the  year 
in  which  the  minimum  takes  place,  is  especially  evident. 

A  few  discrepancies  will,  however,  be  noted,  which  show  that  this  wave  is  not  quite 
free  from  the  influence  of  local  circumstances.  Thus,  for  example,  in  some  regions  the 
amplitude  has  two  maxima  and  two  minima  in  the  year,  while  in  others,  Manila  among 
them,  this  double  variation  disappears.  In  stations  like  Batavia,  Mauritius,  and  Aden 
which  have  the  absolute  maximum  in  October  or  November,  there  is  another  relative 
maximum  in  March  or  April. 

The  reason  of  these  variations  is  fairly  obvious.  The  second  wave,  as  it  is  given 
by  analysis,  must  be  a  complex  wave  resulting  from  the  interference  of  two  with 
the  same  period  of  twelve  hours.  One  of  these  (which  we  will  call  principal  semi- 
diurnal wave)  is  the  product  of  the  variations  of  temperature  in  the  higher  strata  of 
the  atmosphere,  as  was  explained  above;  the  other  (secondary  semidiurnal  wave)  rep- 
resents the  oscillations  which  the  semidiurnal  component  of  the  temperature  in  the 
lower  strata  must  produce  in  the  .pressure.  The  first  of  these  two  presents,  according 
to  Hann  and  Angot,  a  double  variation  during  the  year,  with  maxima  at  the  equinoxes 
and  minima  at  the  solstices.  The  second,  which  is  more  irregular  and  of  less  amplitude 
than  the  first,  owing  to  local  influences,  has  one  maximum  only  (during  the  dry  season) 
and  one  minimum  (in  the  wet  season)  . 

It  is  easy  to  see  that  the  interference  of  these  two  waves  can  cause  the  complete  or 
almost  complete  disappearance  of  the  double  annual  variation  in  the  wave  we  are  dis- 
cussing; for  the  autumn  equinox  occurs  usually  in  the  wet  season  in  the  great  majority 
of  the  stations,  and  the  winter  solstice  generally  occurs  about  the  beginning  of  the  dry 
season,  so  that  the  maximum  and  minimum  corresponding  to  these  epochs  in  the  principal 
wave  are  held  back  either  entirely  or  in  great  part,  because  they  coincide  with  the 
epochs  of  the  minimum  and  maximum  amplitude  of  the  secondary  wave.  The  opposite 
happens  in  the  summer  solstice  and  the  spring  equinox,  because  they  coincide  with 
the  wet  and  dry  seasons  respectively. 

M.  Angot,  in  the  publication  already  cited  (p.  B311  sq)  proposes  an  ingenious 
method  for  the  separation  of  these  two  waves.  It  would  take  up  too  much  space  to  explain 
the  whole  proceeding,  so  we  simply  give  in  the  following  table  the  results  obtained  by  his 
method  for  Manila.1 

1  The  formula  we  finally  obtain  by  this  method  of  Angot  is 

O.mm926  ?  cos4X  cos  (2m+64°—  2e) 


where  h  represents  the  mean  pressure;  8  the  declination  of  the  sun;  r  the  radius  of  the  earth's 
orbit;  X  the  latitude  of  the  place;  m  the  different  hours  expressed  in  degrees;  e  the  equation  of 
time.  The  values  which  we  have  taken  for  these  quantities  are: 

fc=759.10  X=14°35'  cog25 

2e  ~r2~^ 

January  ....................................  +4°  48'  ........................................................  0.9033 

February  ..................................  +7°  01'  ........................................................   0.9722 

March    ................  :  .....................  +4°  19'  ........................................................  1.0043 

April    ........................................  +0°  02'  ........................................................   0.9597 

May  .............  :.  ...........................  -1°  46'  ........................................................  0.8734 

June  ..........................................  +0°  08'  ........................................................  0.8197 

July   ..........................................  +2°  41'  ........................................................  0.8424 

August  ......................................  +1°  50'  ........................................................  0.9208 

September  ................................  -2°  32'  ........................................................  0.9852 

October   ....................................  -7°  00'  .................................................  0.9806 

November  ........  .  .......................  -7°  08'  ........................................................  0.9193 

December  ...                          ......  -1°  56'  ........................................................  0.8746 
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Hour. 

January. 

Feb- 
ruary. 

March. 

April. 

May. 

June. 

July. 

August. 

Septem- 
ber. 

October. 

• 
Novem- 
ber. 

Decem- 
ber. 

0            

-f38 

+43 

+41 

+34 

+29 

+29 

+33 

+35 

+32 

+26 

+24 

+29 

1        ..   

+  1 

f  4 

0 

—  5 

rj 

-  4 

_.        O 

—  3 

—  9 

—15 

—15 

—  7 

2 

—36 

—36 

—40 

—44 

—41 

-37 

-36 

—40 

-48 

-52 

-49 

-42 

3              

—63 

-66 

-70 

-70 

-65 

—60 

—60 

C* 

Ow 

-73 

—75 

—71 

-65 

4          

—73 

—79 

-82 

-78 

-71 

—  6« 

—68 

—75 

—80 

-78 

—73 

—71 

5 

—64 

—70 

—71 

—65 

—58 

-55 

—58 

—63 

—64 

-«o 

—55 

—58 

6           

—38 

—43 

—41 

—34 

—29 

—29 

—33 

—35 

—32 

—26 

—24 

—29 

7        

—  1 

—  4 

0 

+  5 

+  7 

+  4 

+  2 

+  3 

+  9 

+  15 

+  15 

+  7 

8                   -  - 

+36 

+36 

+40 

-t-44 

+41 

+37 

+36 

+40 

+48 

+  52 

+49 

+42 

9          

+63 

+66 

+  70 

+70 

+  65 

+60 

+60 

+66 

+73 

+75 

+71 

+65 

10 

+73 

+79 

+  82 

+78 

+71 

-t66 

+68 

+75 

+  80 

+78 

+73 

+  71 

11               

+64 

+  70 

+  71 

+65 

+58 

+  55 

+58 

+63 

+  64 

+60 

+  55 

+58 

12            

+38 

+43 

+41 

+34 

+29 

+29 

+  33 

+35 

+32 

+26 

+24 

+29 

13 

+  1 

+  4 

0 

-  5 

—  7 

—  4 

—  2 

—  3 

—  9 

—15 

—15 

-  7 

14 

—36 

—36 

-40 

—44 

—41 

—37 

—36 

—40 

—48 

-52 

—49 

—42 

15        

—63 

—66 

-70 

—70 

-65 

—60 

—60 

—66 

—73 

—75 

—71 

-65 

16 

—73 

—79 

—82 

—78 

—71 

-66 

—68 

—75 

-80 

—78 

—73 

—71 

17 

—64 

—70 

—71 

—65 

—58 

—55 

—58 

-63 

—64 

—60 

—55 

—58 

18            

—38 

—43 

—41 

—34 

—29 

—29 

—33 

—35 

-32 

—26 

—24 

-29 

19 

—  1 

—  4 

0 

+  5 

+  7 

+  4 

+  2 

+  3 

+  9 

+  15 

+15 

+  7 

20 

+36 

+36 

+40 

+44 

+41 

+37 

+36 

+40 

+48 

+52 

+49 

+42 

21          

+63 

+66 

+70 

+70 

+65 

+60 

+60 

+66 

+73 

+75 

+71 

+65 

22 

-  73 

+79 

+82 

+78 

+71 

+  66 

-J-68 

+75 

+  80 

+  78 

+73 

+  71 

23 

+  64 

-70 

+71 

+65 

+58 

+65 

+58 

+  63 

+  64 

+60 

+55 

+58 

NOTE.— Black-face  type  indicates  maximum  or  minimum. 

We  cannot  compare  our  results  with  those  of  other  regions  because  we  do  not  know 
whether  this  analysis  has  been  worked  out  for  others  than  the  nine  calculated  by 
Angot,  and  for  these  Angot  took  only  one  or  two  months;  nevertheless,  we  think  it 
useful  to  have  proved  for  Manila  this  character  attributed  to  the  principal  semidiurnal 
wave. 

Third  wave. — There  are  only  two  possible  hypotheses  with  regard  to  the  origin  of 
this  wave,  viz: 

1.  That  it  is  produced  by  an  independent  action  of  an  eight- hour  period. 

2.  That  it  is  part  of  the  effect  due  to  the  causes  producing  the  first  wave ;  since,  as 
was  indicated  in  the  discussion  of  this  wave,  its  causes  cannot  be  represented  by  a  simple 
harmonic  motion. 

The  first  objection  against  the  first  hypothesis  is  that  no  independent  physical  action 
of  an  eight-hour  period  is  known.  And  it  may  be  added  that  the  amplitude  of  the  third 
component  is  very  small,  resembling  in  character  that  of  a  residual  rather  than  the  total 
effect  due  to  an  independent  cause. 

In  confirmation  of  these  general  reasons  is  the  property  which  is  noted  in  this  wave 
in  all  regions,  namely,  that  the  annual  variation  of  its  amplitude  and  the  epoch  of  its 
critical  phases  follow  a  law  like  to  the  one  for  the  first  wave,  with  this  particular,  that 
the  sign  of  this  variation  is  almost  absolutely  opposite,  which  raises  the  suspicion  that 
the  third  wave  is  a  correction  of  the  values  of  the  first. 

Given  the  slight  influence  which  this  as  well  as  the  fourth  wave  can  have  on  the 
barometric  variation,  we  believe  that  it  is  scarcely  worth  while  discussing  the  matter. 
We  do  not  wish,  however,  to  omit  a  peculiarity  of  the  third  wave,  noted  by  Doctor  Hann, 
which  is,  that  toward  the  equinoxes  the  epoch  of  the  phases  undergoes  a  notable  inversion, 
and  that  this  inversion  coincides  exactly  with  the  epoch  of  the  minimum  amplitudes. 
This  property  appears  well  marked  in  Manila ;  the  inversion,  spoken  of  by  Doctor  Hann, 
takes  place  from  April  (Oh  49m)  to  May  (7h  15m)  and  from  August  (6h  52m)  to  Septem- 
ber (lh  23) m;  the  minima  of  the  amplitude  also  take  place  in  April-May  and  August- 
September,  as  may  be  seen  in  Table  V. 

CONCLUSION. 

From  what  has  been  said  up  to  the  present,  it  is  deduced  that  the  modern  theory 
of  the  double  daily  oscillation  of  the  barometer  supposes  two  causes  of  the  phenomenon. 
The  one  is  quite  free  from  local  influences,  and  is  precisely  the  one  which  introduces  into 
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the  barometric  curve  that  double  oscillation  which  has  always  been  so  mysterious  to 
those  who  have  examined  its  cause.  It  is  represented  by  what  we  have  called  the  prin- 
cipal semidiurnal  wave,  and  is  due  to  the  variation  of  the  temperature  in  the  higher  layers 
of  the  atmosphere.  The  other  varies  according  to  the  local  circumstances  of  each  region, 
in  that  it  is  due  to  the  variation  of  temperature  in  the  lower  layers.  The  curve  represent- 
ing it  is  usually  called  the  thermic  wave.  In  the  following  table  we  give  its  values  for 
all  the  months  of  the  year;  they  are  obtained  by  subtracting  the  principal  semidiurnal 
wave  from  the  observed  variation  (Table  III)  : 

Diurnal  or  thermic  "wave"   (unit  0.01). 


Hour. 

January. 

Febru- 
ary. 

March. 

April. 

May. 

June. 

July. 

August. 

Septem- 
ber. 

October. 

Novem- 
ber. 

Decem- 
ber. 

0 

+    8 

+  .6 

+  18 

+     27 

+  33 

+  39 

+  32 

+  25 

-i  26 

+  24 

+  22 

+  18 

1  

+  16 

T    15 

•f  22 

+    35 

+  42 

+  40 

+  26 

+  25 

+  28 

+  23 

+  .25 

+  25 

2 

+  16 

f  22 

+  28 

+     34 

4  35 

4-  32 

+  26 

+  27 

+  31 

+  27 

+  23 

+  17 

3 

+  18 

-f-30 

+  35 

+     40 

+  38 

+  38 

+  30 

+  29 

+  31 

+  31 

+  16 

+  16 

4 

+  33 

+  41 

+  48 

+     53 

+  45 

+  38 

+32 

+33 

+81 

+  29 

+  22 

+  26 

5 

-f  44 

+  53 

+  58 

+     67 

+  50 

+  37 

-f  29 

+  32 

+  38 

+  30 

+  26 

+  34 

6 

+  54 

+  68 

->-  74 

+     72 

+  53 

+41 

+32 

+  89 

+  27 

+  40 

+  39 

+  38 

7                -     - 

+  64 

+7« 

+  84 

+   82 

+  57 

+  38 

+  30 

+  25 

+  26 

^43 

+  44 

+  48 

8     .     

+  76 

+  76 

+84 

+   82 

+  51 

+  25 

+  17 

+  19 

-r28 

+  39 

+  48 

+62 

9 

+  72 

+  75 

+  74 

+    68 

+  38 

+  11 

+    4 

+    8 

+  16 

+  36 

+  47 

+  57 

10 

+  49 

+  56 

+  48 

+    46 

+  19 

—   3 

-   9 

rt 

+    3 

!   25 

+  29 

+  35 

11 

+  15 

+  20 

•   22 

+    18 

-    1 

—  19 

—  19 

—  14 

£• 

—    1 

+    7 

+    8 

12      

-23 

-  16 

-  11 

-      9 

—  20 

-29 

-27- 

-24 

-25 

—  21 

-24 

-23 

IS 

-  64 

-  55 

-  50 

-    53 

-  49 

-44 

—  36 

—  39 

—  47 

—  46 

-  53 

—  62 

14 

—  89 

—  80 

—  79 

—    84 

—  71 

—  55 

-42 

-42 

-  50 

-  62 

-70 

—  81 

15 

-92 

—  01 

—  05 

-     99 

-  82 

-65 

—  46 

—42 

—53 

-65 

-67 

-82 

16          

-  71 

—  77 

-  88 

-104 

—88 

—  67 

—48 

—  41 

-  42 

—  49 

-  54 

-  63 

17 

—  45 

—  57 

—  74 

-    % 

—  77 

—  57 

—  41 

—  36 

-34 

—  36 

-37 

—  37 

18                 

—  30 

-  44 

—  64 

—    76 

—  55 

-40 

-  30 

—  26 

-24 

-36 

—  32 

—  26 

19          

—  20 

-  39 

-  57 

—    59 

—  39 

—  24 

—  18 

—  19 

—  20 

—  31 

—  20 

—  13 

20 

—  11 

—  31 

—  45 

—    48 

—  25 

—  11 

—  10 

—  12 

-  14 

—  15 

-    6 

-    4 

21                 -     - 

—    7 

—  23 

-27 

—    25 

—    5 

-r     5 

+    4 

+    1 

+    2 

—    4 

—    2 

+    1 

22                 

Q 

10 

—  12 

+      5 

+  20 

+  28 

+  26 

+  22 

14 

+    7 

+    5 

+    1 

23        

—    5 

—    8 

+    6 

+     22 

+  34 

-t-44 

+87 

+  29 

+  19 

+  17 

+  13 

+    5 

NOTE.— Black-face  type  indicates  maximum  or  minimum. 

•In  Plates  XIII  and  XIV  may  be  seen  the  curves  which  represent  the  thermic  wave 
and  the  principal  semidiurnal  wave  for  several  months. 

These  are  doubtless  the  two  waves  mentioned  by  M.  Angot  in  his  "Traite  elementaire 
de  Meteorologie,"  page  104,  in  his  explanation  of  the  diurnal  variation  of  the  barometer. 

Without  knowing  exactly  of  what  waves  is  being  treated,  it  is  difficult  to  grasp 
thoroughly  this  explanation,  hence  we  have  thought  it  good  to  translate  it  in  this  place ; 
for,  besides  making  it  easier  to  understand  with  the  curves  before  us,  it  will  also  serve  as 
a  discussion  of  our  curves. 

"The  diurnal  variation  of  the  barometer  has  not  as  yet  been  completely  explained.  The 
most  recent  researches  seem  to  indicate  that  it  is  the  result  of  the  superposition  of  two 
distinct  oscillations — one  which  gives  two  maxima  and  two  minima  each  day  and  whose 
period  is  consequently  half  a  day,  and  which  is  called,  therefore,  the  semidiurnal  oscilla- 
tion or  wave;  the  other  which  gives  only  one  maximum  and  one  minimum  in  the  twenty- 
four  hours,  and  which  is  called  the  diurnal  wave. 

"The  two  daily  maxima  of  the  semidiurnal  wave  ought  to  be  equal  to  one  another,  and 
the  two  minima  likewise  equal  to  each  other.  In  all  the  seasons  of  the  year  and  in  all 
parts  of  the  earth  the  two  maxima  take  place  about  10  o'clock,  morning  and  night;  the 
two  minima  about  4  o'clock,  morning  and  afternoon.  The  difference  between  the  maxima 
and  the  minima  of  the  semidiurnal  wave — that  is  to  say,  the  amplitude  of  this  wave  1 
varies  but  little  for  the  same  station  during  the  course  of  a  year;  it  is  greatest  at  the 
equinoxes  and  smallest  at  the  solstices.  This  amplitude  does  not  depend  on  the  topo- 
graphical conditions,  but  solely  on  the  latitude — about  2  mm.  at  the  equator,  it  then 

1  M.  Angot  understands  here  by  the  amplitude  of  the  wave,  the  amplitude  of  a  complete  oscillation ; 
other  authors  and  Angot  himself,  in  the  publication  we  have  quoted  several  times,  takes  as  the  amplitude 
of  the  wave,  the  amplitude  of  a  semioscillation,  or  the  maximum  ordinate  of  the  curve. 
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diminishes  slowly  up  to  the  Tropics,  and  afterwards  very  rapidly  so  that  in  high  latitudes 
it  is  very  small.  It  is  probable  that  this  semidiurnal  wave  is  caused  by  the  action  of 
solar  heat  on  the  whole  mass  of  the  atmosphere,  but  up  to  the  present  time  the  manner 
in  which  it  is  produced  has  not  been  explained  satisfactorily. 

"While  the  semidiurnal  wave  is  the  same  for  all  regions  situated  in  the  same  latitude, 
the  diurnal  wave  has  a  very  variable  amplitude,  which  depends  on  the  topographical  con- 
ditions of  the  region,  and  which  may  be  very  different  even  in  neighboring  stations. 
In  general,  it  is  much  weaker  in  coast  stations  than  in  continental  stations,  and  its 
amplitude  is  in  direct  ratio  to  the  amplitude  of  the  diurnal  variation  of  the  temperature. 
The  hours  of  the  maximum  and  minimum  vary  greatly  according  to  the  countries  and  the 
seasons  of  the  year ;  the  maximum  takes  place  earlier  in  summer  than  in  winter ;  its  mean 
epoch  is  about  9h;  the  minimum,  on  the  contrary,  whose  mean  epoch  is  about  16h,  takes 
place  earlier  in  winter  than  in  summer.  The  period  of  falling,  therefore,  lasts  on  an 
average  only  seven  hours,  while  the  period  of  rising  extends  to  seventeen  hours ;  the  rise 
is  not  -uniform  during  this  time  for  it  is  rapid  immediately  after  the  minimum,  is  slow 
during  the  night,  and  then  increases  rapidly  again  between  sunrise  and  the  maximum. 

"While  the  theory  of  the  semidiurnal  wave  has  not  been  explained  as  yet 1  the  diurnal 
wave,  on  the  contrary,  is  fully  explained  by  the  local  effect  of  the  daily  variation  of  the 
temperature.  During  the  night,  the  lower  layers  of  the  atmosphere  contract  as  they  are 
cooled  and,  in  order  to  fill  the  space  thus  produced  by  the  contraction,  air  from  all  sides, 
especially  from  above,  flows  in.  After  sunrise  the  lower  layers  are  warmed  rapidly  and 
tend  to  expand ;  but  the  heating  is  produced  first  in  the  layers  in  contact  with  the  surface 
of  the  earth,  so  that  the  layers  of  air  between  the  superficies  and  the  higher  layers  which 
have  not  been  heated  as  yet,  do  not  at  first  increase  their  volume  appreciably;  just  as 
happens  when  the  air  in  a  closed  chamber  is  warmed,  the  rise  of  temperature  is  manifested 
principally  by  an  increase  of  pressure  in  the  same  way  the  first  effect  of  the  diurnal 
heating  is  an  increase  of  pressure  in  the  layers  which  are  in  contact  with  the  surface  of 
the  ground,  and  the  barometer  begins  to  rise. 

"A  little  later,  between  7  and  11  o'clock  in  the  morning,  according  to  the  countries  and 
season  of  the  year,  the  heating  of  the  lower  layers  is  sufficient  to  determine  a  general 
ascending  movement  of  the  air  above  the  station;  the  air  expands,  but  as  soon  as  the 
column  of  air  situated  above  the  heated  region  reaches  a  greater  height  than  that  which 
covers  the  neighboring  regions,  the  air  spreads  out  at  the  sides ;  the  quantity  of  air  above 
the  station  diminishes,  therefore,  and  the  barometer  falls.  Toward  the  end  of  the  day, 
when  the  temperature  begins  to  fall  rapidly  in  the  lower  layers,  the  ascending  current  loses 
its  activity,  ceases,  and  is  finally  replaced  during  the  night  by  an  inverse  movement ;  the 
barometer  then  rises,  rapidly  at  first  while  the  temperature  is  falling  rapidly,  and  then 
slowly  till  sunrise." 

It  will  be  noted  that  the  explanation  of  the  thermic  or  diurnal  wave  given  by  M. 
Angot  coincides  in  great  part  with  the  first  theory  of  the  diurnal  variation  of  the  barom- 
eter proposed  by  Raymond,  and  perfected  by  Kreil  and  Renou,  and  which  was  unsatis- 
factory when  it  was  applied  to  the  total  diurnal  variation  of  the  pressure. 

We  may  conclude,  therefore,  that  although  the  analysis  of  the  curves  which  represent 
the  different  meteorological  elements  of  each  region  is  a  long  and  tedious  one,  it  can, 
nevertheless,  be  of  the  greatest  utility,  both  for  the  simultaneous  study  of  various 
phenomena,  the  mutual  relations  of  which  we  wish  to  investigate,  as  well  as  for  the 
deeper  study  of  the  causes  which  may  intervene  in  each  of  them. 

1  We  do  not  know  whether  M.  Angot  had  any  notice  of  the  investigations  of  Doctor  Margules, 
because  neither  here  nor  in  the  work  which  he  published  in  "Annales  du  Bureau  Central  de  France," 
does  he  make  mention  of  them,  notwithstanding  their  prime  importance  in  the  explanation  of  this  wave. 
Doctor  Hann,  in  publications  posterior  to  those  of  Angot,  laments  the  fact  that  the  work  of  Doctor 
Margules  was  so  little  known,  and  that  it  had  not  even  called  the  attention  of  the  "Meteorologische 
Zeitschrift." 


DEDUCTION  OF  BESSEL'S  EQUATIONS  FOR  THE  DETERMINA- 
TION OF  THE  COEFFICIENTS  OF  FORMULA  (I). 


Suppose  that  there  are  any  number  n  of  observations. 
Let  a0,  alt  a,,  etc.,  represent  the  observations  at  0,  1,  2,  etc.,  hours. 
We  saw  on  page  7  that  these  observations  are  represented  by  the  following  equa- 
tions:1 


+P2+PS+  -------- 

cos  x-{-q\  sin  x+pz  cos  2x-\-qz  sin  2x+ps  cos  3#+<?3  sin  3x+ 
a"2=p-f  p\  cos  2z-f</i  sin  2x-\-pz  cos  4o;+<72  sin  4x-\-pa  cos  Gx-\-qa  sin  6*+ 


dn-i=pjrp\  cos  (n— l)x+qisin  (n— I)x+P2cos2(n— 

+pscos3(n— 1)  a; +93  sin  3  (w— 


If  the  number  of  coefficients  which  we  propose  to  determine  be  the  same  as  the 
number  of  observations,  we  would  have  a  system  of  as  many  equations  as  unknowns, 
and  we  could  resolve  it  easily  by  the  ordinary  algebraic  methods.  But  ordinarily  the 
number  of  observations  greatly  exceeds  the  number  of  the  coefficients  we  require;  thus, 
for  example,  in  the  analysis  of  the  barometric  curve,  we  have  twenty-four  observations 
and  we  employ  only  eight  of  these  coefficients. 

In  this  case  the  system  is  more  than  determined;  and  in  the  impossibility  of  finding 
the  values  p  and  q,  which  substituted  in  the  equations,  reproduce  exactly  the  observa- 
tions, we  must  content  ourselves  with  estimating  which  are  the  most  probable  values 
of  p  and  q,  which  give  for  the  second  members  values  compatible  with  the  observations. 

In  the  theory  of  least  squares  it  is  proved  that  the  most  probable  values  are  those 
which  make  the  sum  of  the  squares  of  the  residual  errors  a  minimum.  The  residual 
error  is  the  departure  from  the  observed  value  of  the  value  obtained  in  substituting  in 
the  equations  the  most  probable  values  of  the  unknowns. 

According  to  this  definition,  representing  the  residual  errors  by  the  symbols  v0,  vlt 
vz,  etc.,  and  the  most  probable  values  of  the  coefficients  by  p0,  plt  p2,  etc.,  we  may  write  : 


cos  x+qi  sin  x+pz  cos  2x+qi  sin  2x+ps  cos  Sx+gs  sin  3x+ 
V2=—  «2+p+pi  cos  2x+qi  sin  2x-\-p^  cos  4x-f  92  sin  4x-f/>3  cos  6x+9s  sin  6x4- 
t»3=—  as-fp+pi  cos  Sx-f-gi  sin  3x+p%  cos  6^+92  sin  Qx-\-ps  cos  9x-\-qs  sin  9x+ 


vn-i=— an-i+p+pi  cos  (n— l)x+gi  sin  (n— l)£-j-p2  cos  2(n— \}x-\-qz  sin  2(n— \)x 

+P3  cos  3(n— l)x+qs  sin  3(w— 1 

The  condition  which  makes  the  values  of  p  and  q  most  probable  is  that 

vl-\-v\-\-vl-\-vi-{- —  minimum. 

and  writing  instead  of  vl,  v\,  etc.,  their  values  we  have 

(— ao+p+pi+p2  -I- )2+(— al^-p^  PI  cos  x-\-q\  sin  x-\- )2 

+  (— 02+p+pi  cos  2x+q\  sin  2x+ )z—a  minimum. 

We  know  from  the  differential  calculus  that  when  a  function  is  a  minimum,  its  deriv- 
atives with  respect  to  each  of  the  variables  are  zero.     Hence,  deriving  the  former  ex- 


1  The  quantity  p  which  appears  in  these  equations  represents  the  constant  M  which  we  saw  com- 
pleted BessePs  formula;  it  is  also  to  be  noted  that  the  equations  on  page  7  are  given  here  according  to 

formula  (I). 
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pression  with  respect  to  p,  pif  qlf  etc.,  and  representing  these  derivatives  by  P,  Plt  Ql 
etc.,  we  obtain  the  following  system  — 

P=2(—  aQ-\-p+pi+p2+  ----  )+2(—  ai+p+pi  cos 


sin  x+  ----  ) 

4-2  (—  az+p-{-p\  cos  2x+qi  sin  2x4-  ____  )  +  ____  —  ° 

P\=2(—  ao-\-p+pi-\-p2+  ____  )4-2  cos  x(—  a\+p-\-p\  cos  0:4-91  sin  x+  ____  ) 

4-2  cos  2x(—  az+p+pi  cos  2x+qi  sin  2a-  -•-  ____  )-f  ____  =0 

Qi=2  sin  x(  —  ai+p+pi  cos  x+<7i  sin  x+  ____  )+2  sin  2x(—  a24-j>4-pi  cos  2x+gi  sin  2x4-  ____  )  +  -  __________  =0 

P2=2(—  ao+p+pi-i-pz  ____  )4-2  cos  2x(—  a\-}-p-\-p\  cos  x-\-q\  sin  x4-  ____  ) 

4-2  cos  4x(  —  O2+/>+Pi  cos  2x4-<7i  sin  2x4-  ____  )4-  ______  =0 

Qz—2  sin  2x(—  ai4-/>4-;>i  cos  x4-qri  sin  x4-  ____  )+2  sin  4x(  —  a2+p4-/>i  cos  2cc+gi  sin  2x4-  ____  )4-  ____________  =0 


This  is  the  system  which  must  be  solved  in  order  to  find  the  values  of  p  and  q,  which 
will  be  the  most  probable;  because  in  this  hypothesis  the  preceding  equations  have  been 
obtained. 

After  dividing  by  2  we  can  write  it  in  the  following  form  : 


ip=-< 


—a, 


4-9! 


On  COS  X 

a2  cos  2x 
aa  cos  3x 


1 

+  cos  x 
+  cos  2x 
4-  cos  3x 
+  cos  4x 


4-cos(n-l)x 


+  cosx 
+  cos  2x 
+  cos  3x 


sin  x 
4-  sin  2x 
4-  sin  3x 
4-  sin  4x 


-\-p2 


4-  cos  2x 
4-  cos  4x 
4-  cos  6x 
+  cos  8x 


+  cos  2(w-l)x 


sin  2x 
4-  sin  4x 
4-  sin  6x 
4-  sin  8x 


4-  sin  2  (n-l  )x 


—  a«-icos(w-l)x 


4-  cos2  x 
4-  cos2  2x 
4-  cos2  3x 


4- cos  (n-l )x 

1 

4-  cos  x  cos  2x 
4-  cos  2x  cos  4x 
4-  cos  3x  cos  6x 


4-cos2(n-l)x 


cos  x  sin  x 
4-  cos  2x  sin  2x 
+  cos  3x  sin  3x 


4-  cos  (n-l )  x  cos  2  (n-l )  x 


4-  cos  (n-l )  x  sin  (n-l)  a? 

cos  x  sin  2x 
4-  cos  2x  sin  4x 
+  cos  3x  sin  6x 


+- 


a»  sin  2x 
a3  sin  3x 


— an-isin(n-l)x 


sinx 
+  sin  2x 
4-  sin  3x 


4-  sin  (n-l )x 


sin  x  cos  x 
4-  sin  2x  cos  2x 
4-  sin  3x  cos  3x 


4-  cos  (n-l  )x  sin  2(n-l)x 

sin2  x 
+  sin2  2ar 
4-  sin2  3a- 


+  sin(n-l)x  cos(n-l)x 

sin  x  cos  2x 
+  sin  2x  cos  4x 
4-  sin  3x  cos  6x 


4-  sin(n-l)  cos2(n-l)x 


.=0 


4-  sin2  (n-l  )x 

sin  x  sin  2x 
4-  sin  2x  sin  4x 
4-  sin  3x  sin  6x 


4-  sin  (n-l )x  sin  2  (n-l)  x 
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This  system,  according  to  the  notation  explained  on  page  11,  may  be  written  in 
the  following  manner: 

—  2om+np-f-/>i  2  cosmx+qi  2  sinmx+p?  2  cos  2mo;+<72  2  sin  2mxjr =0' 

—  2amcosmx-rp  2  cosm#+pi  2cos2  mx+qi  2  cos  mx  smmx+pz  2  cos  mx cos 2rax+<72  2  cosmx  sin  2ma;+=0 

—  2omsinmx+?>2  smmx+p\  2  sin  -ma;  cos  mx  -f  91  2  sin2  mx-t-pz'Z  sin -ma;  cos2mx-\-q2  2  sin  ma;  sin  2mx+  =0 

—  2amcos2ma;+p2  cos  2mo;+j0i2  cos  2ma;cosraa;+<7i2  cos  2raojsinma;+p22cos'  2mx 

+92  2  cos  2ma;  sin  2mz-f  =0 

—  2am  sin  2mx-\-p  2  sin  2w»+pi  2  sin  2m#  cos  mx-{-q\  2  sin  2m«  sin  mx+pz  2  sin  2m#  cos  2mx 

+qz  2  sin 


When  nx  (the  product  of  the  number  of  observations  multiplied  by  the  number  of 
degrees  of  the  parts  into  which  the  different  circles  of  reference  are  divided)  is  a  mul- 
tiple of  27T,  the  preceeding  system  is  reduced  to  the  following:1 

— 2om  -}-np  =0 
2  cos2  mx  =0  | 
2  sin2  mx  =0  \  (B) 

cos2  2mx  =0  | 

sin2  2mx  =0  j 


— 2am  cos 
— 2amsin 
— Sam  cos  2mx-\-p2 


But  as  will  be  proved  later 


I  later 

•5  f»o<?2  -wr  —  n 

*H    V^V/O       l/VtAs  ^-v. 

2 

2  sin2  W#=A 

M 

u 

2fV»c2  O'Wy  
vAyo      £j''C/ti(/  ^v 

then  substituting  these  values  in  system  (B)  we  obtain  finally 


=0 


m  cos 
m  sin 
cos 


sin  2mx-{-- 

A 


cos  m# 
sin  ma 


^SamSin  2ma; 


These  are  in  their  general  form  Bessel's  equations  which  we  required;  in  our  par- 
ticular case  of  twenty-four  observations,  ^=12,  which  reduces  these  equations  to  those 

z 

on  page  11.     It  is  also  proved  that  the  constant  p  or  M,  which  completes  Bessel's  for- 
mula, is  the  arithmetical  mean  of  the  observations;  thus  we  have  in  our  case: 

-4023 


v-  24 

It  only  remains  to  prove  the  change  from  (A)  to  (B)  and  the  equations  (C). 
For  the  first  it  is  sufficient  to  prove  that  when  nx  is  a  multiple  of  2w 

2  sin  mx  —  0,  2  cos  mx  =  0,  2  sin  2mx  —  0,  2  cos  2mx  =  0,  etc. 


1  This  will  be  proved  later  on. 
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Let  us  take  the  well  known  trigonometrical  formula: 

cos(z—y)—cos(z-\-y)  ,.,,. 

2   sm  z  sm  y  =  cos(z— y)  —  cos(z+y) .  .  sm  z=—  '   .  -   (1) 

£  sin  // 

By  a  series  of  applications  of  Formula  (1)  in  which  we  write  x,  2x,  3x 


(n—  1)#  instead  of  z  and^  instead  of  y,  we  have: 

a 

cos|-cos3| 

£4  £, 

-- 
2sin 


sin  2x= 


2  slnf 

5x  Ix 


2sin| 


cos*""'  ~7~—  cos" 
sin(n-2)#=— 


v_.,  „)#  (2n-l)x 

cos- — 5^    —  cos- — =-^ 

—  ^ 

2  sin| 

and  taking  the  sum  we  obtain: 

a;          (2n— l)x 
cos  o— cos  - — n-1 

z  ^ 


sin  #sin  2irsin  3x--  ............    -sin(w-l)x=2  sin  mx= 


cos  ^—  cos  na;  cos^—  sin 

MA 


2  sin 


2sin| 

if  n#  is  a  multiple  of  2*,.  we  know  that 

cos  nx=\;  sin  no;=0 
and  hence 

cos|-cos| 

2  sin  mx=  —  =0 

2sin| 

In  like  manner  it  is  proved  that  2  sin  2mx  —  0,  2  sin  3m#  —  0,  etc.,  so  for  this  last 

we  will  write  in  Formula  (1)  3x,  Gx,  9#,  ,  (n—l)  3x  instead  of  z,  and  „    instead 

of  y,  which  will  give: 

•  ( Qw, 3 )  x  ^x  ^x  ^r 

— cos  -  cos   _  —  cos  3nx  cos  -^ — sin  3nx  sin  0 

£t  &  £  £ 

2  sin??  2  sin3^ 

n  o 

and  as  3w#  is  a  multiple  of  2*,  it  follows  that  2  sin  3mz  =  0. 
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To  prove  that  2  cos  mx  —  0,  2  cos  2mx  =  0,  we  make  use  of  the  trigonometrical  for- 
mula 

sin  (z+y)  —sin  (z—y~)  /n\ 
2cos  z  sin  y=sm(z+y)—sm(z—y) .  .cos  z=  —        '  0   .  -(2) 


X 

and  writing  0,  x,  2x, (n—l)x  instead  of  z,  and  -  instead  of  y,  we  have 

sin  |+sin  ~ 

cos  0=— 

2  sin? 


cos  x— 


2sin| 
sin  f  -sin 


cos  3x= 


2  sin  | 


cos(n—  2)a;=  — 


2  sin 


cos(?i—  l)x=  — 


o     • 

2sm2 


Taking  the  sum  we  obtain : 
1+cos 


2  sin 


sin^+sin  wo:  cos ^— cos  nx  sin- 
"  2sin| 

and  remembering  that  nx—2-rr,  it  follows  that: 

x         x 
sin^-sm- 

2  cos  mx=—  —  =0 

2  sin  I 

127182 3 
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This  having  been  proved,  it  is  readily  seen  that  the  other  terms  of  system   (A)   of 
the  form 

2  cos  mx  sin  mx;  2  sin  mx  cos  2mx;  2  sin  mx  sin  2mx;  2  cos  mx  sin  2mx;  2  cos  mx  cos  2mx 

are  also  reduced  to  zero;  because  as  the  product  of  the  sine  and  cosine  of  a  same  arc 
is  equal  to  the  half  of  the  sine  of  the  double  arc,  we  obtain: 

2  cos  mx  sin  mx  =  ±  2  sin  2mx  —  0 

the  same  can  be  proved  for  the  other  four  expressions  by  making  use  of  the  well-known 
trigonometrical  formulas 

sin  x  cos  y=  i  sin(#+  2/)+|  sm(x-y) 
cos  x  sin  y  =  i  sin  (x-\-y)  —  -i-  sin  (x— y) 
cos  x  cos  ?/=i  cos  (#+?/)  +  |  cos(o;— y) 
sin  #  sin  i/  =  ^  cos  (#— ?/)  — ^  cos  (#+?/) 

for  writing  mx  and  2m#  for  #  and  y,  we  have: 

2  sin  mx  cos  2m£=^2  sin  3mx— ^2  sin  m# 
2  cos  ma:  sin  2m#=^2  sin  3m#+^2  sin  ra# 
2  cos  mx  cos  2m#=i2  cos  3m#+^2  cos  mx 
2  sin  ma;  sin  2mx=^  cos  mx— i2  cos  3m# 

and  as  was  proved  before,  the  second  members  are  zero. 

We  now  proceed  to  prove  the  equations   (C)  ;  a  simple  trigonometrical  expansion 
gives  us  that: 

2  cos2  x  =  1  +  cos  2x 
and  according  to  this  we  have: 

2cos20  =  1  +  cos  0 

2  cos2  x  —  1+cos  2x 
2  cos2  2#  =  1+cos  4z 
2  cos2  Bx  =?=  1+cos  Qx 


2  cos2(n—l)x  =  1+cos  2(n— l)x 
taking  the  sum  we  get 

2  2  cos2  mx=n-\-  (1+cos  2#+cos  4#+cos  6x-\- +cos  2(n— l)x  —  n+2  cos  2mx 

but,  as  we  proved  before  2  cos  2mx  =  0,  hence 

2  cos2  mx  =  ™ 

LA 

In  an  analogous  manner  this  same  value  can  be  proved  for  any  term  of  the  form 
2  cos2  mkx. 

fV\ 

To  prove  that  2  sin2  mx  —  ~  we  make  use  of  the  known  equation 

sin2  x  =  1— cos2  x 
and  writing 

sin2  0  =  1— cos2  0 

sin2  x  —  1— cos2  x 
sin2  2x  =  1— cos2  2a: 
sin2  3x  —  1— cos2  3x 

O  i  n  2  §  /vi     1     |   /y»    "I    /*/^Q2  I  /M I    1   /v» 

Olll      y  /(/ J_  }  *fj    J. ^VJO      y^  Iv X  /  *v 
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and  adding,  we  have 


2  sin2  mx  —  n—^  cos2  mx 


n 


but  we  have  just  proved  that  2  cos2  mx  —  _-,  therefore 

*- 

/Yl>      /Y\J 

2  sin2  mx  =  n—~=9 

u         u 

We  do  not  think  it  necessary  to  delay  any  longer  over  this  work  and  give  the  whole 
of  Bessel's  method;  those  who  may  be  interested  in  it  may  consult  Bessel's  own  expla- 
nation, translated  into  English  by  Mr.  Robert  H.  Scott,  and  published  as  an  appendix  to 
the  Quarterly  Weather  Report  of  the  Meteorological  Office,  London,  1870. 

TABLE  I. — Mean  hourly  values  of  atmospheric  pressure  at  Manila  for  each  month  (1890-1909). 


Hour. 

January. 

Febru- 
ary. 

March. 

April. 

May. 

June.    July. 

August. 

Septem- 
ber. 

,**.. 

Novem- 
ber. 

Decem- 
ber. 

0 

61.  382 

61.  704 

60  924 

59  938 

58  872 

58  5?0    57  880 

57  834 

57  934 

58  824 

59  684 

60  880 

1  

61.092 

61.404 

60  556 

59  624 

58  601 

58  193    57  468 

57  459 

57  549 

58  408 

59  330 

60  590 

2 

60  722 

61  072 

60  214 

59  228 

58  187 

57  789    57  128 

57  102 

57  1% 

58  080 

58  970 

60  162 

3 

60.464 

60  854 

59  982 

59  029 

57  982 

57  620    56  936 

56  868 

56  944 

57  883 

58  686 

59  930 

4 

60.  516 

60.833 

59  989 

59  071 

57  991 

57  556    56  870 

56  820 

56  885 

57  842 

58  726 

59  967 

5  

60.718 

61.  047 

60  204 

59  346 

58  168 

57  649    56  944 

56  926 

57  046 

58  030 

58  953 

60  176 

6 

61  069 

61  470 

60  661 

59  704 

58  486 

57  950    57  222 

57  191 

57  306 

58  476 

59  388 

60  512 

7 

61.  546 

61  943 

61  172 

60  193 

58  885 

58  258    57  546 

57  518 

57  705 

58  915 

59  836 

60  966 

8 

62.  028 

62.  342 

61  564 

60  584 

59  166 

58  446    57  762 

57  830 

58  117 

59  248 

60  218 

61  467 

62.264 

62.630 

61  762 

60  693 

59  274 

58  539    57  869 

57  980 

58  246 

59  452 

60'  423 

61  642 

10 

62  128 

62  568 

61  630 

60  559 

59  142 

58  456    57  819 

57  918 

58  182 

59  366 

60  274 

61  480 

11 

61.  701 

62  120 

61  250 

60  143 

58  816 

58  188    57  622 

57  734 

57  930 

58  928 

59  877 

61  084 

12  .. 

61.  057 

61  494 

60  620 

59  562 

58  333 

57  824    57  283 

57  354 

57  428 

58  349 

59  258 

60  485 

13.  .. 

60.  282 

60  712 

59  827 

58  732 

57  686 

57  344    56  846 

56  828 

56  798 

57  732 

58  572 

59  735 

14 

59.656 

60  058 

59  134 

58  033 

57  118 

56  910    56  444 

56  420 

56  378 

57  204 

58  071 

59  198 

15 

59.  360 

59  650 

58  672 

57  615 

56  771 

56  576    56  169 

56  163 

56  087 

56  944 

57  879 

58  952 

16... 

59.  462 

59  666 

58  620 

57  490 

56  656 

56  4%  l   56  066 

56  087 

56  128 

57  078 

57  993 

59  082 

17 

59  814 

69  954 

58  867 

57  700 

56  890 

56  706    56  232 

56  254 

56  368 

57  391 

58  344 

59  480 

18 

60.228 

60  352 

59  266 

58  204 

57  402 

57  130    56  598 

56  640 

56  788 

57  732 

58  708 

59  882 

19  . 

60.  698 

60  798 

59  749 

58  762 

57  917 

57  619    57  062 

57  087 

57  240 

58  1% 

59  218 

60  369 

20 

61  158 

61  274 

60  261 

59  266 

58  397 

58  078    57  488 

57  528 

57  688 

58  728 

59  698 

21 

61  461 

61  654 

60  742 

59  753 

58  835 

58  474    57  861 

57  924 

58  101 

59  071 

59  964 

61  094 

22.-- 

61.544 

61  839 

61  012 

60  124 

59  144 

58  755    58  158 

58  218 

58  294 

59  208 

60  054 

61  156 

23._. 

61.496 

61  852 

61  078 

60  170 

59  157 

58  811    58  172 

58  177 

58  180 

59  130 

59  958 

61  059 

24 

61  365 

61  717 

60  898 

59  902 

58  852 

58  498    57  868 

57  852 

57  921 

58  867 

59  736 

60  900 

1 

TABLE  II. — Same  as  Table  I  corrected  for  noncyclic  variation. 


Hour. 

January. 

Feb- 
ruary. 

March. 

April. 

May. 

June. 

July. 

August. 

Septem- 
ber. 

October. 

Novem- 
ber. 

Decem- 
ber. 

0- 

61.374 

61.  710 

60  911 

59  920 

58  862 

58  509 

57  874 

57  843 

57  928 

58  846 

59  710 

60  890 

1-.. 

61.  084 

61.410 

60.544 

59  608 

58  592 

58  183 

57  463 

57  467 

57  543 

58  428 

59  354 

60  599 

2 

60.715 

61.077 

60  203 

59  213 

58  179 

57  780 

57  123 

57  110 

57  185 

58  098 

58  992 

60  170 

3 

60.  458 

60.859 

59  972 

59  015 

57  974 

57  612 

56  932 

56  875 

56  939 

57  899 

58  706 

59  938 

4 

60.  510 

60.  837 

59  980 

59  059 

57  984 

57  549 

56  866 

56  826 

56  881 

57  856 

58  743 

59  974 

5  

60.  713 

61.  051 

60.196 

59  336 

58  162 

57  643 

56  940 

56  931 

57  042 

58  043 

58  968 

60  182 

6 

61.065 

61.  473 

60  655 

59  695 

58  481 

57  944 

57  219 

57  186 

57  303 

58  487 

59  401 

60  517 

7 

61.  542 

61.  946 

61  167 

60  186 

58  881 

58  253 

57  544 

57  522 

57  702 

58  924 

59  847 

60  970 

8 

62.  025 

62.  344 

61.561 

60  578 

59  164 

58  443 

57  760 

57  833 

58  115 

59  255 

60  227 

61  469 

9 

62  262 

62.632 

61  759 

60  688 

59  272 

58  536 

57  867 

57  982 

58  244 

59  457 

60  430 

61  644 

10. 

62.  127 

62.569 

61  628 

60  556 

59  140 

58  454 

57  818 

57  919 

58  181 

59  370 

60  278 

61  482 

11.. 

61.700 

62.120 

61.  249 

60  141 

58  815 

58  187 

57  621 

57  735 

57  929 

58  930 

59  879 

61  085 

12 

61  057 

61  494 

60  620 

59  562 

58  333 

57  824 

57  283 

57  354 

57  428 

58  349 

59  258 

60  485 

13 

60.283 

60.  712 

59  828 

58  734 

57  687 

57  345 

56  847 

56  827 

56  799 

57  730 

58  570 

59  734 

14... 

59.  657 

60.057 

59.  136 

58  036 

57  120 

56  912 

56  445 

56  419 

56  379 

57  200 

58  067 

59  196 

15 

59  362 

59  648 

58  675 

57  620 

56  774 

56  579 

56  171 

56  161 

56  089 

56  939 

57  872 

58  950 

16 

59  465 

59.664 

58  623 

57  496 

56  658 

56  499 

56  068 

56  084 

56  130 

57  071 

57  984 

59  080 

17..- 

59.  818 

59.  951 

58  872 

57  707 

56  894 

56  711 

56  234 

56  250 

56  371 

57  382 

58  333 

59  476 

18_-_ 

60.232 

60.  349 

59.  272 

58  213 

57  407 

57  136 

56  601 

56  636 

56  791 

57  721 

58  695 

59  877 

19 

60.703 

60.794 

59  757 

58  772 

57  923 

57  625 

57  066 

57  082 

57  244 

58  183 

59  203 

60  363 

20... 

61.  164 

61.  270 

60  270 

59  278 

58  404 

58  085 

57  492 

57  522 

57  692 

58  714 

59  681 

60  804 

21 

61  467 

61  649 

60  752 

59  767 

58  843 

58  482 

57  866 

57  917 

58  106 

59  055 

59  944 

61  086 

22 

61  551 

61.834 

61  023 

60  139 

59  152 

58  764 

58  163 

58  210 

58  299 

59  190 

60  032 

61  148 

23... 

61.504 

61.846 

61  090 

60  186 

59  166 

58  821 

58  177 

58  169 

58  186 

59  110 

59  934 

61  050 

24 

61  374 

61  710 

60  911 

59  920 

58  862 

58  609 

57  874 

57  843 

57  928 

58  846 

59  710 

60  890 

Mean 

60  910 

61  221 

60  323 

59  313 

58  244 

57  828 

57  227 

57  244 

57  354 

58  343 

59  254 

60  424 
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TABLE  III. — Diurnal  variation  of  atmospheric  pressure  at  Manila  (1890-1909), -  corrected  for  noncyclic 

variation  (unit,  0.01). 


Hour. 

January. 

Febru- 
ary. 

March. 

April. 

May. 

June. 

July. 

August. 

Septem- 
ber. 

October. 

Novem- 
ber. 

Decem- 
ber. 

0.. 

+  46 

+  49 

+  59 

+  61 

+  62 

+  68 

+  65 

+  60 

+  57 

+  50 

-f  46 

+  47 

1 

+  17 

+  19 

—  22 

+  30 

+  35 

+  36 

+  24 

+  22 

+  19 

+     8 

+  10 

+  18 

2 

—  20 

—  14 

—  12 

—  10 

—    6 

—    5 

—  10 

—  13 

—  17 

—  25 

—  26 

—  25 

3.-       . 

—  45 

—  36 

—  35 

—  30 

—  27 

—  22 

—  30 

—  37 

—  42 

—  44 

—  55 

—  49 

4  

—  40 

—  38 

—  34 

—  25 

—  26 

—  28 

—  36 

—  42 

—  47 

—  49 

—  51 

—  45 

5 

—  20 

—  17 

—  13 

+     2 

—    8 

—  18 

—  29 

—  31 

—  31 

—  30 

—  29 

—  24 

6— 

+  16 

+  25 

+  33 

+  38 

+  24 

+  12 

—    1 

—    6 

—    5 

+  14 

+  15 

+    9 

7  

+  63 

+  72 

+  84 

+  87 

+  64 

+  42 

+  32 

+  28 

+  35 

+  58 

+  59 

+  55 

8 

+  112 

+112 

-t-124 

+  126 

+  92 

+  62 

+  53 

+  59 

+  76 

+  91 

+  97 

+  104 

9 

+  135 

+  141 

+  144 

+  138 

+  103 

+  71 

+  64 

+  74 

+  89 

+  11] 

+  118 

+122 

10. 

+122 

+  135 

+  130 

+124 

+  90 

+  63 

+  59 

+  68 

+  83 

+103 

+102 

+106 

11— 

+  79 

+  90 

+  93 

1    83 

+  57 

+  36 

+  39 

+  49 

+  58 

+  59 

+  62 

+  66 

12.... 

+  15 

+  27 

+  30 

+  25 

+     9 

0 

+     6 

+  11 

+     7 

-     5 

0 

+     6 

13— 

—  63 

—  51 

—  50 

—  58 

—  56 

—  48 

—  38 

—  42 

—  56 

—  61 

—  68 

—  69 

14 

—125 

—  116 

—119 

—128 

—112 

—  92 

—  78 

—  82 

—  98 

—114 

—119 

—123 

15— 

—155 

—  157 

—165 

—169 

—147 

—125 

—106 

—108 

—126 

—140 

—138 

—147 

16.... 

—144 

—  156 

—170 

—182 

—159 

—133 

—116 

—116 

—122 

—127 

—127 

—134 

17— 

—109 

—127 

—145 

—161 

—135 

—112 

—  99 

—  99 

—  98 

—  % 

—  92 

—  95 

18. 

—  68 

—  87 

—  105 

—110 

—  84 

—  69 

—  63 

—  61 

—  56 

—  62 

—  56 

—  55 

19.. 

—  21 

—  43 

—  57- 

—  54 

—  32 

—  20 

—  16 

—  16 

—  11 

—  16 

—    5 

—    6 

20.... 

+  25 

+    5 

—     5 

—    4 

+  16 

+  26 

+  26 

+  28 

+  34 

+  37 

+  43 

+  38 

21 

+  56 

+  43 

+  43 

+  45 

+  60 

+  65 

+  64 

+  67 

+  75 

+  71 

+  69 

+  66 

22. 

+  64 

+  61 

+  70 

+  83 

+  91 

+  94 

+  94 

+  97 

+  94 

-  85 

+  78 

+  72 

23— 

+  59 

+  62 

+   77 

+  87 

+  92 

+  99 

+  95 

+  92 

+  83 

+  77 

+  68 

+  63 

TABLE  IV. — Constants  of  the  Harmonic  Formula  (I)  computed  from  Table  III. 


Month. 

Pi- 

flfl- 

PZ- 

<72- 

PS- 

9  3- 

P4- 

9  4- 

January  

4  0.  1341 

+0.  5476 

+0.25% 

—0.9431 

+0.0258 

+0.  1340 

+0.  0196 

—0.  0166 

February  

+  .0891 

+  .65% 

+  .3438 

—  .9215 

+  .0233 

+  .1074 

+  .0258 

—  .0187 

March  . 

+     1208 

+     7495 

+  .3949 

—  .9530 

+  .0395 

+  .0512 

+  .0362 

—  .0161 

April  ... 

+  .2022 

+     7948 

+  .3888 

—  .%61 

+  .0216 

+  .0162 

+  .0383 

—  .0158 

May  

+     2980 

+  .5711 

+  .3396 

—  .8508 

+  .0073 

—  .0048 

-1-  .0313 

—  .0007 

June  

+  .3522 

+  .3860 

+  .3259 

—  .7226 

+  .0098 

—  .0361 

+  .0300 

0 

July  . 

+     2895 

+    2829 

+    3371 

—  .6816 

+    0168 

—  .0453 

+   .0192 

—  .0101 

August  . 

+     2579 

+    2783 

+  .3458 

—  .7282 

+  .0162 

—  .0201 

+  .0083 

—  .0245 

September  

+     2475 

+     3152 

+  .3051 

—  .8257 

+  .0194 

+  .0364 

+  .0041 

—  .0274 

October  .  . 

+     1994 

+  .4065 

+  .2465 

—  .9080 

+  .0404 

+  .0535 

+  .0066 

—  .0317 

November  

+  .1841 

+  .3970 

+  .2030 

—  .9236 

+  .0500 

+  .0836 

+  .0125 

—  .0144 

December 

+     1893 

+     4364 

+    2169 

—    9307 

+     0338 

+  .1225 

+  .0188 

—  .0122 

TABLE  V. — Amplitudes  and  phase  angles  of  Fourier  "waves"  computed  from  Table  IV,  using  L.  M.  T. 


Month. 

wt. 

uz. 

W3. 

W4. 

uv 

t/2- 

t/3. 

u*. 

January  

0.5638 

0.9782 

0.1350 

0.0257 

O          ' 

13    46 

o        / 
164    37 

o        / 
11       1 

O           ' 

130    17 

February                                                      _, 

6607 

9836 

.1100 

.0319 

3    20 

159    32 

12    14 

125    58 

March  

.7592 

1.0315 

.0647 

.0392 

9      9 

157    30 

37    40 

112    34 

April  

.8202 

1.0414 

.0271 

.0415 

14    17 

158      5 

53      8 

112    25 

May  

.6442 

.9161 

.0088 

.0313 

27    33 

158    14 

123    42 

91    22 

une  

.5231 

.7927 

.0374 

.0300 

42    20 

155    44 

164    49 

90      0 

July  

.4047 

.7604 

.0483 

.0217 

45    39 

153    41 

159    39 

117    51 

August.  

.3794 

.8062 

.0259 

.0259 

42    49 

154    36 

141    17 

161    15 

September  

.4007 

.8802 

.0412 

.0278 

38      9 

159    43 

28      6 

171    26 

October  . 

4528 

9409 

.0670 

.0324 

26      8 

164    49 

37      7 

168    12 

November  

.4376 

.9457 

.0974 

.0600 

24    53 

167    36 

30    53 

167    59 

December  .     . 

.4757 

.9557 

.1271 

.0707 

23    27 

166    53 

15    26 

164    37 

HARMONIC   FORMULA   OF   FOURIER   AND   BESSEL. 
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JANUARY. 

FEBRUARY. 

Hour. 

Fourier  "waves." 

Result- 
ant. 

Depart- 
ures 
from  ob- 
served 
values. 

Hour. 

Fourier  "waves." 

Result- 
ant. 

Depart- 
ures 
from  ob- 
served 
values. 

First. 

Second. 

Third. 

Fourth. 

First. 

Second. 

Third. 

Fourth. 

0 

-rO.  134 
•f  .271 
+  .390 
+  .482 
+  .541 
-  .564 

-0.  259 
—  .246 
I  —  .687 
—  .943 
:  —  .947 
—  .6% 
—  .259 
+  .246 
i-  .687 
+  .943 
+  .947 
1   T  .696 
+  .259 

:  o  026 

+0.020 
—  .004 
—  .024 
—  .020 
+  .004 
-  .024 
+  .020 

+0.439 
+  .133 
—  .188 
-  .406 
—  .428 
-  .220 
+  .176 
+  .662 
+  1.096 
+  1.328 
+1.242 
+  .807 

-0.02 
—  .04 
-  .01 

;    +  .04 

—  .03 
—  .02 
+  .02 
-  .03 
—  .02 
—  .02 
+  .02 
+  .02 

0  - 

+0.038 
+  .208 
-  .364 
+  .494 
+  .591 
+  .647 
+  .660 
+  .627 

+0.344 
—  .163 
—  .627 
—  .922 
—  .970 
—  .759 
-  .344 

+0.023 
+  .093 
+  .107 
+  .060 
-  .023 

+0.026 
—  .003 
—  .029 
—  .026 

+0.431 
+  .135 
—  .185 
—  .394 
—  .399 
—  .176 
+  .235 
+  .727 
+  1.174 
+  1.429 
+  1.377 
+  .982 
-  .309 
—  .467 
—1.127 
—1.502 
—1.535 
—1.284 
—  .871 
—  .407 
+  .022 
+  .363 
+  .569 
+  .594 

—0.06 
—  .05 
*-  .04 
—  .03 
—  .02 
—  .01 
-  .01 
+  .01 
+  .05 
+  .02 
+  .03 
+  .08 
+  .04 
+  .04 
+  .03 
+  .07 
+  .02 
—  .01 
0 
+  .02 
—  .03 
—  .07 
—  .04 
—  .03 

1 

H    .112 
+  .133 
+  .075 
—  .026 
—  .112 
—  .133 
—  .075 
+  .026 

1  

2 

2 

3 

3  . 

4 

4 

5 

5  

6. 

+  .548 
+  .495 
-    .407 
+  .293 
+  .158 
+  .012 
—  .134 

6  

7 

7  

8 

8  

+  .553 

9 

9  

+  .440 

10... 

10  

+  .297 

11 

11... 

+  .134 

12 

+  .119 
—  .633 
—  1.  234 
—1.  520 
—1.458 
—1.124 
—  .654 
—  .178 
+  .230 
+  .518 
+  .660 
+  .633 

—  .03 
0 
+  .02 
+  .03 
—  .02 
—  .03 
+  .03 
+  .03 
-  .02 
—  .04 
+  .02 
+  .04 

12... 

—  .038 

13.. 

—  .271 
—  .390 
—  .482 
—  .541 
—  .564 
—  .548 
—  .495 
—  .407 
—    293 

—  .246 
—  .687 
—  .943 

13  

14 

14  

15 

15  

16.. 

i  —  .947 
—  .696 
—  .259 

;  +  .246 

+  .687 
+  .943 

16  

17 

17— 

18 

18..-  

19 

19  

20... 

20  

21 

21- 

22. 

—  .158 
—  .012 

+  .947 
+  .696 

22  

23  

23  

MARCH. 

APRIL. 

0 

+0.  121      -f  0.  395 
+  .  310  •  —  .  135 
+  .479     —  .628 
+  .616     —  .954 
+  .710     —1.023 
+  .755  '  —  .818 
+  .749     -  .395 
-r  .692 

+0.040 
+  .064 
+  .051 
+  .008 
—  .040 

+0.036 
+  .005 
—  .031 
—  .036 

+0.592 
+  .234 
—  .129 
-  .366 
—  .358 
—  .096 
-r  .339 
+  .824 
+  1.226 
+1.427 
+1.339 
+  .934 
+  .270 
—  .504 
—1.  189 
—1.  614 
-1.698 
—1.478 
—1.057 
—  .544 
-  .032 
+  .409 
+  .697 
+  .764 

0 
+0.01 
—  .01 
—  .02 
—  .02 
+  .03 
+  .01 
—  .02 
-  .01 
-  .01 
+  .04 
0 
—  .03 
0 
0 
+  .04 
0 
—  .03 
—  .01 
+  .03 
+  .02 
—  .02 
0 
.01 

0.-. 

+0.203 
+  .401 
+  .572 
+  .705 
+  .790 
+  .820 
+  .795 
+  .715 

+0.388 
—  .147 
—  .642 

f\CC 

.  *roO 

—1.031 
—  .819 
—  .388 

+0.022 
+  .027 
+  .016 
—  .004 
—  .022 

+0.039 
+  .006 
—  .033 
—  .039 

+0.652 
+  .287 
—  .067 
—  .304 
-  .269 
+  .007 
+  .430 
+  .872 
+  1.218 
+  1.373 
+1.263 
+  .859 
+  .202 
-  .569 
-1.263 
—1.706 
—1.805 
—1.  579 
—1.128 
-  .566 
0 
+  .481 
+  .787 
+  .845 

+.  04 
—  .01 
+  .03 
0 
—  .02 
-  .01 
+  .03 
0 
—  .04 
—  .01 
+  .02 
+  .03 
—  .05 
+  .01 
+  .02 
—  .02 
+  .01 
+  .03 
—  .03 
-  .03 
+  .04 
+  .03 
—  .04 
—  .01 

1  

1— 

2 

2 

3 

3 

4  

4  

5  

5— 

6 

6 

7 

7 

8  

-  .589 

8— 

+  .587 

9  

+  .445 

9_—  

+  .419 

10 

+  .270 

10 

+  .222 

11. 

+  .077 

11 

+  .011 

12- 

—  .121 

12— 

—  .203 

13... 

13 

14 

14 

15 

15  

16  

16-.-  

17 

17 

18 

18 

19  

19  

20— 

20  

21 

21 

22 

22 

23 

23 

; 

MAY. 

JUNE. 

0 

+0.298 
+  .435 
+  .544 
+  .614 
+  .643 
+  .629 
+  .571 
+  .475 

+0.340 
—  .131 
—  .567 
—  .851 
—  .907 
—    720 

+0.007 
+  .002 
—  .005 
—  .009 
—  .007 

+0.031 
+  .015 
—  .016 
—  .031 

+0.  676 
+  .321 
—  .044 
—  .277 
—  .286 
—  .077 
+  .267 
+  .630 
+  .904 
+1.015 
+  .915 
+  .587 
+  .066 
—  .553 
—1.122 
-1.487 
—1.558 
—1.  331 
—  .885 
—  .338 
+  .198 
+  .625 
+  .869 
•+•  .885 

+0.06 
-  .03 
+  .02 
—  .01 
-  .03 
0 
+  .03 
—  .01 
-  .02 
—  .01 
+  .02 
+  .02 
—  .02 
+  .01 
0 
-  .02 
+  .03 
+  .02 
—  .04 
—  .02 
+   .04 
+  .03 
—  .04 
-  .03 

0    _ 

+0.352 
+  .440 
+  .498 
+  .522 
+  .511 
+  .464 
+  .386 
+  .282 

+0.326 
—  .079 
—  .463 
—  .723 
—  .789 
—  .644 
—  .326 

+0.  010 
—  .018 
—  .036 
-  .032 
-  .010 

+0.030 
+  .015 
—  .015 
—  .030 

+0.  718 
+  .358 
—  .016 
—  .263 
—  .303 
—  .177 
+  .126 
+  .408 
+  .616 
+  .700 
+  .626 
+  .387 
—  .006 
—  .486 
—  .940 
—1.243 
—1.305 
-1.111 
-  .718 
—  .220 
+  .280 
+  .686 
+  .922 
+  .931 

+0.04 
0 
+  .03 
—  .04 
—  .02 
0 
+  .01 
-  .01 
0 
-  .01 
0 
+  .03 
—  .01 
—  .01 
-  .02 
+  .01 
+  .02 
+  .01 
—  .03 
—  .02 
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HARMONIC   FORMULA   OF   FOURIER   AND   BESSEL. 
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INTRODUCTION. 


Juzgo  que  el  presente  folleto  "The  Harmonic  Formula  of  Fourier  and  Bessel,  and 
Its  Aplication  to  the  Study  of  the  Diurnal  Variation  of  the  Atmospheric  Pressure 
in  Manila  During  the  Period  1890-1909"  preparado  por  el  P.  Antonio  Galan,  S.  J., 
agregado  a  la  direccion  de  este  Observatorio  de  Manila,  ha  de  ser  bien  recibido  por 
los  meteorologos,  los  cuales  son  de  parecer  que  semejante  discusion  o  estudio  deberia 
hacerse  en  todos  los  observatories  principales,  en  donde  se  haya  acumulado  material 
suficiente  en  regular  mimero  de  anos,  para  ser  tratado  matematicamente  con  proba- 
bilidades  de  resultado  cientifico  aceptable. 

Ojala  que  pronto,  como  espero,  se  pueda  encontrar  tiempo  y  oportunidad,  para 
tratar  de  semejante  manera  como  ha  hecho  ahora  el  P.  Galan  con  los  datos  barome- 
tricos,  los  reunidos  en  este  Observatorio  sobre  la  temperatura,  humedad  y  otros  elementos 
de  meteorologia  periodicos. 

Ojala  tambien  que  la  estancia  en  el  Observatorio  de  Manila  del  P.  Galan,  destinado 
a  continuar  la  benemerita  labor  de  los  PP.  Vines,  Gangoiti  y  Gutierrez  Lanza  del 
Observatorio  de  Belen  de  la  Habana,  le  sea  tan  provechosa  como  para  mi  lo  fue  mi 
permanencia  en  la  Habana  los  primeros  meses  del  ano  1892. 

JOSE  ALGUE, 
Director  de  la  Oficina  Meteorologica  de  FiUpinas. 

OBSERVATORIO  DE  MANILA,  19  de  Junto  de  1914. 
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FORMULA  ARMONICA  DE  FOURIER  Y  DE  BESSEL. 


La  formula  armonica  de  Fourier  y  de  Bessel  suele  presentarse  en  una  de  las  dos 
f ormas  siguientes : 

(I)  y=M+pl  cos  x+q^  sen  x+pz  cos  2x+q2  sen  2x-\-p3  cos  3x-\-qs  sen  3x-\- 

(II)  y=M+ul  sen  (Ui+x)+uz  sen  (Uz+2x)  +u3  sen  (£73+3z)  + 

Facilmente  se  ve  que  estas  dos  formulas  son  una  misma  bajo  distinta  forma,  pues 
desarrollando  la  (II)  se  tiene: 

y=M+u1  sen  U1  cos  x-\-u±  cos  Z7,  sen  x-\-u2  sen  U2  cos  2x+u2  cos  C72  sen  2#-|- 

y  haciendo  en  esta  expresion 

ul  sen  Ul=p1  u1  cos  U1=qi 

u2  sen  U2=p2  u2  cos  U2—q2 

etc.  etc. 

la  (II)  queda  reducida  a  la  (I). 

Las  aplicaciones  de  esta  formula  son  muy  numerosas  e  importantes. 

El  teorema  de  Fourier  (dicen  Thomson  y  Tait)  ademas  de  ser  uno  de  los  mas  bellos  resultados 
del  analisis  moderno,  puede  decirse  que  es  instruments  indispensable  para  abordar  casi  todas  las 
reconditas  cuestiones  de  la  fisica  moderna.  Decir  que  las  vibraciones  sonoras,  la  propagacion  de  las 
senales  electricas  a  lo  largo  del  hilo  del  telegrafo  y  la  conduccion  del  calor  por  la  corteza  terrestre 
son  cuestiones  en  su  generalidad  imposibles  de  discutir  sin  su  auxilio,  es  dar  unicamente  una  debil 
idea  de  su  importancia. 

En  el  presente  trabajo  nos  limitaremos  a  la  aplicacion  que  de  la  formula  se  hace 
en  meteorologia ;  procuraremos  declarar  lo  que  cada  uno  de  sus  elementos  representa, 
y  haremos  algunas  indicaciones  sobre  su  utilidad  y  el  modo  practice  de  servirse  de  ella.1 

La  mayor  parte  de  los  fenomenos  que  se  estudian  en  meteorologia  son  periodicos ; 
las  curvas  que  los  representan  son,  por  lo  tanto,  curvas  periodicas,  es  decir,  curvas  que 
representan  el  movimiento  de  un  punto,  que  recorre  un  numero  indeterminado  de  veces 
el  mismo  camino  en  el  mismo  tiempo. 

Analogamente  a  lo  que  Helmholtz  demuestra  en  acustica,  que  un  mismo  sonido 
compuesto  se  resuelve  siempre  en  los  mismos  sonidos  elementales,  probo  Fourier  que 
cualquiera  curva  periodica  se  puede  descomponer  en  una  serie  de  curvas  armonicas  de 
periodos,  1,  ±,  ^,  etc.,  de  la  curva  dada;  y  que  solo  es  posible  una  combination  de  estas 
curvas  elementales,  para  reproducir  una  curva  determinada. 

Prescindiendo  por  ahora  de  la  utilidad  practica  que  esta  descomposicion  de  las 
curvas  pueda  tener  en  meteorologia,  vamos  a  ver  como  la  formula  que  estudiamos 
no  es  otra  cosa  que  la  expresion  analitica  de  las  curvas  periodicas  en  terminos  de  sus 
componentes  armonicas. 

Dicese  curva  armonica  la  que  representa  graficamente  el  movimiento  armonico  de 
un  punto. 

1  Uno  de  los  mejores  trabajos  practices  que  conocemos  sobre  el  empleo  de  las  series  de  Fourier  en 
fisica  matematica  es  el  de  W.  E.  Byerly  "An  Elementary  Treatise  on  Fourier's  Series  and  Spherical 
Harmonics"  publicado  por  Ginn  and  Co.,  New  York: 

La  obra  original  del  mismo  Fourier  "Theorie  Analytique  de  la  Chaleur"  es  quiza  de  tanta 
actualidad  como  cualquiera  otra  publicacion  moderna  sobre  esta  materia.  Puede  consultarse  tambien 
Mellor  "Higher  Mathematics  for  Students  of  Chemistry  and  Physics,"  Longmans,  Green  and  Co., 
39  Paternoster  Row,  London. 
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Movimiento  armonico  es  aquel  a  que  estaria  sometida 
la  proyeccion  sobre  un  diametro  de  un  punto  que  recorriera 
la  circunferencia  con  movimiento  uniforme.  Esta  circun- 
ferencia  recibe  el  nombre  de  circulo  de  referenda. 

Si  el  punto  M  (fig.  1)  recorre  la  circunferencia  con  mo- 
vimiento uniforme,  el  movimiento 
de  su  proyeccion  P  sobre  el  diame- 
tro BD  es  armonico;  mientras  M 
describe  una  vez  la  circunferencia 
en  el  sentido  ABCDA,  al  punto 


Fig  1 


P  recorre  el  diametro  en  el  sentido 
OBODO. 

En  gracia  de  aquellos  que  no 
esten  familiarizados  con  el  estudio  de  esta  clase  de  movimientos, 
permitasenos  recordar  algunas  definiciones,  que  hemos  de 
necesitar. 

Periodo  es  el  tiempo  que  emplea  el  punto  M  en  recorrer  la 
circunferencia,  o  el  punto  P  en  su  doble  recorrido  del  diametro. 

Elongation  es  la  distancia  variable  del  centre  0  al  punto  P. 

Amplitud  es  la  elongacion  maxima,  o  distancia  maxima 
a  que  puede  separarse  de  0  el  punto  P  en  su  movimiento ;  como 
se  ve,  es  el  radio  del  circulo  de  referencia. 

Los  puntos  O  y  A  son  el  origen  de  los  movimientos  de 
los  puntos  P  y  M  respectivamente.  El  origen  de  los  tiem- 
pos  es  arbitrario;  puede  tomarse  en  cualquier  punto  de  la 
circunferencia. 

Fase  es  la  distancia  angular  del  punto  movil  al  origen 
del  movimiento.  Tratandose  del  punto  que  se  mueve  sobre 
el  diametro,  hay  que  referirle  al  circulo  de  referencia  para 
determinar  sus  fases ;  asi  en  la  figura  1,  si  el  punto  se  encuentra 
en  P,  su  fase  esta  dada  por  los  grades  del  arco  MA  o  del  angulo 
MOA. 

Epoca  de  una  fase  determinada  es  el  instante  en  que  ocurre 
esa  fase;  asi,  en  nuestro  caso,  si  el  arco  MA  es  de  60°  y  ban 
transcurrido  dos  horas  desde  que  el  punto  se  encontraba  en 
el  origen  de  los  tiempos  diremos  que  las  2h  es  la  epoca  de  la 
fase  de  60°. 

Fase  inicial  es  la  distancia  angular  del  origen  de  los  tiempos 
al  origen  de  los  movimientos.  En  nuestra  figura,  si  tomamos 
al  punto  N  como  origen  de  los  tiempos,  la  fase  inicial  del 
movimiento  esta  dada  por  el  arco  NA. 

Supongamos,  para  fijar  las  ideas,  que  el  punto  M  emplea 
un  dia  en  recorrer  la  circunferencia — es  decir,  que  el  periodo 
de  su  movimiento  es  veinticuatro  horas.  A  partir  de  N,  origen 
de  los  tiempos  (fig.  2),  dividamos  la  circunferencia  en  24  partes 
iguales ;  y  rectifiquemosla  sobre  la  recta  NX;  esta  queda  dividida 
en  24  partes  iguales,  cada  una  de  las  cuales  representa  la 
longitud  de  un  arco  de  15°,  o  en  tiempo,  una  hora.  Hagamos 
que  el  circulo  de  referencia  se  desplace  sobre  la  recta  NX 
paralelamente  al  diametro  vertical,  con  un  movimiento  igual 
al  del  punto  que  la  recorre.  Al  fin  de  cada  hora  encontraremos 
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Fig.  2. 
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el  centro  O  en  los  puntos  sucesivos  en  que  ha  quedado  dividida  la  recta,  y  el  punto  que 
se  mueve  sobre  el  diametro  ira  describiendo  la  sinusoide  de  la  figura. 

Si  quisieramos  trazar  otra  curva  armonica  de  periodo  ^  del  de  la  anterior,  las  partes 
en  que  se  habria  de  dividir  el  circulo  de  referencia,  serian  de  doble  numero  de  grades  que 
las  anteriores,  pues  es  claro  que  para  que  la  circunferencia  quede  recorrida  en  doce  horas, 

el  punto  habra  de  recorrer  30°  cada  hora ;  y  asi  para  curvas  de  periodo  £,  £,  etc., del 

de  la  primera,  los  respectivos  circulos  de  referencia  habrian  de  dividirse  en  partes,  cuyo 
numero  de  grados  f uera  triple,  cuadruple,  etc. ;  por  donde  vemos  que  los  periodos  de  estas 
curvas  estan  en  razon  inversa  del  numero  de  grados  de  las  partes  en  que  se  dividen  sus 
circulos  de  referencia. 

Una  ordenada  cualquiera  en  la  curva,  representa  la  distancia  a  que  se  encuentra  de 
0  el  punto  P  en  la  hora  correspondiente.  Designemos  por  ylf  ?//,  i//',  etc.,  estas  ordenadas. 
Del  triangulo  rectangulo  POM  (fig.  1)  deducimos  como  expresion  general  de  las  ordena- 
das— 

OP  =  y  =  OM  sen  a  =  u  sen  ( U-\-x) 

donde  u  es  el  radio  del  circulo  de  referencia,  o  amplitud  de  la  curva;  U  la  fase  inicial,  y 
x  representa  en  cada  caso  el  mimero  de  partes  de  la  circunferencia  que  para  ese  momento 
ha  recorrido  el  punto  M. 

Partiendo  ahora  de  estas  consideraciones,  obtenemos  para  las  distintas  ordenadas  de 
las  curvas  armonicas  caracteristicas  de  una  curva  periodica  cualquiera  las  expresiones  si- 
guientes : 

Primera  Segunda  Tercera 

ordenada.  ordenada.  ordenada. 

Primera  curva  ut  sen  Ui  Ui  sen  (Ui  +  x)         ut  sen  (Ui-\-2x)         

Segunda  curva  u2  sen  U»  w2  sen  (C72+2x)        %  sen  (t/2+4o;)         

Tercera  curva  u3  sen  Us  u3  sen  (U3+3x)         u3  sen  (U3+6x) 


Sumando  algebraicamente  las  ordenadas  del  mismo  nombre,  obtenemos  para  las  orde- 
nadas correspondientes  de  la  curva  periodica  resultante  las  expresiones  siguientes: 

Resultante. 

Primera  ordenada  =  Ui  sen  Ui  +  Ui  sen  E72  +  u3  sen  U3  -f 

Segunda  ordenada  =  Ui  sen  (Ui+x)  +  tfcsen  (J72  +  2cc)  +  u3  sen  (U3+3x)  + 

Tercera  ordenada  =  u\  sen  (Ui-\-2x)  +  u2  sen  (t/z+4o;)  +  u3  sen  (U3+6x)  + 


que,  como  facilmente  se  ve,  no  son  mas  que  expresiones  particulares  de  la  formula  gene- 
ral— 

Ui  sen  (t/i  +  tf)   -f  Ui  sen  (C72  +  2x)    +  u3  sen  (U3  +  Bx)  + _ 

en  la  que  el  elemento  variable  x  va  tomando  los  valores  que  le  corresponden  segun  lo  expli- 
cado. 

Deducida  por  este  procedimiento  la  formula  de  Bessel,  aparece  claramente  lo  que  al 
principio  dijimos,  que  es  la  expresion  analitica  de  una  curva  periodica  en  terminos  de  sus 
componentes  armonicas ;  cada  uno  de  sus  terminos  representa  una  de  estas  componentes. 

Las  constantes  u  y  U  suelen  denominarse :  la  primera  coeficiente  armonic  lineal;  la 
segunda  coeficiente  armonico  angular  (o  angulo  de  posicion  del  punto  de  partida  del 
elemento  variable),  pues  ellas  son  las  que  caracterizan  el  movimiento  armonico  que  se 
estudia ;  la  u  da  la  amplitud,  y  por  tanto  la  importancia  del  movimiento ;  pues,  como  facil- 
mente se  ve  en  las  figuras  precedentes,  a  medida  que  decrece  u,  el  movimiento  del  punto 
P  sobre  el  diametro  va  siendo  mas  lento;  la  U  (fig.  2)  fija  la  epoca  de  las  distintas  fases; 
un  aumento  en  su  valor  hace  adelantar  la  epoca  de  todas  las  fases,  y  al  contrario  una  dis- 
minucion  la  retrasa. 

Es  de  suma  importancia  tener  presente  esta  propiedad  de  la  fase  inicial ;  si  en  una 
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estacion  cuya  longitud  fuera,  por  ejemplo,  15°  W.  de  Greenwich  se  quisiera  calcular  en 
tiempo  medio  de  Greenwich  en  vez  de  en  tiempo  medio  local,  se  adelanta  una  hora  el  origen 
de  los  tiempos,  y  por  consiguiente  hay  que  disminuir  el  valor  de  las  fases  iniciales,  15° 
en  I/,,  30°  en  U2,  45°  en  U3,  etc.  Por  esta  razon  debe  hacerse  constar  en  las  publicaciones 
el  tiempo  en  que  se  calcula. 

La  cantidad  M  que  aparece  en  la  formula  representa  la  media  aritmetica  de  las 
observaciones  que  entran  en  la  discusion.  Lo  que  se  desarrolla  en  la  serie  trigonometrica 
de  la  formula  es  la  diferencia  entre  cada  observacion  y  la  media  aritmetica  de  todas  ellas ; 
por  consiguiente,  para  que  nuestras  observaciones  (o  las  ordenadas  de  las  curvas  de 
observacion)  queden  representadas  segun  esta  formula,  hemos  de  agregar  la  media  arit- 
metica al  desarrollo  de  estas  diferencias. 

Creemos  que  lo  dicho  hasta  aqui  sera  suficiente  para  formarse  una  idea  general  del 
mecanismo  de  la  formula ;  en  la  aplicacion  que  haremos  de  ella  al  estudio  de  la  variacion 
diurna  de  la  presion  atmosferica  en  Manila  procuraremos  concretar  mas  estas  ideas. 

APLICACION  DE  LA  FORMULA  AL  ESTUDIO  ANALfTICO  DE  LA  VARIACI6N  DIURNA 
DE  LA  PRESI6N  ATMOSFERICA  EN  MANILA  EN  EL  PERlODO  1890-1909. 

En  la  preparation  de  los  dates  que  habian  de  entrar  directamente  en  el  calculo  hemos 
procedido  de  la  manera  siguiente: 

De  las  publicaciones  oficiales  de  este  Observatorio  tomamos  las  medidas  horarias  de 
cada  mes,  calculando  ademas  por  separado  la  media  correspondiente  a  la  hora  cero,  o 
primera  medianoche  del  dia  representative  de  cada  mes;  esta  media  es  el  valor  medio 
de  las  lecturas  de  la  primera  medianoche  de  cada  dia  del  mes.  Como  la  hora  cero,  o 
primera  medianoche  de  un  dia,  es  la  ultima  medianoche  del  anterior,  en  un  mes  cual- 
quiera  (Enero  por  ejemplo)  las  lecturas  de  las  Oh  seran:  Para  el  dia  1.°  la  lectura  de 
las  12  p.  m.  del  ultimo  dia  de  Diciembre  del  ano  anterior;  la  del  dia  2,  la  de  las  12  p. 
m.  del  dia  1.°  del  mismo  Enero,  etc.,  de  modo  que,  los  valores  que  entran  en  el  promedio 
de  esta  hora  son  los  mismos  que  entran  en  el  de  las  12  p.  m.  a  exception  de  la  lectura 
de  las  12  p.  m.  del  ultimo  dia  del  mes;  que  es  sustituida  por  la  correspondiente  del  mes 
anterior. 

Facilmente  se  ve  que,  teniendo  ya  calculada  la  media  de  las  12  p.  m.  basta,  para 
calcular  la  de  las  Oh,  restar  la  lectura  de  las  12  p.  m.  del  ultimo  dia  del  mes  de  que  se 
trata,  de  la  lectura  correspondiente  en  el  mes  anterior;  dividir  esta  diferencia  por  el 
numero  de  dias  del  mes  en  cuestion,  y  sumar  algebraicamente  el  resultado  a  la  media  de 
las  12  p.  m. 

Obtenidas  estas  25  medias  horarias,  se  han  calculado  las  medias  horarias  generales 
de  los  doce  meses  para  todo  el  periodo,  agrupando  los  meses  del  mismo  nombre,  y  pro- 
mediando  las  medias  obtenidas  para  cada  hora  de  un  mismo  mes  en  los  distintos  anos. 
Estas  medias  horarias  generales  se  han  aproximado  hasta  las  milesimas,  y  son  las  que 
forman  la  Tabla  I. 

La  Tabla  II  es  la  misma  Tabla  I  en  la  que  se  ha  eliminado  el  elemento  aperiodico. 
Vamos  a  explicar  brevemente  en  que  consiste  esta  correction,  y  cual  es  su  fin. 

Variacion  no  ciclica,  o  aperiodica. — En  la  explication  de  la  formula  vimos  que  las 
curvas  armonicas  elementales  constan  en  su  integridad  de  25  ordenadas,  correspondientes 
a  las  distintas  horas  del  dia,  de  Oh  a  24h,  ambes  inclusive.  Las  ordenadas  extremas  son 
iguales  para  una  misma  curva;  por  tanto  las  ordenadas  correspondientes  a  las  Oh  y  24h 
en  la  curva  resultante  han  de  ser  iguales,  por  ser  la  suma  algebraica  de  las  ordenadas 
correspondientes  en  las  distintas  curvas  elementales. 

Ahora  se  entendera  por  que  es  necesario  calcular  la  media  de  las  Oh ;  es  precise  ase- 
gurarse  de  que  la  curva  que  se  va  a  analizar  como  periodica,  cumple  con  el  caracter  esen- 
cial  de  periodicidad  de  que  sus  ordenadas  extremas  scan  iguales.  En  caso  de  no  serlo, 
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debemos  deducir  que  se  ha  introducido  en  la  marcha  del  fenomeno  algun  elemento  que 
perturba  su  periodicidad.  Esta  perturbacion  es  la  que  suele  denominarse  variacion  no 
ciclica  o  aperwdica. 

Cual  sea  de  ordinario  el  origen  de  esta  variacion,  se  entendera  facilmente  con  un 
ejemplo  bien  familiar. 

En  las  distintas  estaciones  del  ano  varia  la  temperatura  de  modo  que  va  subiendo 
en  primavera  y  verano,  y  bajando  durante  el  otono  e  invierno;  asi  que  en  Mayo,  por 
ejemplo,  aunque  no  puede  decirse  en  general  que  la  temperatura  de  un  dia  sea  mas 
elevada  que  la  del  precedente,  pero  si  queremos  calcular  la  variacion  media  diurna  de 
la  temperatura  en  ese  mes  para  un  ano  cualquiera,  tomando  para  cada  hora  del  dia  el 
valor  medio  de  aquella  hora,  derivado  de  los  treinta  y  un  dias,  notaremos  que  la  tempe- 
ratura de  la  primera  medianoche  del  dia  representative  de  Mayo  es  algo  mas  baja  que 
la  de  la  ultima  medianoche  del  mismo  dia. 

En  el  caso  de  observaciones  barometricas  nos  encontramos  con  un  fenomeno  ana- 
logo.  Es  bien  sabido  que  la  presion  atmosf  erica  esta  sujeta  a  una  oscilacion  anual, 
que  se  sobrepone,  como  en  el  caso  de  la  temperatura,  a  la  oscilacion  diurna.  Al  estu- 
diar  pues  los  caracteres  propios  de  la  variacion  diurna  mensual  de  la  presion  atmosferica, 
debemos  comenzar  por  la  eliminacion  del  elemento  aperiodico. 

Esta  eliminacion  suele  hacerse  aplicando  a  cada  hora  la  correccion  — 

D(l2-h) 
24 

en  la  que  D  se  obtiene  restando  del  valor  de  la  hora  24  el  de  la  hora  cero,  y  h  representa 
la  hora  en  la  que  se  va  a  aplicar  la  correccion.  No  estara  demas  alguna  discusion  sobre 
la  legitimidad  de  esta  formula.  Es  logica  esta  correccion  unicamente  si  se  admite  la 
hipotesis  de  que  la  diferencia  entre  los  valoies  de  las  dos  horas  extremas  se  ha  intro- 
ducido por  una  variacion  uniforme  durante  las  veinticuatro  horas. 

El  probar  si  es,  o  no,  legitima  esta  suposicion  no  carece  de  dificultad;  sin  embargo, 
cuando  la  investigacion  se  limita  a  las  veinticuatro  horas,  se  la  da  por  buena,  pues  parece 
la  unica  hipotesis  aceptable. 

Veamos  como  la  formula  supone  implicitamente  esta  hipotesis.  Si  la  variacion  se 
supone  uniforme  durante  las  veinticuatro  horas,  ha  de  estar  representada  por  una  recta. 
La  formula  de  correccion  sera  por  lo  tanto  la  ecuacion  de  la  recta  — 

y  =  b+aX  (1) 

Tomemos  como  origen  de  coordenadas  rectangulares  la  hora  cero,  y  llevemos,  como 
de  ordinario,  las  horas  al  eje  X.  En  la  formula  (1)  la  x  representara  las  horas  suce- 
sivas;  la  y,  la  correccion  que  hay  que  aplicar  en  cada  una.  Aplicada  la  correccion  a  las 
dos  horas  extremas,  sus  valores  han  de  venir  a  coincidir  ;  luego  hay  que  anadir  a  la  una 
y  restar  de  la  otra  la  semidiferencia  de  sus  valoies.  La  correccion  para  estas  dos  horas, 
llamando  D  su  diferencia  sera  por  consiguiente  — 

para  la  hora  24  ______  ±0  =6-f  24<* 

c 

para  la  hora  0  _______  +    =b 

& 

Eliminando  a  b  entre  estas  dos  ecuaciones  obtenemos  — 


Llevando  a  la  ecuacion   (1)   los  valores  obtenidos  para  a  y  b,  queda  — 

-D,D      D(l2-x) 

1J  =.-{  ---  -+-    -    /y-  —  -  S  --  — 

2~24X~        24 
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Los  signos  quedaran  determinados  por  el  que  tenga  D  en  cada  caso  particular. 
Vemos  pues  que,  para  la  deduccion  de  la  formula  de  correccion,  es  necesaria  y  suficiente 
la  hipotesis  de  que  la  variacion  aperiodica  se  verifique  de  un  modo  uniforme  durante 
las  veinticuatro  horas. 

Un  ejemplo  hara  ver  que  influencia  ejerce  esta  correccion  en  las  curvas. 

Tomemos  el  mes  de  Noviembre,  que  es  en  nuestro  caso,  el  que  tiene  mayor  co- 
rreccion. La  curva  llena  (fig.  3)  representa  la  variacion  sin  corregir;  la  punteada 
representa  la  misma  variacion  corregida. 

Por  la  forma  y  posicion  de  estas  curvas  se  ve  claramente,  que  la  eliminacion  del 
elemento  aperiodico  se  reduce  a  hacer  girar  la  curva  un  angulo  pequenisimo,  tomando 
las  12  a.  m.  como  centre,  sin  que  en  nada  varien  sus  ondulaciones  caracteristicas. 

Hecha  esta  correccion,  se  ban  calculado  para  la  Tabla  II  las  medias  diurnas,  divi- 
diendo  por  24  la  suma  de  las  24  medias  horarias  de  cada  mes,  de  Oh  a  23h ;  y  restando 
luego  en  cada  mes  la  media  diurna,  de  cada  media  horaria,  se  ha  formado  la  Tabla 
III,  que  es  la  que  contiene  los  datos  de  que  nos  vamos  a  servir  en  el  calculo. 

En  esta  ultima  evaluacion  hemos  aproximado  unicamente  hasta  las  centesimas; 
pero,  como  en  los  calculos  anteriores  llevamos  la  aproximacion  hasta  las  milesimas, 
la  inexactitud  de  los  datos  de  la  ultima  tabla  no  puede  pasar  de  media  centesima,  y 
la  suma  de  los  valores  positives  no  se  diferencia  en  valor  absolute  de  la  de  los  nega- 
tives mas  que  en  dos  o  tres  centesimas  a  lo  sumo,  resultado  que  no  se  hubiera  obtenido 
si  en  los  calculos  anteriores  se  hubiera  aproximado  unicamente  a  centesimas. 


Figs 

En  el  estudio  analitico  que  con  los  datos  de  esta  tabla  vamos  a  emprender  hemos 
empleado  el  tiempo  medio  local. 

Sabemos  que  lo  unico  que  necesitamos  para  la  determinacion  de  cada  una  de 
las  ondas  elementales,  en  las  que  vamos  a  resolver  la  variacion  barometrica  diurna  son  : 
Los  coeficientes  armonicos  lineales  u,  y  los  angulares  U  —  es  decir  las  constantes  de  la 
Formula  II,  por  ser  ellas,  como  al  principio  indicamos,  las  que  caracterizan  estas  ondas. 

Estas  constantes  no  pueden  hallarse  directamente;  hay  que  determinarlas  valien- 
dose  de  las  constantes  p  y  q  de  la  Formula  I,  que  son  las  que  se  calculan  directamente. 

Las  ecuaciones  que  relacionan  estas  constantes  son,  como  vimos  antes  : 

p  —  itsen  U\/-,\ 
q  —  ucos  U> 
De  aqui  se  obtienen  las  relaciones: 


s«L=tan  U=P.  sen  U=P  cos  u=g 
cos  U  q  u'  u 

De  la  primera  de  estas  tres  ecuaciones  podemos  deducir  la  constante  angular  U  en 
funcion  de  p  y  q;  y  de  las  dos  ultimas  en  funcion  de  p  o  q  y  u;  esta  ultima  se  obtiene 
por  la  ecuacion: 
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que  se  deduce  facilmente  de  las  ecuaciones  (1)  elevandolas  al  cuadrado,  y  sumando 
ordenadamente.  Toda  la  dificultad  esta  en  determinar  las  constantes  p  y  q. 

El  metodo  que  dio  Bessel  para  determinarlas,  y  que  explicaremos  al  fin  de  nuestro 
trabajo,  es  largo  y  penoso,  y  sobre  todo,  requiere  una  tension  constante  en  el  trabajo 
para  evitar  los  errores  en  que  facilmente  se  puede  incurrir  a  causa,  principalmente, 
del  considerable  numero  de  angulos  que  hay  que  reducir  y  por  la  f  recuencia  con  que  cam- 
bian  de  signo  las  lineas  trigonometricas  al  pasar  de  uno  a  otro  cuadrante.  De  este 
metodo  dice  al  P.  Dechevrens,  S.  J.  :  "La  methode  naturelle  n'est  que  pour  les  habiles." 
Este  Padre  ha  publicado  dos  folletos,  el  uno  titulado  "Methode  simplifiee  dite  des 
facteurs  pour  le  calcul  des  series  de  Fourier  et  de  Bessel  appliquees  a  la  Meteorologie," 
en  el  que  por  medio  de  tablas,  se  facilita  notablemente  el  calculo.  El  otro  folleto  lleva 
el  titulo  "Note  complementaire  a  la  methode  simplifiee  du  calcul,  etc."  l  donde  explica 
el  procedimiento  que  ha  seguido  para  la  formacion  de  sus  tablas. 

Las  tres  primeras  de  estas  tablas  ^irven  para  la  determinacion  de  los  coeficientes 
p  y  q,  en  series  de  4,  6,  8,  10,  12,  24  y  36  observaciones. 

Creemos  no  carecera  de  interes,  para  quienes  no  esten  avezados  a  esta  clase  de 
analisis,  resolver  practicamente  un  ejemplo  por  el  metodo  de  Bessel,  y  ver  en  que  esta 
la  utilidad  de  las  tablas  del  P.  Dechevrens  ;  por  esto,  aun  a  riesgo  de  parecer  nimios  en 
detallar  el  calculo,  vamos  a  seguir  paso  por  paso  el  metodo  de  Bessel  en  la  determinacion 
de  las  constantes  p  y  q  para  Enero. 

Las  ecuaciones  que  en  nuestro  caso  dan  los  valores  de  estas  constantes  para  los 
ocho  primeros  terminos  de  la  formula  (I)  son:  - 


=2  am  cos  ra#  q=^  am  sen  mx 

=2  am  cos  m(2x)  12q2=2  am  sen  m(2x) 

=5  am  cos  m(3x)  \2q.=$  am  sen  m(3x) 

12p4=2  am  cos  m(4x)  12^=2  am  sen  m(4x) 


El  signo  2  indica  que  la  cantidad  a  que  afecta  se  ha  de  desarrollar  en  una  serie 
de  terminos  de  su  misma  forma,  que  se  suceden  segun  la  serie  de  los  numeros  natu- 
rales  0,  1,  2,  etc.,  asi  que  escribiremos  : 


2  am  cos  mx=an  cos  0-j-c^  cos  x-\-a2  cos  2x-\-a3  cos  3x-{-  ................ 

x  representa  el  valor  de  una  de  las  partes  en  que  esta  dividido  el  circulo  de  referencia 
de  la  primera  curva  elemental.  Como  el  periodo  de  esta  curva  es  de  veinticuatro  horas, 
su  circulo  de  referencia  esta  dividido  en  24  partes,  y  por  tanto  #=15°.  Los  periodos 
de  las  otras  curvas  son  ^,  ^,  £  del  de  la  primera;  por  consiguiente  las  divisiones  de  sus 
respectivos  circulos  de  referencia  valen 

2x;  3x;  4x;  o  30°;  45°;  60°; 

am  representa  los  valores  de  la  Tabla  III  correspondientes  a  las  distintas  horas;  asi  en 
Enero,  por  ejemplo  — 

a0=46,  a1==17,  a2=—  20,  etc. 

Segun  esto,  para  la  determinacion  de  p^  en  Enero,  tendremos  — 

12p1==46  cos  0+17  cos  15°—  20  cos  30°—  45  cos  45°  ................ 

Solo  falta  hallar  en  las  tablas  trigonometricas  naturales  el  valor  de  estos  cosenos,  y 
ejecutar  las  operaciones  indicadas.  El  trabajo,  que  seria  improbo  si  se  hubiera  de  hacer 
la  multiplicacion  en  cada  termino,  puede  simplificarse  mucho  del  modo  siguiente: 

1  Tipografia  della  Pace  di  Filipo  Cuggiani,  Roma,  Via  della  Pace  num.  35. 

2  Al  fin  de  este  trabajo  daremos  la  prueba  matematica  de  estas  ecuaciones. 
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El  valor  de  los  cosenos  en  las  distintas  horas  es : 


0"  

cos  0°  — 

1 

.966 
.866 
.707 
.5 
.259 
0 
-   .259 
-   .5 
-   .707 
-   .866 
-   .966 

12" 

—      1 

1  

cos  15°  — 

13 

—        966 

2  .... 

cos  30°  — 

14 

—        866 

3  .... 

cos  45°  — 

15 

—        707 

4  .... 

cos  60°  — 

16 

5  .... 

cos  75°  — 

17 

—        259 

6  .... 

cos  90°  — 

18 

—     0 

7  .... 

....  cos  105°  —  —cos  75°  — 

19 

—        259 

8  .... 

....  cos  120°  —  —cos  60°  — 

20 

—        5 

9  

21 

—        707 

10  



22 

—        866 

11  . 

— 

23  .. 

..  -        .966 

Cada  uno  de  estos  numeros  aparece  como  factor  comun  de  varias  horas. 
mos  por  tanto  : 

12  P   =  (0- 


Tendre- 


-10-14  +  22).866 
+  (3-9-15+21)  .707+  (4-8-16+20)  .5+  (5-7-17+19)  .259 

Sustituyendo  en  vez  de  las  horas  que  aparecen  en  los  parentesis,  los  valores  correspon- 
dientes  tornados  de  la  Tabla  III  tenemos  por  fin: 

12  Pl  =(46-15)  1+  (17-79+63+59)  .966+  ......... 

Para  la  determinacion  de  ql  obtenemos  por  el  mismo  procedimiento  : 


0"  

sen  0° 

—  0 

12" 

—  o 

1  

sen  15° 

—  .259 

13 

—    259 

2  

sen  30° 

-  .5 

14 

—  —  5 

3  

sen  45° 

—  .707 

15 

—  -  .707 

4  

sen  60° 

—  .866 

16 

—  -  866 

5  

sen  75° 

—  .966 

17 

—  -  966 

6  

sen  90° 

—  1 

18 

—  —1 

7  

sen  105°  —  sen  75° 

—  .966 

19 

—  -  966 

8  

—  .866 

20 

—  -  866 

9  

—  .707 

21 

—  -  .707 

10  

—  .5 

22 

—  —  5 

11  . 

=  .259 

23  .. 

.259 

Teniendo  en  -cuenta  los  f actores  comunes  de  varias  horas,  escribiremos : 

12  q,  =  (1+11-13-23)  .259+  (2  +  10-14-22)  .5+  (3+9-15-21)  .707 

+  (4+8-16-20).866+(5+7-17-19).966+ (6-18)1 

y  sustituyendo  en  vez  de  las  horas  sus  valores  queda  determinada  q^. 

Para  la  determinacion  de  p2  y  q.2  procedemos  del  mismo  modo,  teniendo  en  cuenta 
que  las  divisiones  del  circulo  de  referencia  valen  30° ;  diremos  pues : 


0" cosO' 

1  cos  30C 

2  cos  60C 

3  cos  90C 

4  cos  120°  =  — cos  60C 

5  cos  150°  =  — cos  30C 

6  

7  

8  

9 

10  

11  . 


1 

12"  

^ 

.866 

13  

-   .866 

.5 

14  

-   .5 

0 

15 

.  -  0 

.5 

16  

—  —  .5 

.866 

17  

--  .866 

1 

18 

—  _  i 

.8C6 

19 

.  .  ..  --  .866 

.5 

20     

--  .5 

0 

21   . 

=  0 

.5 

22   

=   .5 

.866 

23  .. 

..  =   .866 
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12  p2  =(0-6  +  12-18)1+  (1-5-7  +  114-13-17-19+23).  866 

+  (2-4-8+10+14-16-20+22).5 


Oh sen  Oc 

1  sen  30' 

2 sen  60' 

3  sen  90C 

4  sen  120°=sen  60' 

5  

6  

7  

8  

9 

10  

11  .. 


=     0 


=     1 


12". 
.5          13 
.866     14  . 

15 

.866      16 
.5          17 


=     0 


0 

18  

-   .5 

19  

-   .866 

20  

—  1 

21 

—   .866 

22 

-   .5 

23  .. 

=       .5 
=       .866 
=     1 
=       .866 
=       .5 
=     0 
=  -   .5 
=  -   .866 

=:  —  1 

=  —   .866 
=  -   .5 


12  g,  =  (1  +  5-7-11+13  +  17— 19-23)  .5+  (2+4-8-10+14 

+16-20-22). 866+  (3-9  +  15-21)1 

El  circulo  de  referenda  de  la  tercera  curva  esta  dividido  de  45°  en  45° ;  por  consi- 
guiente : 


Ou cos  0°=     1         |  12" 

1  cos  45°=       .707     13 

2  ..  ..  cos  90°=     0  14 


3 

4 
5 
6 
7 
8 
9 

10 
11 


cos  135°  =  — cos  45°  =  -   .707 

-t 


15 
16 


=  -   .707 

=     0 
=       .707 

—     1 


=  -   .707  I  17  ..  =       .707 


=     0 


18 


=     0 


=       .707     19  ..  =-   .707 


-< 


20 


=        .707     21  ..  ..  =-   .707 


=      0 


22 


=     0 


=-   .707     23  =       .707 

12  p,=  (0-4+8-12  +  16-20)1+  (l-3-5+7+9-ll-13  +  15  +  17-19-21  +  23).707 


o" 


sen 


.707 


13 


14 


1  sen  45°  = 

2  sen  90°  = 

3  sen  135°=sen45°=       .707  |  15 

4  =     0  16 

5  =-   .707     17 

6  =-1          I  18 

7  

8  

9 

10  

11  . 


'=     0  12" =     0 

=  -   .707 

—  _i 

=  -   .707 

=     0 

=       .707 

=     1 

—  _   .707     19  =       .707 

=     0  20  =     0 

=       .707     21 =-   .707 

85      1  22    =-1 

=       .707     23  ..  ..  =  -   .707 


12  qz=  (1+3-5-7+9  +  11-13-15+17+19-21-23)  .707+  (2-6+10-14+18-22)  1 
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Las  divisiones  del  circulo  de  referenda  de  la  cuarta  curva  valen  60° ;  tendremos 


pues: 


ou cos  oc 

1  cos  60C 

2  cos  120°  = -cos  60C 

3  

4  

5  

6  

7  

8  

9  

10 

11  . 


-  .5 
.5 

1 
.5 

-  .5 

-  .5 
.5 


12" 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


12   p4=  (0-3  +  6-9+12-15+18-21)1+ (1-2-4+5+7-8-10 

+  11+13-14-16+17+19-20-22+23)  .5 


oh  

sen  0°—  0 

12h 

—  0 

1  

sen  60°—   .866 

13 

—    866 

2  

sen  120°—  sen  60°—   .866 

14 

—    866 

3  

—  0 

15 

—  0 

4  

—  -  .866 

16 

—  —  866 

5  

—  -  .866 

17 

—  —  866 

6  

—  0 

18 

—  0 

7  

—   .866 

19 

—    866 

8  

—   .866 

20 

—   .866 

9  

—  0 

21 

—  o 

10  

—  —  .866 

22 

—  —  866 

11  . 

--  .866 

23  .. 

..  --  .866 

I2q4=  (1+2-4-5+7+8-10-11  +  13  +  14— 16-17+19+20-22-23).  866 

Hemos  descendido  a  tantos  pormenores  porque  pretendemos  hacer  ver  que  las  tablas 
del  P.  Dechevrens,  para  la  determination  de  estos  coeficientes,  son  de  entera  confianza, 
por  estar  fundadas  en  el  metodo  mismo  de  Bessel,  y  tienen  la  ventaja  de  dar  verificadas 
ya  las  operaciones  mas  expuestas  a  errores.  Las  formulas  que  da  el  P.  Dechevrens  son : 

p1=  (0-12)  .083+  (1+23-11-13)  .08+  (2+22-10-14)  .072+  (3+21-9-15)  .059 

+  (4+20-8-16)  .042+ (5+19-7-17). 022 

q,=  (6-18)  .083+  (1+11-13-23)  .022+  (2+10-14-22)  .042+  (3  +  9-15-21)  .059 

+  (4+8-16-20)  .072+ (5+7-17-19).08 

P2=  (0+12-6-18).  083+  (1  +  11  +  13+23-5-7-17-19).  072+  (2  +  10+14 
n  +22-4-8-16-20)  .042 

qa=  (3+15-9-21)  .083+  (1  +  5+13  +  17-7-11-19-23)  .042+  (2+4 

+  14  +  16-8-10-20-22)  .072 
p3=  (0+8+16-4-12-20)  .083+  (1+7+9  +  15+17 

+23-3-5-11-13-19-21)  .059 
qa=  (2+10+18-6-14-22)  .083+  (1+3+9+11  +  17 

+  19-5-7-13-15-21-23)  .059 
[04=  (0+6+12+18-3-9-15-21).  083+  (1+5+7+11  +  13  +  17+19+23-2 
IV  -4-8-10-14-16-20-22)  .042 

q4=  (1+2+7+8+13  +  14+19+20-4-5-10-11-16-17-22-23).  072 
Comparandolas  con  las  halladas  por  el  metodo  de  Bessel,  vemos  que  son  identicas; 
los  coeficientes  de  los  parentesis  son  los  mismos  de  Bessel  divididos  por  12. 
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Por  medio  de  estas  tablas  determinamos  las  constantes  que  forman  la  Tabla  IV. 

Con  estas  constantes  se  pasa  a  determinar  los  coeficientes  armonicos,  por  medio  de 
las  formulas  que,  segun  antes  vimos,  los  relacionan;  damos  a  continuacion  las  que  he- 
mos  empleado,  y  el  procedimiento  que  hemos  seguido  en  el  calculo. 

=  tg  U  .'.  log  p—  log  q  =  log  tg  U 


=  logp—  log  sin  U 
o 

En   el  ejemplo  que  estamos   resolviendo  para  Enero  tendremos,  segun  esto: 

logp,  ==1.  127429 
log  q,  =  1.  738463 


log  pt—  log  tfx  ==  1.  888966  =  log  tg  U,  .'.  U,  =  13°  45'  37" 

logp,  ==  1.  127429 
log  sin  U,  --  =  f  .  376320 


logp,  -  log  sin  Ul  =  1.  751109  =  log  M,  .'.  u,=  0.5638 

Hay  que  tener  presente  en  que  cuadrante  se  halla  U.  Puede  seguirse  esta  regla, 
para  averiguarlo  facilmente. 

Si  p  y  q  son  positivas,  U  pertenece  al  primer  cuadrante;  si  p  es  positiva,  y  q  negati- 
va,  al  segundo  ;  si  ambas  son  negativas,  al  tercero  ;  finalmente,  si  p  es  negativa,  y  q  posi- 
tiva, al  cuarto. 

Asi,  en  este  mismo  mes,  al  determinar  U2  se  hallara  directamente  en  las  tablas 
trigonometricas  el  angulo  15°  23';  pero  como  p2  es  positiva,  y  q2  negativa,  hay  que  re- 
ducirle  al  segundo  cuadrante,  lo  que  da 

U2  =  164°  37' 

De  este  modo  se  ban  ido  deduciendo  para  cada  mes  los  valores  de  u  y  U,  que  damos 
en  la  Tabla  V. 

Aunque  en  rigor  basta  este  calculo  para  discutir  las  curvas  elementales,  nos  ha  pare- 
cido  calcular  por  separado  para  cada  mes  estas  cuatro  curvas. 

Los  valores  de  las  ordenadas  correspondientes  a  cada  hora,  juntamente  con  los  de  la 
resultante,  aparecen  en  la  Tabla  VI,  en  la  que  hemos  anadido  la  ultima  columna  con  las 
diferencias  entre  la  variacion  resultante  calculada  y  la  observada.  Lo  insignificante  de 
estas  diferencias  nos  parece  una  buena  comprobacion  de  los  calculos. 

Una  vez  calculadas  estas  ondas,  juzgamos  no  careceria  de  interes  anadir  las  12 
laminas,  en  cada  una  de  las  cuales  aparecen  las  curvas  armonicas  elementales  de  un 
mes  con  su  resultante.  Asi  se  podra  apreciar,  mas  facilmente  que  por  los  valores  de 
las  tablas,  la  importancia  de  cada  una  de  las  curvas.  En  ellas  se  ve  claramente  por  que 
basta  calcular  los  cuatro  primeros  terminos  de  la  formula;  la  serie  es  rapidamente  con- 
vergente,  de  modo  que  la  cuarta  sinusoide  se  confunde  casi  con  una  recta. 

Aunque,  para  la  verification  de  estos  ultimos  calculos,  da  tambien  tablas  de  P. 
Dechevrens,  juzgamos  que  seria,  si  no  mas  facil,  mucho  mas  breve  calcular  por  el  metodo 
natural. 

Por  via  de  ejemplo,  vamos  a  calcular  las  cuatro  componentes  de  Enero.  La  expre- 
sion  general  de  las  ordenadas  de  la  componente  de  veinticuatro  horas  es 


u  sen 

por  consiguiente,  en  nuestro  caso,  tendremos  como  expresion  general 

0.564  sen  (13°  46'+z) 

Para  determinar  en  las  distintas  horas  el  agulo  expresado  por  el  parentesis,  basta  re- 
cordar  que,  en  esta  primera  onda,  el  punto  que  recorre  el   circulo  de  referencia,   ha 

127182  -  5 
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7". 

.sen  118°  46'  = 

sen  61°  14'  = 

.877 

8 

.sen  133°  46'  = 

sen  46°  14'  = 

722 

9 

.  sen  148°  46'  = 

sen  31°  14'  = 

.519 

10  . 

.  sen  163°  46'  = 

sen  16°  14'  = 

.280 

11  . 

.  sen  178°  46'  = 

sen  1°  14'  = 

.022 

12  . 

.  sen  193°  46'  = 

—  sen  13°  46'  =  — 

.238 

descrito  un  arco  de  15°  al  fin  de  la  primera  hora,  30°  al  fin  de  la  segunda;  45°  al  fin 
de  la  tercera,  etc.,  de  modo  que  se  obtendra  para  las  distintas  horas : 

0" sen  13°  46'=   .238 

1  sen  28°  46'=   .481 

2  sen  43°  46'=   .692 

3  :.  sen  58°  46'=   .855 

4  sen  73°  46'=   .960 

5  sen  88°  46'  =  1 

6  sen  103°  46'  =  sen  76°  14'=   .971 

De.las  12h  en  adelante  se  van  repitiendo  los  valores  en  el  mismo  orden,  pero  con 
signo  contrario.     Basta  ahora  multiplicar  cada  una  de  las  funciones  naturales,  que  apa- 
recen  en  ultimo  miembro,  por  el  coeficiente  0.564,  para  obtener  los  valores  que  repre- 
sentan  las  distintas  ordenadas  de  esta  onda. 
Segunda  onda. — La  expresion  general  es: 

uz  sen    (U2-\-2x) 
que  en  nuestro  caso  se  reduce  a 

0.978  sen   (164°  37'  +  2x) 
El  elemento  variable  2x  aumenta  en  30°  cada  hora;  por  tanto 


0" sen  164C 

1  sen  194C 

2  sen  224C 

3  ..  ..  sen  254° 


37'=  sen  15C 
37'= -sen  14C 
37'= -sen  44 c 
37'= -sen  74° 


23'=  .265 
37'= -.252 
37'= -.702 
37'= -.964 


4"  ....  sen  284C 

5  sen  314C 

6  sen  344C 


37'=— sen  75C 
37'= -sen  45C 
37'= -sen  15C 


23' =-.968 
23'= -.712 
23'= -.265 


Desde  las  6h  comienzan  a  repetirse  los  valores,  pero  con  signo  contrario ;  de  modo  que, 
para  obtener  las  ordenadas  de  esta  curva,  basta  multiplicar  estos  valores  por  0.978,  e 
ir  repitiendo  esta  serie  de  productos  hasta  la  hora  23,  teniendo  cuidado  de  cambiar  de 
signo  en  cada  repeticion. 

Las  dos  ondas  restantes  se  obtienen  por  analogo  procedimiento  teniendo  en  cuenta 
el  valor  del  elemento  variable  en  cada  una. 


Tercera  onda. 

Cuarta 

onda. 

Oh 

sen  11 

0  1'  = 

.191 

0". 

.  sen  130° 

17'=- 

sen  49° 

43'=- 

.763 

1 

sen  56 

0  1'  = 

.829 

1  . 

.  sen  190° 

17'=- 

sen  10° 

17'=- 

.179 

2  .... 

..  sen  101 

0  1'  = 

sen  78° 

59'  = 

.982 

2  . 

.  sen  250° 

17'=- 

sen  70° 

17'=- 

941 

3  .... 

..  sen  146 

°  1'  = 

sen  33° 

59'  = 

.559 

3  . 

.  sen  310° 

17'=- 

sen  49° 

43'=- 

763 

4  .... 

..  sen  191 

0  1'  = 

—  sen  11° 

1'=- 

.191 

Este  es  el  modo  natural  de  proceder  en  el  calculo;  en  realidad  es  mucho  mas  sen- 
cillo,  si  nos  fijamos  en  que  los  angulos  van  siendo  complementarios  de  los  precedentes, 
desde  las  6h  en  la  primera  onda ;  desde  las  3h  en  la  segunda,  y  desde  las  2h  en  la  tercera. 
De  modo  que  basta  determinar  los  angulos  hasta  las  5'1  inclusive  en  la  primera  ohda ; 
hasta  las  2h  en  la  segunda,  y  hasta  la  lh  en  la  tercera,  y  sin  cuidarse  de  los  restantes, 
tomar  a  la  vez,  en  las  tablas  de  las  funciones  trigonometricas  naturales,  el  valor 
del  seno  de  cada  uno  de  estos  angulos,  y  el  de  su  complementario,  que  se  encuentra 
al  extreme  de  la  misma  linea.  Por  ejemplo,  para  la  hora  cero,  en  la  primera  onda  de 
Enero,  hallamos  el  valor  del  seno  de  13°  46',  que  es  0.238,  y  al  otro  extreme  de  la 
misma  linea,  encontramos  el  valor  0.971  correspondiente  al  seno  de  su  complementario 
76°  14',  que  escribiremos  al  frente  de  la  hora  6. 

Por  medio  de  las  tablas  del  P.  Dechevrens,  para  hallar  la  ordenada  correspondiente 
a  las  Oh  de  Enero,  por  ejemplo,  hay  que  aplicar  la  siguiente  formula: 

(0-12)  .083+  (1-13)  .08+  (2-14)  .072+  (3-15)  .059+  (4-16)  .042+  (5-17)  .022 

+  (7-19)  X  -.022+  (8-20)  X  -.042+  (9-21) 
X-.059+(10-22)X-.072+ (11-23)  X-.08 
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sustituir  en  vez  de  las  horas  que  aparecen  en  los  parentesis,  sus  valores,  y  ejecutar  las 
operaciones  indicadas.  En  nuestro  caso,  hechas  ya  las  restas  indicadas  en  los  parentesis, 
tendriamos 

.3lX.083+.8X.08+1.05x.072+l.lX.059-fl.04x042 

+.89X.022— .84X.022— .87X.042-.79X.059-.58X.072-.2X.08=0.1341 

For  el  metodo  natural,  una  vez  hallados  los  valores  de  las  lineas  trigonometricas, 
obtenemos  el  mismo  resultado  con  una  sola  multiplicacion : 

0.238X0.564=0.1342 

Claro  esta  que  el  metodo  del  P.  Dechevrens  es  menos  expuesto  a  errores,  por 
reducirse  a  sencillas  operaciones  aritmeticas ;  pero  creemos  que  basta  alguna  practica 
en  el  calculo  trigonometrico,  para  ganar  tiempo  por  el  metodo  natural. 

Epocas  de  las  fases  criticas  (mdximos  y  minimos). — La  determinacion  de  las 
epocas  de  las  fases  criticas  no  ofrece  dificultad.  Teniendo  presente  que  una  ordenada 
cualquiera  tiene  por  expresion 

u  sen  (U-\-x) 

la  condicion  necesaria  y  suficiente,  para  que  este  valor  sea  maximo,  es  que  el  seno  que 
aparece  en  la  expresion  general,  valga  la  unidad,  lo  que  deja  finalmente  como  condicion 
para  el  maximo. 

U+x=90°  6  450° 

segun  que  U  sea  menor  o  mayor  que  90°.     Esto  supuesto,  recordemos  que  el  elemento 

variable  x  tiene  por  valores  para  las  horas  sucesivas  0,  1,  2,  3, h,  en  la  primera 

curva;  0°,   15°,   2x15°,   3X15°,   &Xl5°,    del    mismo    modo    en    la    segunda, 

tercera  y  cuarta  tiene  respectivamente  por  valor  para  una  hora  cualquiera  h,  las  ex- 
presiones  /iX30°,  ftx45°,  /&X60°.     Segun  esto,  vamos  a  determinar  las  epocas  del  pri- 
mer maximo  para  las  cuatro  curvas  de  Enero. 
Para  la  primera  tenemos  como  condicion : 


15 


15 


6 


para  la  segunda: 


164 


=9h  31r 


para  la  tercera: 

11°   r-ffex450=90°.'./i=lh  46m 
cuarta : 

130°  17'-f-fex60°=4500.'.to=5h  20m 

Una  vez  determinada  la  epoca  del  primer  maximo,  quedan  determinadas  las  restantes, 
por  distar  entre  si  doce  horas  en  la  primera  curva,  seis  en  la  segunda,  cuatro  en  la 
tercera,  y  tres  en  la  cuarta.  Las  epocas  del  primer  maximo,  para  las  cuatro  curvas 
en  los  distintos  meses,  aparecen  en  la  siguiente  tabla: 

fipoca  del  primer  maximo. 


Meses. 

Primera 
onda. 

Segrunda 
onda. 

Tercera 
onda. 

Cuarta 
onda. 

Meses. 

Primera 
onda. 

Segunda 
onda. 

Tercera 
onda. 

Cuarta 
onda. 

Enero 

H.  m. 

5        5 

H.   m. 
9      31 

H.  m. 
1      46 

H.   m. 

5      20 

Julio  .  

H.   m. 
2      57 

H.  m. 
9      52 

H.    m. 
6      27 

H.  m. 
5    32 

5      47 

9      41 

1      44 

5      24 

3       9 

9      51 

6      52 

4    49 

Marzo 

5      23 

9      45 

1      10 

5      37 

Septiembre  

3      25 

9      41 

1      23 

4    38 

Abril 

5      3 

9      44 

0      49 

5      38 

Octubre  

4      16 

9      30 

1      11 

4    42 

Mayo 

4      10 

9      44 

7      15 

5      59 

Noviembre  

4      20 

9      25 

1      19 

4    42 

3      11 

9      49 

6      20 

0      0 

Diciembre      

4      26 

9      26 

1      40 

4    46 
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DISCUSI6N. 

En  lo  que  precede,  hemos  resuelto  la  oscilacion  diurna  de  la  presion  en  Manila  en 
una  serie  de  oscilaciones  armonicas  simples,  las  dos  primeras  de  las  cuales,  de  periodps  de 
veinticuatro  y  doce  horas,  son  mas  importantes  sin  comparacion  que  las  de  periodos  de 
ocho  y  seis  horas;  de  modo  que  los  principales  caracteres  de  la  oscilacion  diurna  de  la 
presion,  estan  representados  casi  exclusivamente  por  las  dos  primeras. 

Pero  este  analisis  de  las  curvas  no  es  de  suyo  otra  cosa  que  un  artificio  de  calculo; 
toda  curva  periodica  puede  someterse  a  el.  Por  lo  tanto,  las  curvas  armonicas  componentes 
pudieran  tener  una  signification  puramente  matematica,  sin  corresponder  a  nada  real. 

Helmholtz  al  tratar  de  esta  clase  de  analisis  dice:  "No  es  mas  que  una  fiction  mate- 
matica, admirable  por  lo  que  facilita  el  calculo,  pero  no  guarda  relation  necesaria  con  algo 
real."  Sin  embargo,  hablando  en  general,  no  puede  negarse  la  utilidad  practica  de  este 
analisis. 

Un  fenomeno  periodico  no  puede  resultar  sino  de  causas  periodicas,  y  de  periodos 
submultiplos  del  periodo  del  fenomeno;  de  lo  contrario,  la  periodicidad  de  este  quedaria 
perturbada.  Por  medio  de  este  analisis,  se  descompone  la  curva  en  las  unicas  elementales, 
que  son  posibles,  de  periodos  submultiplos  del  suyo ;  por  tanto,  es  sumamente  probable  que 
las  causas  parciales  que  influyen  en  el  fenomeno  aparezcan,  o  totalmente,  o  en  sus  caracteres 
principales,  en  la  manera  de  variar  de  las  curvas  elementales.  Si  desearamos  averiguar  si 
dos  f  enomenos  estan,  o  no,  relacionados  entre  si,  el  estudio  de  las  armonicas  caracteristicas 
de  ambos,  y  especialmente  su  modo  de  variar  durante  el  ano,  es  el  unico  procedimiento 
por  donde  podemos  llegar  a  cerciorarnos  de  si  existe,  o  no,  esa  relacion ;  el  que  en  las  cur- 
vas que  los  representan  no  aparezcan  rastros  de  relacion  mutua,  no  nos  dice  mucho ;  pues 
la  influencia  que  el  uno  pudiera  tener  en  el  otro,  puede  facilmente  estar  contrarrestada  por 
otras  causas  que  influyan  a  la  vez. 

Aunque  de  la  comparacion  de  dos  f  enomenos  por  medio  de  este  analisis  no  sacaramos 
mas  que  la  persuasion  de  que  no  estan  relacionados  entre  si,  bien  se  ve  que  no  se  podria 
dar  por  perdido  nuestro  trabajo. 

Viniendo  al  caso  particular  del  barometro,  vamos  a  traducir  lo  que  el  Dr.  Hann  siente 
de  este  estudio:1 

Por  mi  parte,  estaba  persuadido  de  que  todos  los  esfuerzos  que  se  hicieran  para  explicar  la 
oscilacion  diurna  del  barometro  por  medio  de  las  variaciones  diarias  de  los  elementos  meteorologicos 
en  una  region  cualquiera,  como  lo  intentaron  Kreil,  Blanford,  Renou  y  otros,  no  podia  llevar  a  con- 
clusion alguna;  he  publicado  una  serie  de  escritos,  en  los  que  doy  una  description  precisa  del  fenomeno 
Begun  se  manifiesta  en  toda  la  tierra,  lo  mismo  al  nivel  del  mar  que  a  todas  las  alturas  para  las 
que  existen  observaciones,  y  me  he  esforzado  en  presentar  los  resultados  en  una  forma  tal,  que  pudiera 
servir  de  base  a  una  teoria  fisico-matematica.  Con  este  objeto,  he  representado  todos  los  resultados 
de  las  observaciones  por  funciones  periodicas,  y  he  calculado  las  amplitudes  y  epocas  de  las  fases  de 
las  ondas  individuales,  con  sus  periodos  propios,  las  cuales,  al  combinarse,  producen  el  resultado  com- 
plejo  de  la  curva  barometrica  diaria,  tal  como  se  nos  presenta  a  la  observation  directa  *  *  * 
Este  es  el  unico  medio  con  que  nos  ha  sido  posible  hacer  progresos  notables  en  la  comprension  y  expli- 
cation del  fenomeno. 

Lo  que  principalmente  nos  proponemos  en  este  trabajo  es  comparar  los  resultados 
obtenidos  para  Manila  con  los  que  se  ban  obtenido  para  otras  regiones;  seguiremos  pues 
en  la  discusion  una  forma  tal,  que  se  preste  a  esta  comparacion. 

Primera  onda. — Uno  de  los  hechos  mas  f  undamentales  que  se  establecen  en  meteoro- 
logia  es,  que  las  capas  inferiores  de  la  atmosfera  reciben  su  calor  de  la  superficie  terrestre 
por  conduction,  y  lo  transmiten  luego  por  convection  a  las  inmediatas.  Este  movimiento 
de  conveccion  es  principalmente  ascendente ;  el  aire  se  expansiona  lateralmente  en  la  alta 

1  "Further  contributions  to  the  foundation  of  a  theory  of  the  daily  barometric  oscillation."  Q.  Jour- 
nal of  the  Meteorological  Society.  London,  1899,  p.  40. 
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atmosfera,  y  determina  en  consecuencia  un  descenso  en  el  barometro.  Lo  contrario  sucede 
con  el  enf riamiento  de  la  superficie  terrestre ;  el  aire  en  contacto  con  ella  se  enf ria,  y  con- 
trae,  dando  lugar  a  corrientes  descendentes  y  ascenso  en  el  barometro.  For  lo  tanto, 
en  virtud  de  esta  variacion  de  temperatura,  han  de  producirse  en  la  presion  atmosferica 
un  maximo  y  un  minimo,  hacia  las  horas  del  minimo  y  del  maximo  de  temperatura. 

El  cambio  de  calor  entre  la  superficie  terrestre  y  las  capas  inferiores  de  la  atmosfera 
es  de  caracter  periodico,  siendo  su  principal  periodo  de  veinticuatro  horas.  Decimos  su 
principal  periodo,  porque  esta  variacion  no  es  simplemente  armonica;  de  modo  que  si  se 
analizara  por  separado,  no  quedaria  representada  por  un  sencillo  movimiento  armonico, 
sino  que,  ademas  de  una  componente  fundamental  de  periodo  de  veinticuatro  horas,  resul- 
tarian  otras  menos  importantes  de  periodos  de  doce,  ocho,  seis,  horas. 

La  primera  onda  en  el  analisis  de  la  curva  barometrica  es  de  periodo  de  veinticuatro 
horas,  y  presenta  invariablemente  el  maximo  y  el  minimo  hacia  las  horas  de  minima  y 
maxima  temperatura.  Este  caracter  general  de  la  curva  hace  ya  sospechar,  que  representa 
el  caracter  principal  de  la  variacion  de  presion,  causada  por  el  cambio  de  calor  entre  la 
superficie  terrestre  y  las  capas  inferiores  de  la  atmosfera. 

Tanto  la  amplitud  como  las  epocas  de  las  fases  criticas  de  esta  onda  presentan  en 
Manila  caracteres  especiales,  que  no  podemos  menos  de  hacer  notar.  La  amplitud  es 
menor,  y  la  epoca  del  maximo  ocurre  mas  temprano  que  en  la  mayor  parte  de  las  estaciones, 
en  las  que  se  ha  verificado  este  analisis.  Hann  hizo  notar  estas  mismas  propiedades  para 
Mauricio  (isla)  Singapore  y  Batavia,  en  su  publicacion  "Untersuchungen  iiber  den 
Taglichen  Gang  des  Barometers." 

Mr.  Eliot  incluye  estas  tres  estaciones  en  la  tabla  siguiente,  que  copiamos  del  t.X  de 
"Indian  Meteorological  Memoirs,"  p.  73,  con  varias  de  las  estudiadas  por  el  en  la  India, 
para  comparar  con  ellas  a  Trivandrum,  unica  estacion  de  la  India  que  presenta  estos 
caracteres.  A  estas  cuatro  estaciones  puede  anadirse  Manila  como  aparecera  por  los  va- 
lores  de  la  tabla. 


Estacion. 

Latitud 
norte. 

Valor 
medio 
anual  de 
la  ampli- 
tud Ml. 

Epoca  del 
maximo. 

Valor 
medio 
anual  de 
la  fase 
inicial 
Ui. 

Estacion. 

Latitud 
norte. 

Valor 
medio 
anual  de 
la  ampli- 
tud Ml. 

Epoca  del 
mdximo. 

Valor 
medio 
anual  de 
la  fase 
inicial 
Ui. 

Trivandrum    _              _ 

o    / 
8    31 

0.413 

H.    m. 
4    36 

O      ' 

21      2 

o     / 
25    26 

0  774 

H.    m. 

1    38 

O       1 

335    33 

Madras    

13      4 

.595 

6      1 

359    38 

21      9 

767 

7      4 

343    52 

Bombay  

18    54 

.472 

7    56 

330    53 

Lahore 

31    34 

667 

8    32 

322      2 

Rangoon  

16    46 

.685 

6    20 

355      3 

20      6S 

310 

2    54 

46    35 

Chittagong  

22    21 

.582 

6    50 

347    45 

1    17 

531 

4     18 

25    35 

Calcuta  (Alipore)  

22    32 

.667 

6    56 

346      5 

Batavia 

6    US 

626 

4    19 

25    17 

Kurrachee                       

24    47 

.497 

7      5 

343    49 

Manila 

14    35 

544 

4    15 

25    58 

Aden 

12    45 

.779 

7    26 

338    24 

Mr.  Eliot  termina  la  discusion  de  esta  onda  con  las  siguientes  palabras : 

La  primera  componente  de  la  oscilacion  diurna  en  Trivandrum  presenta,  por  tanto,  propiedades 
caracteristicas  que  difieren  grandemente  de  las  otras  estaciones  de  la  India,  pero  que  se  asemejan  de 
una  manera  notable  a  las  de  Batavia  y  Singapore,  y  tambien  (aunque  no  tanto)  a  Mauricio.  Los 
principales  caracteres  son:  su  pequena  amplitud,  y  la  temprana  epoca  del  maximo.  Esto  probable- 
mente  se  debe  en  gran  parte,  sino  del  todo,  a  su  posicion  maritima  y  su  baja  latitud. 

Esta  semejanza,  que  hace  notar  Mr.  Eliot  entre  Trivandrum,  Batavia  y  Singapore, 
aparece  mucho  mas  marcada  en  Manila,  tanto  en  la  amplitud,  como  en  la  fase  inicial  y 
epoca  del  maximo ;  es  pues  Manila  una  confirmation  mas  de  estas  propiedades  que  se  sos- 
pechan  para  estaciones  costeras  de  baja  latitud. 

Variacion  anual  de  u^ — La  variacion  anual  de  u1  y  la  variacion  anual  de  la  absorcion 
del  calor  en  las  capas  inferiores  de  la  atmosfera  guardan  entre  si  relation  de  causa  a 
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efecto,  como  lo  han  hecho  notar  Hann  ("Untersuchungen,  etc.")  y  Angot  ("Annales  du 
Bureau  central  Meteorologique  de  France,"  1887,  p.  B.  297),  y  puede  verse  confirmado 
en  los  trabajos  de  Mr.  Eliot  ("Indian  Meteorological  Memoirs,"  vol.  XII,  p.  347). 

Breves  consideraciones  haran  ver  que  se  verifica  lo  mismo  en  Manila.  Como  puede 
verse  en  la  Climatologia  publicada  por  el  P.  Coronas  en  el  t.  II  de  "El  Archipielago 
Filipino,"  los  meses  de  Febrero,  Marzo  y  Abril  son  los  mas  secos,  y  de  nubosidad 
minima;  Julio,  Agosto  y  Septiembre  por  el  contrario,  son  los  mas  lluviosos,  y  de  nubo- 
sidad maxima.  Las  circunstancias  pues,  son  sumamente  favorables  en  los  tres  primeros 
meses  para  la  absorcion  del  calor  en  las  capas  inferiores  de  la  atmosfera,  y  desfavorables 
en  los  tres  ultimos,  por  interceptar  la  gran  cantidad  de  nubes  buena  parte  de  la  energia 
calorifica,  y  por  haberse  de  convertir  parte  de  la  que  llega  a  la  superficie  terrestre  en 
calor  latente  de  evaporacion. 

Una  ojeada  a  los  valores  de  ulf  hace  ver  que  los  tres  meses  en  que  alcanza  sus  valores 
maximos  son  Febrero,  Marzo  y  Abril,  en  los  que  la  absorcion  de  la  energia  calorifica  en 
las  capas  inferiores  de  la  atmosfera  ha  de  ser  maxima;  y  llega  a  sus  valores  minimos 
en  Julio,  Agosto  y  Septiembre,  epoca  en  que  la  absorcion  es  minima. 

Mr.  Eliot  hace  notar  en  las  publicaciones  antes  citadas,  que  en  el  transito  de  la 
estacion  lluviosa  a  la  estacion  seca  suele  verificarse  un  maximo  relative  en  los  valores 
de  Wj.  Este  maximo  relative  aparece  efectivamente  en  Manila  en  Octubre,  mes  que  puede 
considerarse  como  divisorio  entre  las  dos  estaciones. 

Segunda  onda. — Dos  son  los  caracteres  que  principalmente  llaman  la  atencion  en 
esta  onda;  primero,  ser  notablemente  mas  importante  que  las  otras  tres  en  todos  los 
meses,  y  en  todas  las  regiones  para  las  que  se  ha  verificado  este  analisis;  segundo,  la 
constancia  en  la  epoca  de  sus  fases  criticas. 

Veanse  en  el  siguiente  cuadro  las  epocas  del  primer  maximo  halladas  por  Mr.  Eliot 
para  varias  estaciones  de  la  India: 


Latitud 
Norte. 

Estaci6n. 

Enero. 

Febre- 
ro. 

Marzo. 

Abril. 

Mayo. 

Junio. 

Julio. 

Agosto. 

Sep- 
tiem- 
bre. 

Octu- 
bre. 

No- 
viem- 
bre. 

Diciem- 
bre. 

o        / 
8      31 
10      50 
13        4 
18      54 
24      47 
34      10 
Valor  m 
14      35 

Trivandrum 

H.  m. 
9    30 
9    50 
9    45 
9    43 
9    41 
9    48 
9    42 
9    31 

H.   m. 
9    38 
9    59 
9    52 
9    51 
9    58 
10      3 
9    53 
9    41 

H.  m. 
9    34 
9    59 
9    50 
9    53 
9    58 
10      6 
9    50 
9    45 

H.  m. 

9    36 
9    58 
9    49 
9    54 
9    57 
9    58 
9    47 
9    44 

H.  m. 
9    37 
9    52 

9    48 
9    54 
9    58 
9    52 
9    47 
9    44 

H.   m. 
9    42 
9    52 
9    53 
10      2 
10    10 
9    50 
9    54 
9    49 

H.  m. 

9    47 
9    53 
9    57 
10      6 
10      9 
9    58 
9    58 
9    52 

H.   m. 
9    41 
9    53 
9    53 
9    57 
9    51 
10    22 
9    55 
9    51 

H.   m. 

9  29 
9  48 
9  39 
9  45 
9  49 
9  47 
9  41 
9  41 

H.   m. 

9  16 
9  38 
9  28 
9  25 
9  28 
9  49 
9  29 
9  30 

H.   m. 

9    14 
9    39 
9    27 
9    23 
9    29 
9    49 
9    28 
9    25 

H.   m. 

9    22 
9    47 
9    33 
9    36 
9    35 
9    49 
9    36 
9    26 

Trichinopoly 

Madras 

Bombay    _        

Leh 

edio  para  la  India  Tropical-  .- 
Manila                       .-  

Esta  onda  es,  por  su  doble  oscilacion  diurna,  la  que  produce  las  mareas  caracteris- 
ticas  del  barometro.  Es  por  tanto  de  sumo  interes  averiguar  su  origen. 

Las  discusiones  en  que  entran  Hann,  Angot  y  Lord  Kelvin  para  averiguar  el  origen 
de  esta  onda,  son  demasiado  largas  para  que  podamos  reproducirlas  aqui  en  todos  sus 
detalles.  Vamos  a  intentar  resumirlas  de  la  manera  mas  breve  y  clara  que  nos  sea 
posible. 

El  fenomeno,  que  este  representado  por  esta  onda,  ha  de  tener  estas  dos  propiedades : 
Primera,  influir  mas  que  ninguno  otro  en  el  barometro;  segunda,  ser  independiente  de 
circunstancias  locales. 

Desde  luego,  no  se  ve  que  otro  origen  pueda  tener  la  causa  de  este  fenomeno,  que 
la  influencia  del  sol  sobre  la  atmosfera  terrestre.  Pero  no  puede  deberse,  dice  Lord 
Kelvin,  a  la  atraccion  del  sol  que  origina  las  mareas,  pues  en  ese  caso  encontrariamos 
una  influencia  mucho  mayor  de  la  misma  especie,  debida  a  la  luna,  siendo  asi  que,  en 
realidad  la  marea  barometrica  lunar  es  nula,  o  apenas  perceptible.  Parece  pues  cierto 
que  ha  de  atribuirse  a  la  variacion  de  la  temperatura.  Ahora  bien ;  esta  variacion  en  las 
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capas  inferior es  de  la  atmosfera  esta  sujeta  a  todas  las  variaciones  posibles  a  que  dan 
lugar  las  circunstancias  locales  de  cada  estacion,  mientras  que  la  accion  periodica  de 
los  rayos  solares  que,  dia  tras  dia  se  verifica  en  las  capas  superiores,  ha  de  producir 
necesariamente  movimientos  periodicos  de  gran  regularidad. 

Esto  explica  facilmente  la  regularidad  de  esta  onda.  Lo  que  no  aparece  tan  claro 
es  el  porque  de  su  doble  oscilacion  diaria  y  su  notable  amplitud.  La  oscilacion  de  la 
temperatura  en  las  capas  superiores  de  la  atmosfera  ha  de  ser,  en  su  forma  general, 
analoga  a  la  oscilacion  de  la  temperatura  en  las  inferiores,  pues  las  causas  que  la  pro- 
ducen  son  las  mismas ;  la  insolacion  durante  el  dia,  y  la  radiacion  durante  la  noche.  La 
primera  hace  que  el  maximo  se  verifique  poco  despues  que  el  sol  adquiere  su  altura 
maxima ;  la  segunda  produce  el  minimo  poco  antes  de  la  salida  del  sol.  i  Como  pues  esta 
oscilacion  termica,'  de  un  solo  maximo  y  un  solo  minimo,  ha  de  poder  dar  origen  a  la 
doble  oscilacion  barometrica  que  estudiamos? 

Para  solucionar  la  dificultad,  Hann  acude  al  unico  medio,  que  como  indicamos  antes, 
existe  para  averiguar  la  relacion  que  guardan  entre  si  dos  fenomenos.  Resuelve  la 
variacion  diurna  de  la  temperatura  en  sus  componentes  armonicas,  y  halla  que,  de  las 
dos  ondas  principales  que  completan  el  fenomeno,  la  primera  (de  periodo  de  veinticuatro 
horas)  aparece  constantemente  con  una  amplitud  mayor  que  la  segunda  (de  periodo  de 
doce  horas)  inversamente  a  lo  que  sucede  con  la  presion.  He  aqui  algunos  ejemplos: 

Ui  U2 

Latitud    6°    0',  Pacifico  ecuatorial 0°.87  0°.ll 

Latitud  16°    0',  Pacifico  tropical 1°.39  0°.31 

Latitud  33°    0',  Paciflco  subtropical 0°.97  0°.25 

Latitud  48°  12',  Zona  templada  (Viena) 2°.86  0°.51 

A  pesar  de  esta  discordancia  en  las  componentes  armonicas  de  los  dos  fenomenos, 
Hann  no  duda  en  asignar  como  causa  de  la  doble  oscilacion  diaria  del  barometro  esta 
doble  oscilacion  diaria  de  la  temperatura. 

Al  aparecer  esta  afirmacion  del  eminente  meteorologo,  varies  cientificos  opusieron  a 
su  hipotesis  estas  dos  dificultades :  Primera,  la  variacion  termica  que  se  analiza  es  la  que 
se  observa  en  las  capas  inferiores ;  la  que  se  trata  de  estudiar  es  la  de  las  superiores, 
y  puede  dudarse  de  si  la  marcha  de  esta  segunda  es  precisamente  del  mismo  caracter 
que  la  primera.  Segunda,  la  doble  oscilacion  diaria  de  la  temperatura  es  unicamente  un 
resultado  del  calculo,  una  mera  ficcion  matematica,  que  no  puede  servir  de  base  a  la  ex- 
plicacion  de  un  fenomeno  real.  En  realidad  no  existe  variacion  diaria  de  temperatura 
con  dos  maximos  y  dos  minimos;  y  si,  a  pesar  de  esto,  obtenemos  una  doble  onda  para 
la  variacion  de  la  temperatura,  por  empefiarnos  en  representarla  por  una  serie  de  senos,  los 
resultados  que  obtenemos  son  matematicamente  obligados,  y  no  pueden  servir  para  ex- 
plicar  un  fenomeno  de  observacion ;  la  causa  de  este  debe  buscarse  en  algo  real  y  natural. 

En  cuanto  a  la  primera  dificultad,  queda  resuelta  con  lo  que  antes  indicamos,  acerca 
de  la  forma  de  la  variacion  de  la  temperatura  en  las  capas  superiores;  atendiendo  a  las 
causas  que  la  producen,  ha  de  ser  de  la  misma  forma  general  que  la  que  observamos  en 
las  inferiores. 

La  solucion  de  la  segunda  pide  alguna  mayor  explicacion.  Es  sabido  que,  si  una 
masa  limitada  de  fluido,  liquido  o  gaseoso,  entra  en  oscilaciones  simples  o  pendulares, 
las  amplitudes  de  estas  oscilaciones  obedecen  a  las  condiciones  propias  del  fluido  (sus 
dimensiones,  temperatura,  etc.),  de  modo  que,  si  la  masa  del  fluido  recibe  un  impulse 
unico,  se  movera  siempre  con  aquella  clase  de  ondulaciones,  que  le  son  posibles  en  virtud 
de  sus  condiciones  propias.  Estas  son  las  que  suelen  denominarse  oscilaciones  libres. 

Si  el  impulse  se  verifica  periodicamente,  se  fuerza  a  la  masa  a  entrar  en  oscilaciones 
del  mismo  periodo  que  el  del  impulse,  aunque  estas  nuevas  oscilaciones  no  coincidan  con 
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la  forma  de  oscilacion  a  que  pertenecen  las  ondas  libres.  Si  este  impulso  periodico  no 
fuera  de  caracter  simplemente  armonico — es  decir,  si  no  pudiera  quedar  representado 
por  una  simple  sinusoide,  podria  resolverse,  por  el  teorema  de  Fourier,  en  sus  compo- 
nentes  armonicas  de  periodos  1,  ^,  i,  etc.,  del  suyo,  de  modo  que  el  efecto  debido  a  este 
impulso  seria  el  mismo  que  producirian  obrando  a  la  vez  varies  impulses  simplemente 
armonicos,  cada  uno  de  los  cuales  estuviera  representado  por  una  de  las  componentes 
armonicas.  Notese  bien  que  esta  completamente  fuera  de  proposito  la  cuestion  de  si 
se  dan,  o  no,  en  la  practica  esos  impulses  componentes;  basta  que  el  impulso  real  pueda 
considerarse  como  resultante  de  ellos. 

Si  el  periodo  de  alguno  de  estos  impulses  componentes  coincide  con  el  de  las  oscila- 
ciones  libres  de  la  masa  del  fluido,  la  oscilacion  forzada  resultante  adquiere  una  ampli- 
tud  desproporcionalmente  grande;  de  modo  que,  aunque  la  amplitud  de  este  impulso  sea 
relativamente  pequefia  comparada  con  la  de  los  otros  impulses  componentes,  el  efecto 
por  el  producido  es  mucho  mayor;  por  donde  se  ve  claramente  que,  si  se  analizara  por 
medio  de  la  serie  de  senos  la  oscilacion  compleja  resultante,  la  curva  armonica  corres- 
pondiente  a  las  oscilaciones  libres  seria  de  gran  amplitud,  mientras  que  la  correspon- 
diente  en  el  analisis  del  impulso  seria  de  amplitud  pequena,  y  no  por  esto  dejaria  de  ser 
la  segunda,  causa  de  la  primera. 

Viniendo  ahora  a  nuestro  caso:  La  curva  que  representa  la  variacion  diurna  de  la 
temperatura  en  las  capas  superiores  de  la  atmosfera,  ha  de  ser  necesariamente  asime- 
trica,  por  serlo  sus  dos  causas,  la  insolation  y  radiacion,  que  dependen  de  la  duration 
del  dia  y  de  la  noche.  Si  pues,  al  analizarla,  vemos  que  queda  reproducida  casi  exclu- 
sivamente  por  la  combination  de  las  dos  componentes  armonicas  de  periodos  de  veinti- 
cuatro  y  doce  horas,  podemos  considerar  la  variacion  real  que  observamos,  como  la  re- 
sultante de  dos  variaciones  armonicas  de  la  temperatura;  una  de  periodo  de  veinti- 
cuatro  horas,  y  otra  de  doce;  no  importa  que  no  existan  de  hecho  dos  causas,  que  den 
origen  a  estas  dos  variaciones;  basta,  para  nuestro  proposito,  que  el  fenomeno  se  pro- 
duzca  del  mismo  modo  que  si  en  realidad  existieran. 

Ahora  bien;  el  Dr.  Margules,  siguiendo  las  indicaciones  de  Lord  Kelvin,  emprendio 
el  dificil  y  laborioso  trabajo  de  calcular  las  oscilaciones  que  en  la  atmosfera  terrestre 
ha  de  causar  su  calefaccion  y  enfriamiento  periodico.1  En  la  primera  de  sus  publica- 
ciones  probo  que  el  periodo  de  oscilacion  libre  de  nuestra  atmosfera  es  de  doce  horas 
proximamente. 

Estos  trabajos  y  conclusiones  del  Dr.  Margules  han  sido  acogidos  con  entusiasmo, 
y  han  dado  ocasion  a  tan  interesantes  publicaciones  como  la  de  Paul  Jaerisch  "Zur 
Theorie  der  hydrodynamischen  Gleichungen  in  spharischen  Koordinaten"  (Meteorologis- 
che  Zeitschrift,  1907,  pp.  481-498)  ;  la  de  W.  Trabert  "Die  Theorie  der  tag-lichen  Luft- 
druckschwankung  von  Margules,  und  die  tagliche  Oszillation  der  Luftmassen"  (Met.  Zeit., 
1903,  pp.  481-501)  ;  la  de  W.  J.  Humphreys  "On  the  diurnal  variations  of  atmospheric 
pressure"  (Bulletin  of  the  Mount  Weather  Observatory,  1912,  vol  5,  part  2,  pp.  132-156) 
y  la  mas  completa  quiza  de  todas  de  Horace  Lamb  "Atmospheric  Oscillations"  (Proc.  Roy. 
Soc.,  1911,  vol.  84,  pp.  551-572)  .2 


1  "tlber  die  Schwingugen  periodisch  erwarmter    Luft."     Sitzungsb.     Wien.,    Akad.    Marzo,     1890. 
"Luftbewegungen  in  einer  rotirenden   Spharoidschale."     Ibid.    Abril,    1890;    Enero   y    Diciembre,    1893. 
Puede  verse  una  traduccion   al  ingles  en  "The  Mechanics   of  the   Earth's   Atmosphere"   by   Cleveland 
Abbe.   Smithsonian  miscellaneous   collections — 843 — pag.  296. 

2  Estando  ya  en  la  imprenta  nuestro  trabajo  se  ha  recibido  en  este  Observatorio  una  publicacion 
de    Prof.    Frank    H.    Bigelow    "La    termodinamica  de  la  circulacion,  y   la   radiacion   de   la   atmosfera 
terrestre"    (Boletin  de  la  Oficina  Meteorologica  Argentina  No.  3).     En  la  pagina  15  propone  dos  di- 
ficultades   contra   la   teoria:    (1)    Las   ondas   semidiurnas   de  presion   desaparecen   dentro   de   los   2,000 
metros  de  la  superficie;    (2)    el  termino  de  calor    (Qi  — Qo)    no  estaba   incluido   en   la  ecuacion  general 
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Resumiendo  pues  brevemente,  podemos  explicar  la  doble  oscilacion  y  notable  am- 
plitud  de  la  segunda  onda  del  modo  siguiente:  Al  recibir  nuestra  atmosfera  un  impulso 
termico  (si  es  licita  esta  denomination)  se  produce  en  ella  el  mismo  efecto  que  si  re- 
cibiera  simultaneamente  dos,  que  pudieran  estar  representados  por  sinusoides  de  perio- 
dos  de  veinticuatro  y  doce  horas  respectivamente.  La  amplitud  o  intensidad  del  segun- 
do  es  menor  que  la  del  primero,  pero  la  oscilacion  que  en  la  atmosfera  causa,  ha  de  ser 
de  gran  amplitud,  por  coincidir  su  periodo  con  el  de  las  oscilaciones  libres,  propias  de 
nuestra  atmosfera;  mientras  que  las  oscilaciones  producidas  por  el  primero  quedan  no- 
tablemente  reducidas,  por  no  coincidir  su  periodo  con  el  de  las  oscilaciones  libres. 

La  hipotesis  de  que  la  segunda  onda  se  deba  a  la  absorcion  y  radiacion  de  calor  en 
la  alta  atmosfera,  esta  confirmada  por  la  comparacion  de  la  variation  anual  de  estas,  y 
la  de  la  amplitud  u.,.  Bastaran  a  nuestro  proposito  breves  indicaciones.1 

Uno  de  los  factores  mas  importantes  en  la  influencia  que  la  absorcion  y  radiacion 
pueden  tener  en  los  cambios  de  presion,  es  la  presencia  o  ausencia  de  nubes  en  la  atmos- 
fera. La  cantidad  de  calor,  que  absorbe  una  atmosfera  seca  y  despejada,  se  emplea  en 
dilatarla,  dando  lugar  a  movimlentos  ascendentes  y  laterales,  que  se  traslucen  en  el 
barometro  por  cambios  de  presion.  Al  contrario;  la  energia  calorifica  que  llega  a  una 
atmosfera  cargada  de  nubes,  se  convierte  en  gran  parte  en  calor  latente  de  evaporation,  lo 
que  no  puede  afectar  en  nada  a  la  presion.  Por  tanto,  en  la  estacion  seca,  los  cambios  de 
presion  debidos  a  estas  causas  seran  mayores  que  en  la  estacion  lluviosa,  por  la  notable 
diferencia  de  nubosidad  entre  estas  dos  estaciones. 

La  Tabla  V  hacer  ver  claramente  esta  diferencia  en  los  valores,  de  u2,  que  son  cons- 
tantemente  mayores  en  los  meses  de  la  estacion  seca  que  en  los  de  la  lluviosa,  y  en  los 
que  resalta  sobre  todo  la  separation  perfectamente  bien  definida  de  las  dos  estaciones 
en  el  salto  de  Mayo  a  Junio  que  en  estos  valores  se  verifica. 

Vamos  a  terminar  la  discusion  de  esta  onda  comparando  sus  principales  caracteres 
en  Manila  con  los  que  presenta  en  las  mismas  estaciones  que  comparabamos  en  la  primera 
onda. 

El  cuadro  que  sigue  esta  copiado  del  Tomo  X  de  "Indian  Meteorological  Memoirs," 

de  movimiento  empleada  por  Margules,  Jaerisch  y  Gold.  No  es  la  primera  vez  que  Mr.  Bigelow 
impugna  esta  teoria;  vease  por  ejemplo  su  publicacion  "Studies  on  the  diurnal  periods  in  the  lower 
strata  of  the  Atmosphere"  (Monthly  Weather  Review,  1905,  p.  93)  y  el  capitulo  9,  pagina  458  de  su 
"Report  on  the  internacional  cloud  observations,"  1898-99.  No  parece  se  pueda  suponer  que  cientificos 
como  Lord  Kelvin  y  Hann  ignoraran  que  la  oscilacion  semidiurna  va  disminuyendo  con  la  altura 
hasta  hacerse  nula  o  casi  imperceptible;  y  a  pesar  de  esto  no  dudaron  en  establecer  su  teoria. 
Deberian  pues  suponer  que  la  oscilacion  semidiurna  reaparece  en  la  alta  atmosfera.  Algo  de  esto 
parece  indicar  Mr.  Ernest  Gold  en  varies  pasajes  de  su  iriteresante  memoria,  "The  International 
Kite  and  Balloon  ascents"  que  obtuvo  el  primer  premio  ofrecido  en  1908  por  la  Sociedad  Meteorologica 
Alemana.  Para  abreviar  nos  limitaremos  a  traducir  dos  lineas  de  la  pagina  120  "De  aqui  que  cam- 
bios de  presion  en  la  superficie  terrestre  van  unidos  a  cambios  en  direccion  opuesta  en  la  estratoesfera, 
que  dan  lugar  a  corrientes  entre  la  estratoesfera  y  la  tropoesfera." 

La  segunda  dificultad  pudiera  inducir  a  creer  que  a  los  tres  eminentes  matematicos  que  nombra 
se  les  habia  pasado  por  alto  la  influencia  que  el  termino  de  calor  libre  pudiera  tener  en  la' ecuacion 
general  de  movimiento,  siendo  asi  que  el  Dr.  Margules  despues  de  introducir  en  la  ecuacion  diferencial 
de  movimiento  el  termino  de  calor  libre  termina  asi:  "Para  nuestro  proposito  sera  mas  conveniente 
considerar  las  variaciones  de  presion  como  consecuencia  de  las  variaciones  de  temperatura,  no  como 
consecuencia  del  flujo  variable  de  calor.  Por  tanto  volvemos  a  la  ecuacion  (4a) " — donde  ciertamente 
ya  no  entra  el  termino  de  que  habia  Mr.  Bigelow. 

1  Vease  una  larga  y  detallada  discusion  sobre  el  particular  en  el  Tomo  XII  de  "Indian  Meteorologi- 
cal Memoirs,  Calcutta,"  paginas  328-347. 

Pueden  consultarse  tambien  Langley:  "Researches  on  solar  heat  and  its  absorption  by  the  earth's 
atmosphere"  publicado  en  el  No.  15  de  "Professional  Papers  of  the  Signal  Service,"  y  en  el  Vol. 
I  de  "Annals  of  the  Astrophysical  Observatory  of  the  Smithsonian  Institution,"  Cleveland  Abbe 
"Preparatory  Studies  for  Deductive  Methods  in  Storm  and  Weather  Predictions." 
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pagina  74,  en  el  que  hemos  reducido  a  milimetros  los  valores  que  en  el  original  estan 
dados  en  pulgadas  para  facilitar  la  comparacion. 


Estaeion. 

Amplitud  media 

U2.  « 

Fase  media  [7z. 

•» 

TJ 

!i 
6s 

o  c 

s 

<W 

Relaci6n  del  maxi- 
mo  al  minimo 
de  W2. 

Mes  en  que 
ocurre  el 
mfiximo. 

Mes  en  que 
ocurre  el 
minimo. 

fipocas  extre- 
mas  del  maxi- 
mo  de  W2. 

Mes  de  la 
6poca  mas 
temp  ran  a. 

Mes  de  la 
epoca  mas 

tardia. 

Mas 
tem- 
prano. 

Mas 
tarde. 

Trivandrum  .-          -     

mm. 
1.100 
1.115 
.977 

1.023 

.946 

1.028 
.797 
.908 
.918 
1.010 

.628 
.733 
.995 
.959 

.919 

0           / 

164     14 
158      5 
156    51 

153      3 

148      5 

147    50 
155    40 
159    59 
150    13 
154    44 

138    18 
163    44 
156    18 
159    56 

160      5 

H.  m. 
9    32 

9    44 
9    46 

9    54 

10      4 

10      4 
9    49 
9    40 
9    59 
9    51 

10    23 
9    32 

9    47 
9    40 

9    40 

1.4 
1.2 

1.5 

1.4 

10 
.  0 

1.3 
1.6 
1.5 
1.2 
1.4 

1.3 
1.2 
1.3 
1.2 

1.4 

Febrero 

H.  m. 
9    14 
9    27 
9    15 

9    37 

9    39 

9    41 
9    28 
9    20 
9    37 
9    28 

9    56 
9    20 
9    31 
9    25 

9    25 

H.    m. 

9    47 
9    57 
9    58 

10    12 

10    19 

10    24 
10    10 
9    58 
10     17 
10    10 

10    36 

9    45 
9    58 
9    54 

9    52 

Noviembre  _ 
id 

Julio. 
Id. 
Id. 

JAgosto. 

|  Marzo. 
\Agosto. 
Julio. 
Junio. 
Abril. 
Agosto. 
Julio. 
jFebrero. 
\Julio. 
Febrero. 

Julio. 
Febrero. 

iJulio. 
JAgosto. 

Madras        

Marzo 

id 

Bombay  ..  

Febrero 

id 

id 

id 

{Noviembre  . 
Diciembre.- 

Noviembre  _ 
id 

Chittagong                       

Calcutta  

Marzo  

Julio 

Febrero  
Noviembre  - 

Junio  
id 

Aden  

tuore 

Allahabad 

id 

Noviembre  _ 
Octubre  

Magpur                                 -     

id 

Lahore  

Febrero 

id 

id 

Noviembre  . 

Singapore  

jMarzo  

I        id 

Batavia  

lAbril  
Octubre 

1 
id 

id 

Manila    

/Marzo  

1       id 

{Noviembre  . 

lAbril  

r 

Diciembre 

Por  este  cuadro  puede  apreciarse  de  un  modo  general  la  constancia  en  los  caracteres 
principales  de  esta  onda  en  las  distintas  regiones.  Resalta  sobre  todo  la  concordancia 
casi  perfecta  de  la  relacion  entre  el  maximo  y  minimo  de  la  amplitud,  y  la  epoca  del  ano 
en  que  se  verifica  el  minimo. 

Notanse  sin  embargo  algunas  discrepancias,  que  indican  no  estar  del  todo  libre  esta 
onda  de  la  influencia  de  circunstancias  locales.  Asi  por  ejemplo,  en  algunas  regiones, 
la  amplitud  presenta  dos  maximos  y  dos  minimos  en  el  ano,  mientras  que  en  otras,  Ma- 
nila entre  ellas,  desaparece  esta  doble  variacion.  En  estaciones  que,  como  Batavia,  Mau- 
ricio  y  Adden,  presentan  en  Octubre  o  Noviembre  el  maximo  absolute,  se  verifica  otro 
relative  en  Marzo  o  Abril. 

La  razon  de  estas  variaciones  es  bastante  obvia.  La  onda  segunda,  tal  como  resul- 
ta  del  analisis,  ha  de  ser  una  onda  compleja,  resultante  de  la  interferencia  de  dos,  con  el 
mismo  periodo  de  doce  horas.  Una  de  estas  (que  denominaremos  onda  semidiurna  prin- 
cipal) es  la  producida  por  las  variaciones  de  temperatura  en  las  capas  superiores  de  la 
atmosfera,  segun  lo  anteriormente  explicado;  la  otra  (onda  semidiurna  secundaria)  re- 
presenta  las  oscilaciones  que  en  la  presion  ha  de  producir  la  componente  semidiurna  de 
la  temperatura  en  las  capas  inferiores.  La  primera  de  estas  dos  presenta,  segun  Hann 
y  Angot,  una  doble  variacion  durante  el  ano,  siendo  maxima  en  los  equinoccios,  y  mini- 
ma en  los  solsticios.  La  segunda,  mas  irregular  y  de  menor  amplitud  que  la  primera, 
por  obedecer  a  influencias  locales,  ha  de  pasar  por  un  solo  maximo  (en  la  estacion  seca) 
y  un  solo  minimo  (en  la  lluviosa). 

Facilmente  se  ve  como  de  la  interferencia  de  estas  dos  ondas  puede  resultar  la 
desaparicion  completa,  o  casi  completa,  de  la  doble  variacion  anual  en  la  onda  que  estu- 
diamos;  pues  el  equinoccio  de  otono  tiene  lugar  en  la  mayor  parte  de  las  regiones  en  la 
estacion  lluviosa,  y  el  solsticio  de  invierno  hacia  el  comienzo  de  la  estacion  seca,  de  modo 
que  el  maximo  y  minimo  correspondiente  a  estas  epocas  en  la  onda  principal,  quedan  con- 
trarrestados,  o  en  todo  o  en  gran  parte,  por  coincidir  con  las  epocas  de  amplitud  mini- 
ma y  maxima  de  la  onda  secundaria.  Lo  contrario  sucede  con  el  solsticio  de  verano  y 
equinoccio  de  primavera,  por  coincidir  el  primero  con  la  estacion  lluviosa,  y  el  segundo 
con  la  seca. 

M.  Angot,  en  la  publicacion  ya  citada,  p.  B311  y  siguientes,  propone  un  metodo  bas- 
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tante  ingenioso  para  la  separacion  de  estas  dos  ondas.  No  juzgamos  oportuno  detener- 
nos  en  la  explication  de  todo  el  procedimiento,  que  nos  ocuparia  demasiado  espacio. 
Nos  limitaremos  a  dar  en  la  tabla  de  la  pag.  25,  los  resultados  que  por  este  medio  he- 
mos  obtenido  para  Manila.1 

No  podemos  comparar  nuestros  resultados  con  los  de  otras  regiones,  pues  no  sabe- 
mos  que  se  haya  verificado  este  analisis  mas  que  por  el  mismo  Angot,  y  este  unica- 
mente  ha  tornado  uno  o  dos  meses  para  las  nueve  regiones  que  somete  a  este  calculo; 
juzgamos  sin  embargo  util  haber  comprobado  para  Manila  este  caracter  que  se  atri- 
buye  a  la  onda  semidiurna  principal. 

Tercera  onda.  —  Dos  son  las  unicas  hipotesis  posibles  acerca  del  origen  de  esta  onda  : 

l.a  Que  provenga  de  una  accion  independiente  de  periodo  de  ocho  horas. 

2.a  Que  sea  parte  del  ef  ecto  producido  por  las  causas  que  originan  la  primera  onda  ; 
pues  como  indicabamos  en  la  discusion  de  esta,  las  causas  de  la  primera  onda  no  pue- 
den  representarse  por  un  movimiento  simplemente  armonico. 

Lo  primero  que  ocurre  contra  la  primera  hipotesis  es,  que  no  se  conoce  accion  fisica 
independiente  de  periodo  de  ocho  horas.  A  esto  se  anade  que  la  amplitud  de  la  tercera 
componente  es  muy  pequena,  lo  que  dice  mejor  con  el  caracter  de  un  .residue,  que  de  un 
efecto  total  debido  a  una  causa  independiente. 

En  confirmation  de  estas  razones  generates  viene  la  propiedad  que  en  todas  las  re- 
giones se  nota  en  esta  onda  ;  la  variacion  anual  de  su  amplitud,  y  de  la  epoca  de  sus  f  ases 
criticas  sigue  una  ley  pareclda  a  la  primera,  con  la  particularidad  de  que  el  sentido  de 
esta  variacion  es  casi  perfectamente  inverso,  lo  que  hace  sospechar  que  la  tercera  onda 
es  una  correction  de  los  valores  de  la  primera. 

Dada  la  pequena  influencia  que,  tanto  esta  onda  como  la  cuarta,  han  de  tener  en  la 
variacion  barometrica,  creemos  no  vale  la  pena  detenerse  mas  en  su  discusion.  No 
queremos  sin  embargo  omitir  una  particularidad  de  la  tercera  onda  que  hace  notar  el  Dr. 
Hann;  esto  es,  que  hacia  los  equinoccios,  la  epoca  de  las  fases  experimenta  una  inver- 
sion notable,  y  que  esta  inversion  concuerda  exactamente  con  la  epoca  de  los  minimos  de 
la  amplitud.  Esta  propiedad  aparece  bien  marcada  en  Manila;  la  inversion  de  que  ha- 
bla  el  Dr.  Hann  se  verifica  de  Abril  (OM9m)  a  Mayo  (7h  15m)  y  de  Agosto  (6h  52ra)  a 
Septiembre  (lh  23m)  ;  los  minimos  en  la  amplitud  se  verifican  tambien  de  Abril  a  Mayo, 
y  de  Agosto  a  Septiembre,  como  puede  verse  en  la  Tabla  V. 

1  La  formula  a  que  finalmente  se  llega  por  este  metodo  de  Angot,  es 
O.mm  926fef'f  cos4  X  cos  (2m+64°-2«) 


en  la  que  h  representa  la  presion  media;  8  la  declinacion  del  sol;  r  el  radio  de  la  orbita  terrestre;  X 
la  latitud  del  lugar;  m  las  distintas  horas  expresadas  en  grades,  de  modo  que  sera  cero  para  las  Oh;  15" 
para  la  lh;  30°  para  las  2h,  etc.;  e  la  ecuacion  del  tiempo.  Los  valores  que  hemos  tornado  para  estas 
cantidades  son: 

//  =  759.10  X  =  14°  35' 

9,  cos3  8 

r    ' 

Enero  ........................................    +4°  48'  ........................................................  0.9033 

Febrero    ....................................   +7°  01'  ........................................................  0.9722 

Marzo  ........................................    +4°  19'  ........................................................  1.0043 

Abril  ..........................................   +0°  02'  ........................................................  0.9597 

Mayo  ..........................................    -1°  46'  ........................................................  0.8734 

Junio  ..........................................   +0°  08'  ........................................................  0.8197 

Julio  ..........................................    +2°  41'  ........................................................  0.8424 

Agosto   ......................................   +1°  50'  ........................................................  0.9208 

Septiembre    ..............................    —2°  32'  ........................................................  0.9852 

Octubre    ....................................    -7°  00'  ........................................................  0.9806 

Noviembre  ................................    -7°  18'  ........................................................  0.9193 

Diciembre    .  ..  —  1°  56'  ..  ..  0.8746 
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CONCLUSI6N. 

De  todo  lo  hasta  aqui  expuesto  se  deduce,  que  la  teoria  moderna  de  la  doble  osci- 
lacion diaria  del  barometro  supone  ser  dos  las  causas  del  fenomeno.  La  una  esta 
completamente  libre  de  influencias  locales,  y  es  precisamente  la  que  introduce  en  la  curva 
barometrica  esa  doble  oscilacion,  siempre  tan  misteriosa  para  los  que  ban  tratado  de 
investigar  su  origen.  Esta  representada  por  la  que  hemos  denominado  onda  semidiurna 
principal,  y  es  debida  a  la  variacion  de  la  temperatura  en  las  capas  superiores  de  la 
atmosfera.  La  otra  varia  segun  las  circunstancias  locales  de  cada  region,  por  deberse 
a  la  variacion  de  la  temperatura  en  las  capas  inferiores.  La  curva  que  la  representa 
suele  denominarse  onda  termica.  En  la  tabla  de  la  pagina  26  damos  sus  valores  para 
todos  los  meses;  estan  deducidos  restando  de  la  variacion  observada  (Tabla  III)  la  onda 
semidiurna  principal. 

En  las  Laminas  XIII  y  XIV  pueden  verse  las  curvas  que  representan  la  onda 
termica,  y  la  onda  semidiurna  principal,  para  varies  meses. 

Estas  son  sin  duda  las  dos  ondas  de  que  habla  M.  Angot  en  su  "Traite  elementaire 
de  Meteorologie"  (p.  104)  al  explicar  la  variacion  diurna  del  barometro. 

Sin  saber  de  que  ondas  se  trata  no  es  facil  darse  perfecta  cuenta  de  esta  expli- 
cacion.  Hemos  juzgado  pues  sera  de  interes  traducirla  aqui,  pues  al  mismo  tiempo  que 
se  entendera  mejor  teniendo  a  la  vista  las  curvas,  servira  como  discusion  de  estas. 

"La  variacion  diurna  del  barometro  no  esta  aun  explicada  de  una  manera  completa. 
Las  investigaciones  mas  recientes  parecen  indicar  que  resulta  de  la  superposicion  de  dos 
oscilaciones  distintas ;  una,  que  presenta  dos  maximos  y  dos  minimos  por  dia,  cuyo  pe- 
riodo  es,  por  consiguiente,  de  medio  dia,  y  que  se  denomina,  por  esta  razon,  oscilacion 
u  onda  semidiurna;  otra,  que  presenta  un  solo  maximo  y  un  solo  minimo  en  veinticuatro 
horas,  y  que  se  denomina  onda  diurna. 

"Los  dos  maximos  diarios  de  la  onda  semidiurna  seran  iguales  entre  si,  lo  mismo 
que  los  dos  minimos.  En  toda  estacion  del  ano,  y  para  toda  la  tierra,  los  dos  maximos  se 
verificaran  hacia  las  10h  de  la  manana  y  de  la  tarde;  los  dos  minimos  hacia  las  4h  de 
la  manana  y  de  la  tarde.  La  diferencia  entre  los  maximos  y  los  minimos  de  la  onda 
semidiurna — es  decir,  la  amplitud  de  esta  onda,1  varia  poco,  para  una  misma  estacion, 
durante  el  ano;  es  maxima  en  los  equinoccios,  y  minima  en  los  solsticios.  Esta  ampli- 
tud no  dependera  de  condiciones  topograficas,  sino  unicamente  de  la  latitud — cerca  de 
2  mm.  en  el  ecuador,  disminuye  primero  lentamente  hasta  los  tropicos,  luego  muy  ra- 
pidamente,  de  modo  que  llega  a  ser  muy  pequena  en  latitudes  altas.  Es  probable  que 
esta  onda  semidiurna  este  causada  por  la  accion  del  calor  solar  en  toda  la  masa  de  la 
atmosfera;  pero  no  se  ha  podido  explicar  todavia  convenientemente  el  modo  como  se 
produce. 

"Mientras  que  la  onda  semidiurna  es  la  misma  para  todas  las  regiones  situadas 
en  la  misma  latitud,  la  onda  diurna  tiene  una  amplitud  muy  variable,  que  depende  de 
las  condiciones  topograficas  de  la  region,  y  que  puede  ser  muy  diferente  aun  para  esta- 
taciones  vecinas.  En  general  es  mucho  mas  debil  en  estaciones  costeras  que  en  estaciones 
continentales,  y  su  amplitud  esta  en  razon  directa  de  la  amplitud  de  la  variacion  diurna 
de  la  temperatura.  Las  horas  del  maximo  y  del  minimo  varian  notablemente  segun 
los  paises  o  las  estaciones  del  ano;  el  maximo  se  verifica  mas  temprano  en  verano  que 
en  invierno;  su  epoca  media  se  acerca  a  las  9h;  el  minimo,  al  contrario,  cuya  epoca 
media  es  cerca  de  las  16h,  se  produce  mas  temprano  en  invierno  que  en  verano.  El 
periodo  de  bajada  dura  pues,  por  termino  rnedio  solamente  siete  horas,  mientras  que 

1  M.  Angot  toma  aqui  como  amplitud  de  la  onda,  la  amplitud  de  una  oscilacion  completa;  otros 
autores,  y  el  mismo  Angot  en  la  publicacion  que  hemos  citado  varias  veces,  toma  como  amplitud  de  la 
onda  la  amplitud  de  una  semioscilacion,  o  la  ordenada  maxima  de  la  curva. 
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el  periodo  de  subida  se  extiende  a  diez  y  siete  horas;  la  subida  no  es  uniforme  durante 
este  tiempo;  rapida  inmediatamente  despues  del  minimo,  se  detiene  durante  la  noche,  y 
crece  de  nuevo  con  mas  rapidez  entre  la  salida  del  sol  y  la  hora  del  maximo. 

"Mientras  que  la  teoria  de  la  onda  semidiurna  no  esta  hallada  aun  1  la  onda  diurna 
se  explica,  por  el  contrario,  completamente  por  el  efecto  local  de  la  variacion  diurna 
de  la  temperatura.  Durante  la  noche,  a  medida  que  se  produce  el  enfriamiento  en  las 
capas  inferiores  de  la  atmosfera,  estas  se  contraen,  y  para  llenar  el  vacio  que  pro- 
duciria  esta  contraccion,  el  aire  afluye  por  todos  lados,  sobre  todo  por  la  parte  superior ; 
la  masa  de  aire  situada  sobre  la  region  que  se  considera  se  hace  por  lo  tanto  mayor, 
y  la  presion  aumenta.  Desde  la  salida  del  sol  las  capas  inferiores  se  calientan  rapi- 
damente  y  tienden  a  elevarse;  pero  la  calefaccion  se  produce  primero  en  las  capas 
que  estan  en  contacto  con  la  superficie  terrestre;  las  capas  de  aire,  comprendidas 
entre  la  superficie  y  las  capas  superiores  que  no  se  calientan  aun,  no  aumentan  al  principio 
de  volumen  de  una  manera  apreciable;  lo  mismo  que  sucede  cuando  se  calienta  el 
aire  en  una  camara  cerrada,  la  elevacion  de  temperatura  se  manifiesta  entonces  prin- 
cipalmente  por  un  aumento  de  presion;  el  primer  efecto  de  la  calefaccion  diurna  es  por 
consiguiente  un  aumento  de  presion  en  las  capas  que  estan  en  contacto  con  el  suelo ;  el 
barometro  comienza  pues  a  subir. 

"Un  poco  mas  tarde  entre  7h  y  llh  de  la  manana,  segun  los  paises  y  las  estaciones 
del  ano,  la  calefaccion  de  las  capas  inferiores  es  suficiente  para  determinar  por  fin  un 
movimiento  ascendente  general  del  aire  sobre  la  estacion;  el  aire  se  dilata;  pero  desde 
que  la  columna  de  aire  situada  sobre  la  region  caliente  alcanza  mayor  altura  que  la 
que  cubre  las  regiones  circunvecinas,  el  aire  se  extiende  hacia  los  lados,  desde  la  region 
mas  caliente  a  los  meridianos  mas  distantes  del  sol ;  la  cantidad  de  aire  que  insiste  sobre 
la  estacion  disminuye  por  lo  tanto,  y  el  barometro  baja.  Hacia  el  fin  del  dia,  cuando 
la  temperatura  comienza  a  disminuir  rapidamente  en  las  capas  inferiores,  la  corriente 
ascendente  pierde  actividad,  cesa  y  acaba  por  ser  reemplazada  durante  la  noche  por  un 
movimiento  inverse;  el  barometro  sube  entonces,  rapidamente  primero  mientras  la  tem- 
peratura disminuye  con  rapidez;  despues  mas  lentamente  hasta  la  salida  del  sol." 

Notese  que  la  explicacion  que  de  la  onda  termica  o  diurna  da  M.  Angot  coincide 
en  gran  parte  con  la  primera  teoria  de  la  variacion  diurna  del  barometro  propuesta  por 
Raimond,  y  perfeccionada  luego  por  Kreil  y  Renou,  y  que  era  insuficiente  cuando  se  la 
aplicaba  a  la  variacion  diurna  total  de  la  presion. 

Concluyamos  pues,  que  el  analisis  de  las  curvas  que  representan  los  diversos  ele- 
mentos  meteorologicos  de  cada  region,  aunque  un  tanto  largo  y  penoso,  puede  ser  de 
gran  utilidad,  tanto  para  el  estudio  simultaneo  de  varies  fenomenos,  cuyas  relaciones 
mutuas  queramos  investigar,  como  para  profundizar  mas  en  las  causas  que  en  cada  uno 
de  ellos  puedan  intervenir. 


1  No  sabemos  si  M.  Angot  tendria  noticia  de  las  investigaciones  del  Dr.  Margules,  pues  ni  aqui 
ni  en  el  trabajo  que  publico  en  los  Anales  del  Euro  Central  de  Francia,  hace  mencion  de  ellas, 
siendo  de  tan  capital  importancia  en  la  explicacion  de  esta  onda.  El  Dr.  Hann  en  publicaciones  pos- 
teriores  a  las  de  Angot,  se  lamenta  de  que  los  trabajos  del  Dr.  Margules  no  sean  mas  conocidos,  y  de 
que  no  hubieran  llamado  aun  la  atencion  de  la  "Meteorologische  Zeitschrift." 


DEDUCCION  DE  LAS  ECUACIONES  DE  BESSEL,  PARA  LA 

DETERMINACION  DE  LOS  COEFICIENTES 

DE  LA  FORMULA  (I). 


Supongamos  que  se  nos  da  un  numero  cualquiera  n  de  observaciones. 

Representemos  por  a,,,  alf  a,,  etc.,  las  observaciones  correspondientes  a  las  horas  0, 
1,  2,  etc. 

Vimos  en  la  pagina  43  que  estas  observaciones  estan  representadas  por  las  ecua- 
ciones  siguientes:1 


cos  x  +  qi  sen  x+p«  cos  2x-\-q^  sen  2x+p3  cos  3x-}-q3  sen  3x  + 
cos  2x  +  qi  sen  2x-srp?  cos  4x  +  q2  sen  4z+/>3  cos  6x  +  g3  sen  6x  + 
cos  Stf  +  ^i  sen  3x+p>  cos  6x  +  g-j  sen  6x+p3  cos  9x  +  g3  sen  9a;  + 


cos  (w  —  l)x  +  qi  sen  (n  —  I)z-fp2  cos  2(n—\)x-\-qi.  sen  2(w  —  I)cc+p3  cos  3(n—  l)x 

+  q3  sen  3(w  —  l)x  +  ................ 

Si  el  numero  de  los  coeficientes  que  pretendemos  determinar,  fuera  el  mismo  que 
el  de  nuestras  observaciones,  tendriamos  un  sistema  de  tantas  ecuaciones  como  incog- 
nitas, y  podriamos  resolverle  facilmente  por  los  metodos  ordinaries  del  algebra.  Pero 
ordinariamente  el  numero  de  las  observaciones  excede  con  mucho  al  numera  de  coefi- 
cientes que  necesitamos  ;  asi,  en  el  analisis  de  la  curva  barometrica,  teniamos  veinticuatro 
observaciones,  y  empleabamos  unicamente  ocho  de  estos  coeficientes. 

En  este  caso  el  sistema  es  mas  que  determinado;  y  en  la  imposibilidad  de  hallar 
los  valores  de  p  y  q  que,  sustituidos  en  las  ecuaciones,  reproduzcan  exactamente  nues- 
tras observaciones,  hemos  de  contentarnos  con  averiguar  cuales  son  los  valores  mas  pro- 
bables de  p  y  q,  que  dan  para  los  segundos  miembros  valores  compatibles  con  las 
observaciones. 

En  la  teoria  de  los  minimos  cuadrados  se  prueba,  que  estos  valores  mas  probables 
son  aquellos  que  hacen  minima  la  suma  de  los  cuadrados  de  los  errores  residuos. 

Error  residue  es  lo  que  se  aparta  del  valor  observado  el  valor  obtenido  al  sustituir 
en  las  ecuaciones  los  valores  mas  probables  de  las  incognitas. 

Segun  esta  definicion,  designando  por  VM  vlf  v2,  etc.,  los  errores  residuos,  y  supo- 
niendo  que  p,  plf  </,,  etc.,  son  los  valores  mas  probables  de  estos  coeficientes,  podemos 
escribir  : 


cos  x  +  qt  sen  x+p2  cos  2x  +  q2  sen  2x+ps  cos  3x-\-qs  sen 
cos  2x  +  qt  sen  2x+p^  cos  4*  +  g2  sen  4x+ps  cos6x  +  q3  sen 
cos  3x  +  q>  sen  3x+p~  cos  &x  +  q~  sen  6x+p3  cos  9x  +  q»  sen 


vn— i— an— i+p  +  pi  cos  (n  —  \)x -\-q-i  sen  (n  — I)x+p2  cos  2(n  —  I)x  +  q2  sen  2(n  —  l)x 

+ps  cos  3(n— I)x  +  q3  sen  3(w 

La  condicion  que  hace  mas  probables  los  valores  de  p  y  q  es 

vl+v'i+vl+vl  + +  v~n-i  =minimo. 

y  escribiendo  en  vez  de  v%,  v\,  etc.,  sus  valores,  queda 

(—art+p+p1+p.2  + )2+  (-^H-P+P,  cos  x+q^  sen  x+ )2 

-f  (— a2-hp-hPi  cos  2x+ql  sen  2#+ )2=minimo. 

Ahora  bien ;  sabemos  por  el  calculo  dif erencial  que,  cuando  una  funcion  esta  en  un 
minimo,  sus  derivadas  con  respecto  a  cada  una  de  las  variables  son  cero.     Por  consi- 

1  La  p,  que  aparece  en  estas  ecuaciones,  representa  la  constante  M  que  vimos  completaba  la  for- 
mula de  Bessel;  notese  ademas  que  las  ecuaciones  de  la  pagina  43  estan  dadas  aqui  segun  la  Formula  (I). 
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guiente,   derivando  la  expresion   anterior  con  respecto  a  p,  plf  qlf  etc.,  y  representando 
por  P,  Pj,  Qlf  etc.,  estas  derivadas,  obtenemos  el  sistema  siguiente: 

P  =  2(  —  af>Jrp-irp1-}-p2  +  ........  ...)+2(  —  ch+p  +  pi  cos  a;-|-gi  sen  x~\-  ............  ) 

+2  (—  .ck+p+jOi   cos  2x-\-q,  sen  2x-f  ............  )+  ............  =  0 

Pi  =  2(—a.o+p+p1+p-i  +  ............  )  +2  cos  x(  —  OI+P+PI  cos  x-  +  <?i  sen  x  +  ............  ) 

+2  cos  2x(—  ck+P+pi  cos  2x-\-qi  sen  2x-f-  ............  )+  ............  =  0 

Qi  =  2  sen  x(  —  o-j+p+pi  cos  X-\-QI  sen  #  +  ............  ) 

+  2  sen  2cc  (  —  a2-)-P+pi  cos  2x-\-q^  sen  2a5  +  ............  )  _|_  ............  =  0 

P.  =  2(—  a«  -\-p-\-  PI  +p^  +  ______  )+2cos  2x(—al-\-p-\-p1  cos  X-\-QI    sen  x-\-  ............  ) 

+2  cos  4x(—  dz+p+pi  cos  2£  +  q'i  sen  2x-\-  ............  )+  ............  =  0 

Q  =  2  sen  2x(—  ck  +  p+pi  cos  X  +  Q!  sen  »  +  ............  ) 

+2  sen  4x(  —  (k  +  p+pi  cos  2a;-f  tfi  sen  2x  +  ............  )+  ______  —  0. 


Este  es  el  sistema  que  hemos  de  resolver  para  hallar  los  valores  de  p  y  q,  que  seran' 
los  mas  probables;  pues  en  esta  hipotesis  se  ban  obtenido  las  ecuaciones  precedentes. 
Despues  de  dividir  por  2,  podemos  darle  la  siguiente  forma: 


IP—  —  o«     +  np  +  pi  1 

+  cos  x 
+  cos  2x 
+  cos  3x 
+  cos-4x 


cos(n-l)x 


-f  sen  2x 
+  sen  3x 
+  sen  4cc 


+  cos  2x 
+  cos  4or 
+  cos  6x 
+  cos  8» 


+  cos  2(n-l)x 


sen  2a- 
+  sen  4x 
+  sen  6x 
+  sen  80; 


+P 


ai  cos  x 

ttzCOS  2x 

a3  cos  3x 


—  an-i  cos  (n-1)  it- 


+  cos  x 
+  cos  2x 
+  cos  3x 


+  cos(n-l)a; 


-f-  cos2  a: 

+  COS'  2.T 

+  cos"  Bx 


cos  a;  sen  x 
+  cos  2x  sen  2x 
+  cos  3x  sen  3x 


+  cos  x  cos  2# 
+  cos  2x  cos  4x 
+  cos  Bx  cos  6x 


+  cos(n-l)x  sen( 

cos  x  sen  2a; 
+  cos  2x  sen  4a; 
+  cos  3*  sen  6x 


+  cos(n-l)£c  cos  2(n-l)x 


-ai sen  x 
-a?  sen  2x 
-as  sen  3x 


+P 


sen  x 
-f-  sen  2cc 
+  sen  Bx 


sen  x  cos  x 
+  sen  2x  cos  2* 
4-  sen  Bx  cos  3# 


+  sen(n-l)x 

sen  x  cos  2x 

-  sen  2x  cos  4x 

-  sen  Bx  cos  6x 


+  sen(«-l)x  cos  2(n-l)x 


+  €08(^-1)0;  sen  2(n-l)x 

-\-q\  I        sen2  a; 
+  sen2  2x 
+  sen2  Bx 


+  sen2  («-!)« 
sen  a;  sen  2x 
+  sen  2a;  sen  4x 
+  sen  3o;  sen  6x 


+  sen(w-l)a;  sen 


.  =  0 
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Este  sistema  segun  la  notacion  explicada  en  la  pagina  47,  puede  escribirse  del 
modo  siguiente: 

—2  am  +  np+pi  2  cos  mx  +  qi  2  sen  mx+p*  2  cos  2mx  +  q*  2  sen  2mx-\- =  0 

—2  am  cos  mx+p  2  cos  mx  +  pi  2  cos2  mx  +  qi  2  cos  mx  sen  mx+p*  2  cos  mx  cos  2mx 

+  92  2  cos  mx  sen  2mx-\- =0 

—2  am  sen  mx-\-p  2  sen  mx  +  pi  2  sen  mx  cos  mx-\-qi  2  sen2  mx+p2  2  sen  mx  cos  2mx 

+qi  2  sen  mx  sen  2mx+ =  0 

—2  am  cos  2mx+p  2  cos  2mx+pi  2  cos  2mx  cos  wx  +  gi  2  cos  2mx  sen  mx-\-p2  2  cos2  2m«  '     ' 

+  g2  2  cos  2mx  sen  2mx-\- 

—2  am  sen  2mx+p  2  sen  2mx+pi  2  sen  2mx  cos  ma;  +  gi2  sen  2?/u;  sen  mx-\-p2  2  sen  2mx  cos  2mx 

+  q2  2  sen2  2mx  + 

Cuando  wx  (producto  del  numero  de  observaciones  por  el  numero  de  grades  de  las 
partes  en  que  estan  divididos  los  distintos  circulos  de  referencia)  es  multiple  de  2?r,  el 
sistema  anterior  se  reduce  al  siguiente:  : 

—2  am  cos  mx-\-p1  2  cos2  mx=0 
—2  am  sen  mx-\-q^  2  sen2  mx=0   (^ 
—2  am  cos  2mx-{-p2  2  cos2  2mx—0 
— 2  am  sen  2mx+g9  2  sen2  2wx=0 

Pero  como  luego  probaremos 

2  cos2 

2  sen2 

V  (C) 
2  cos2  2mx=!£, 

M 

2  sen2 


luego,  sustituyendo  estos  valores  en  el  sistema   (B)   se  obtiene  finalmente 

,—  2  a» 

—2  am 

— 2  am  cos 


^^1=2  alu  cos 

~  gi=2  am  sen  mx 
2 

^P2=2  am  cos 
?  ^=2  am  sen 


Estas  son,  en  su  forma  general,  las  ecuaciones  de  Bessel  que  buscabamos ;  en  nuestro 
caso  particular  de  veinticuatro  observaciones,  ^=12,  lo  que  reduce  estas  ecuaciones  a 
las  que  escribimos  en  la  pagina  47.  Queda  probado  tambien  que  la  cantidad  cons- 

1  Hemos  juzgado  conveniente  dejar  para  el  fin  esta  demonstracion,  para  no  interrumpir  el 
raciocinio. 
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tante  p  o  M,  que  completa  la  formula  de  Bessel  es,  como  indicamos  antes,  la  media 
aritmetica  de  las  observaciones  ;  asi  en  nuestro  caso,  tenemos: 


24 

Faltanos  unicamente  probar  el  paso  del  sistema  (A)  al  sistema  (B)  y  las  igual- 
dades  (C). 

Para  lo  primero,  nos  basta  probar  que,  cuando  nx  es  multiple  de  2-n-,  2sen  mx=Q, 
2  cos  mx—Q,  2  sen  2mx=0,  2  cos  2mx=0,  etc. 

Tomemos  la  formula  bien  conocida  en  trigonometria  : 

2  sen  z  sen  y=cos  (z-y)  -cos  (z+y)  /.sen  z=cos  (*-^)-cos  OH-  1/)  ^\ 

2  sen  y 

For  una  serie  de  aplicaciones  de  la  formula  (1)  en  la  que  escribiremos  sucesiva- 
mente  en  vez  de  z,  x,  2x,  3x,  ............  ,  (n-l)x,  y  en  vez  de  y,~  obtenemos: 


sen  x= 


2sen 


sen  2x= 


2  sen 


sen  3x= 


2  sen? 


— — — 

sen  (n— 2)x= — — 


2Sen 


sen  (n— \)x=  — 
y  sumando  ordenadamente  queda : 

o 

sen  #-|-sen  2x-\-sen.  3x-\- -{-sen  (n—l)x=^senmx=— 


2  sen- 


^—  cos  nx  cos^—  sen  nx  sen- 


si  nx  es  multiple  de  2-rr,  sabemos  que 

cos  nx=l;  sen  nx=Q 
y  queda  finalmente  : 

cos|—  cos| 

2  sen  w#  =  -  =0 

2  sen? 
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For  analogo  procedimiento  se  probaria  que  2  sen  2mx=Q,  2  sen  3mx=0,  etc.,  asi  para 
este  ultimo,  escribiriamos  en  la  formula    (1)   Bx,  Qx,  9x,  ............  ,  (n—  1)3#,  en  vez  de 

O/£ 

z,  y     -  en  vez  de  y,  lo  que  daria  : 


cos^-cos  (Qn~~  ^x-    cos^-cos  3nx  cos^-sen 
2  sen  Zmx= 


2  sen^ 


y  como  3nx  es  multiple  de  2™,  resulta:  2  sen  3ra#=0 

Para  probar  que  2  cos  mx=Q,  2  cos  2ma;=0,  etc.,  nos  serviremos  de  la  formula  trigo- 
nometrica 


2  cos  z  sen  ?/=sen  («+2/)  -sen  (z-y).'.cos  z= 


sen 


i  sen  y 

/y* 

y  escribiendo  en  vez  de  z,  0,  x,  2x,  ............  ,  (n—  l)x,  y  en  vez  de  y,  „,  obtendremos  por 

£ 

analogo  procedimiento  : 


sen- -I- sen - 
cos0= — A_ 

2  sen? 


—  sen- 

2  2 

-  * 


2  sen? 


senf-senf 

cos  2x= =- 

2  sen? 


senT|-senf 
cos  3x= — 


cos  (n-2)x= 


X 


2sen| 
Sumando  ordenadamente  obtenemos : 

spn_ 


1+cos  a^+cos  2a;+cos  3ic+ +cos  (n— 1)#=2  cos  mo;=- 

2sen| 

sen^+sen  no?  cos^— cos  wa;  sen| 


x 
2 
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y  teniendo  en  cuenta  que  nx=2^,  queda: 

x  x 

sen-  —  sen  -= 

2  cos  mx=  —  —  0 

2VV 
sen  g 

Probado  lo  que  precede,  facilmente  se  ve,  que  los  demas  terminos  del  sistema   (A) 
de  la  forma 

2  cos  mx  sen  mx  ;  2  sen  mx  cos  2mx  ;  2  sen  mx  sen  2mx  ;  2  cos  wz  sen  2mx  ; 

2  cos  mx  cos 


se  reducen  tambien  a  cero;  pues  recordando  que  el  producto  del  seno  por  el  coseno  de 
un  mismo  arco  es  igual  a  la  mitad  del  seno  del  arco  duplo,  se  obtiene: 

2  cos  mx  sen  mx  =^2  sen  2mx=0 

lo   mismo   podemos   probar   para   las   cuatro  expresiones  restantes  valiendonos  de  las 
bien  conocidas  formulas  de  trigonometria  : 


sen  x  cos  #=|sen  (#+?/)  +isen  (x—y) 
cos  x  sen  7/=isen  (#-|-2/)  —  ^sen  (x—y) 
cos  ic  cos  y—^cos  (x-\-y)  -|-|cos  (x—y) 
sen  x  sen  i/=icos  («—  ?/)  —  ^cos  (x+y) 

pues  escribiendo  en  ellas  en  vez  de  x  e  y,  mx  y  2m#,  obtenemos  : 


sen  m#  cos  2mx—%$  sen  Smcc—  ^2  sen  mx 
cos  w#  sen  2mx=%2  sen  3mx+^2  sen  mo; 
cos  w#  cos  2mo?=^2  cos  3w#4-^2  cos  mo- 
sen  mx  sen  2m£=|2  cos  mx—  12  cos 


y  por  lo  probado  anteriormente,  los  segundos  miembros  son  cero. 

Pasemos  a  probar  las  ecuaciones  (C)  ;  un  facil  desarrollo  trigonometric©  hace  ver  que: 

2  cos2  #=l-fcos  2x 
segun  esto  tenemos: 

2  cos2  0=1+  cos  0 
2  cos2  o;=l+cos  2x 
2  cos2  2ic=l+cos  4# 
2  cos2  3ic=l+cos  Gx 


2  cos2  (n—  l)z=l+cos  2(n—l)x 
y  sumando  ordenadamente 

2  2  cos2  mx=n+l+cos  2x-\-  cos  4a;-fcos  6ic-f  ............  -fcos  2(n—  l)x=n+2  cos 

pero  segun  lo  anteriormente  demonstrado, 

2  cos  2m#=0 
luego 

2  cos2  w#=- 

A 

Por  analogo  procedimiento  probariamos  este  mismo  valor  para  un  termino  cualquiera 
de  la  forma  2  cos2 


Para  probar  que  2  sen2  mx=™,  nos  valdremos  de  la  conocida  igualdad. 

z 
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sen2  x=l— cos2  x 


y  escribiremos 

sen2  0=1— cos-  0 

sen2  x=l— cos2  x 

sen2  2x=l— cos2  2x 

sen-  3#=1— cos2  3x 


sen-  (n— l)x=l— cos2  (n— l)x 
sumando  ordenadamente  tendremos : 

2  sen2  mx=n— 2  cos2  mx 

pero  acabamos  de  probar  que2  cos-  mx=^,  luego: 

5  sen2  mx=n— 

No  juzgamos  oportuno  detenernos  mas  en  dar  a  conocer  todo  el  metodo  de  Bessel; 
quien  en  ello  estuviera  interesado,  podria  consultar  la  explication  original  del  mismo  Bes- 
sel, traducida  al  ingles  por  Mr.  R.  H.  Scott  en  un  apendice  al  "Quarterly  Weather  Report 
of  the  Meteorological  Office,  London,  1870." 
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